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ABSTRACT 


Title of Thesis: An X-Ray Survey of the Sky From Balloon Altitudes 

Gunter R. Riegler, Doctor of Philosophy, 1969 

Thesis directed by: Frank B. McDonald, Visiting Professor 


A series of balloon-borne X-ray experiments was carried out between 
December 1965 and July 1967. X-radiation in the energy band from 12 to 60 keV 
was detected with a combination of proportional and scintillation counters . 

A clear relation between the background counting rate at a residual at- 
mospheric depth of 3 g/crn and the geomagnetic latitude was found. The 
secondary atmospheric X-ray background is assumed to be a component of at- 
mospheric electromagnetic cascades which originate in the collisions of high- 
energy cosmic-ray particles with atmospheric nuclei; other X-ray production 
processes near the top of the atmosphere are shown to be insufficient to explain 
the observed secondary background. The correct prediction of the geomagnetic- 
latitude effect is shown to constitute a stringent test for models of atmospheric 
X-ray production. 


Two observations of the Crab Nebula (Tau XR-1) resulted in net source 
spectra which can be described by a power- law photon spectrum with a spectral 
index of n = 1. 93 Qg . A comparison with the results of measurements by 
other experimenters at lower and higher energies indicates that the X-ray spec- 
trum of Tau XR-1 is a power-law of index n ~ 2 . 0 at energies between 2 and 
400 keV. 
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A measurement of Tau XR-1 with good statistical accuracy was used to 
demonstrate two distinctly different methods of data analysis . Although the 
spectral parameters remained stable, the statistical uncertainty in the indices 



showed some dependence on the analysis method. 

Three observations were made of the source Cyg XR-1, which has been 
reported by other observers to have variable intensity. The measurement with 
the highest statistical significance yields a spectrum similar to that obtained 
from Tau XR-1 and a flux level that is consistent with the results of an observa- 
tion by another group of experimenters, conducted three weeks prior to the 
measurement reported here. 

An experiment on 15 December 1966 showed that the spectrum of Sco X-l 

between 15 and 35 keV is consistent with thin-source bremsstrahlung from a 

7 

plasma at a temperature of approximately 5 • 10 °K. However, the photon 

flux at energies between 35 and 60 keV was found to be distinctly above the ex- 
trapolation of the lower-energy spectrum and was significantly higher than that 
observed a few months later by other groups. Similar high-energy deviations 
from a simple exponential spectrum have been observed by other experimental 
groups on two occasions. The existence of this high-energy feature and the 
variability of the X-ray spectrum are discussed in terms of two models for 
Sco X~1 which have been proposed recently. 

An observation of the vicinity of the galactic center on 18 December 1966 
detected a significant flux of hard X-radiation. Subsequent experiments by other 
groups have confirmed this finding, but to date a clear identification of the source 
or sources of this hard X-radiation has not been given. 
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Three observations of the vicinity of the north galactic pole brought 
contradictory results. An increase in the counting rate during the first experi- 
ment on 5 December 1965 was interpreted as the X-ray flux from a source in 
the direction of the Coma cluster of galaxies, which is near the direction of 
the north galactic pole. Subsequent observations failed to detect a similar 
feature and set upper limits to the photon flux from sources in the vicinity of 


the north galactic pole. 
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CHAPTER I 


INTRODUCTION 

Within the past 6 years, X-rays have become an important new tool for the 
astronomer. X- and gamma-radiation have extended the observed energy range 
of electromagnetic radiation from discrete celestial objects to encompass more 
than 12 orders of magnitude. Observations of the diffuse cosmic photon back- 
ground have been carried out over a total energy interval of 15 orders of magni- 
tude. 

The experimental methods of X-ray astronomy are characteristic of both 
cosmic- ray measurements and optical and radio astronomy. On one hand, 
X-rays are treated like corpuscular particles which are counted and energy- 
analyzed; however, unlike charged particles, X-rays traverse magnetic fields 
without deflection so that the source position may be deduced from the direction 
of their arrival at the earth. 

X-ray astronomy is expected to furnish new information for our under- 
standing of high-energy processes in stellar atmospheres as well as in inter- 
stellar space. We will briefly review some processes which are involved in the 
production of photons in stellar atmospheres. The following mechanisms are of 
greatest interest: 

(a) atomic and molecular processes 

(b) Coulomb scattering 

(c) interaction between a charged particle and the magnetic field 


1 



(d) collisions between photons and electrons 

(e) nuclear reactions 

(a) Atomic and molecular processes result in the absorption and emission 
of both lines and the continuum. Line spectra, which have been studied most 
thoroughly at optical wavelengths, are generated in stellar atmospheres and 
provide information on ion abundances, temperature, velocity, and the atmo- 
spheric pressure. For example, the observed optical lines may be members of 
the Balmer series for hydrogen and, similarly, of differing series for other 
atoms, molecules, and ions. Recent refinements in radio engineering tech- 
niques have made it possible to observe recombination lines from atomic levels 
with very high total quantum numbers, e.g. , H 109a, H 158a, and He 158a. 

A hyperfine transition in the ground level of neutral hydrogen gives rise to the 
21 cm line. Radio observations of this line allow the determination of the spin 
temperature and the integral number of neutral hydrogen atoms along the line of 
sight. Emission lines in the low-energy X-ray region of the spectrum are 
expected to be generated in hot plasmas at temperatures of more than 10 °K. 
The most prominent X-ray lines from a Maxwellian gas with cosmic abundance 

7 

distribution at a temperature of approximately 10 °K are emitted at energies up 

+25 

to 7. 15 keV (line emission of Fe ). However, detection of such lines in the 
X-ray spectrum has not been successful to date. 

Edges in the continuum spectrum are caused by photoelectric absorption 
and by radiative recombination. They are produced in the source of electro- 
magnetic radiation or in the medium between the source and the observer. The 
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emission spectrum of a Maxwellian gas of the type described above has recom- 

+25 

bination edges at energies up to 9. X keV (the ionization energy of Fe ). 

(b) Coulomb scattering of one charge in the field of another charge with 
accompanying emission or absorption of photons is called bremsstrahlung or 
free-free emission, respectively. Although the terms are sometimes inter- 
changed, these processes are the inverse of each other and their cross-sections 
are related. Electron bremsstrahlung is an important emission mechanism for 
hot ionized gases. The shape of the spectrum depends on the optical depth, the 
temperature, the elemental abundance distribution, and the degree of ionization 
of the emitting plasma. Thermal bremsstrahlung emission has been detected 
from a variety of sources and over the full range of the observable photon 
spectrum, for example from interstellar HII regions at radio frequencies and 
from some sources of low-energy X-radiation. If the plasma is at a uniform 
temperature, then this temperature can be derived from the form of the observed 
spectrum. The amplitude can be related to the product of the electron and ion 
densities multiplied by the source volume. 

The collision of a non- relativistic proton with a low-energy electron can 
also give rise to bremsstrahlung at X-ray energies. A simple transformation 
converts this process to the more familiar case of electron scattering in the 
field of a stationary proton. If the electron is initially bound, then one speaks 
of "internal bremsstrahlung. " 

(c) Interaction between a charged particle and the magnetic field. Mag- 
netic bremsstrahlung is caused by the interaction of charged particles with 
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magnetic fields. At a given velocity, the emitted power is inversely proportional 
to the square of the mass of the particle so that emission by electrons is much 
more important than emission by protons. 

Synchrotron radiation is generated by relativistic particles and is assumed 
to be the source of non-thermal optical and radio emission from the Crab Nebula. 
This hypothesis is supported by the detection of polarization which is a necessary 
feature of synchrotron radiation from ordered magnetic fields. Non-relativistic 
magnetic bremsstrahlung from electrons, called cyclotron radiation, has been 
observed at radio frequencies from Jupiter and the sun. 

(d) Collisions between photons and electrons. Elastic collisions between 
photons and free electrons are called Compton collisions. When an energetic 
electron transfers some of its energy to a photon, then one speaks of the 
"inverse Compton effect." 

The single observation that is perhaps most intimately correlated with 
the study of galactic and extragalactic cosmic rays is an observation of the 
diffuse cosmic X- and gamma-ray background. Recent measurements have 
indicated that this background is- approximately isotropic. It has been suggested 
(Hoyle, 1965; Gould, 1965a; Felten and Morrison, 1966a) that this diffuse X-ray 
flux is generated by inverse Compton scattering of electrons by the universal 
3. °K black-body spectrum. If the Compton effect is assumed to dominate, the 
appropriate part of the electron spectrum can be deduced from the spectra of 
the low- and high-energy photon distributions. However, the real situation is 
considerably more complicated because other energy loss mechanisms influence 
the equilibrium electron spectrum. Felten and Morrison (1966a) have proposed 
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the existence of an extragalactic electron component in addition to the more 
intense galactic component. At present there is no consistent theoretical model 
which also agrees with the results of measurements on high-energy electrons. 

(e) Nuclear reactions generate gamma radiation either directly, or 
indirectly through positron and neutron formation. Among the more important 
reactions that produce line emission are the neutron capture by H , deuteron 
production in proton-proton collisions with subsequent pair annihilation, and 
tritium production. Of course, the positron in the pair-annihilation reaction 
above can also be supplied by many other reactions. 

Proton-proton collisions produce high-energy gamma rays via pair 
annihilation or tt ° -decay. The subsequent electromagnetic cascade involves 
pair-production and -annihilation, bremsstrahlung, and Compton scattering and 
will produce a substantial flux of lower-energy gamma radiation. 

The r-process of nuclear synthesis (Burbidge et al. , 1957) produces 
several heavy radioactive isotopes. The iowest-energy emission line with high 
intensity is generated at 60 keV by Am-241 (Jacobson* 1968). Decay via internal 
conversion gives rise to the emission of characteristic K-shell radiation from 
various heavy isotopes at energies between 100 and 130 keV. 

It is apparent that the observation of high-energy photons provides 
valuable information about the astrophysical processes and the energy balance 
in the source. For example, the total luminosity in the X-ray region is larger 
than that in the optical region by a factor of about 3 for a power-law energy 
spectrum of index 1. A spectrum of this form is assumed to Q rn the optical 
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and X-ray emission of the Crab Nebula. On the other hand, a thin-source 

7 

thermal spectrum at a temperature of 5*10 °K can explain some of the obser- 
vations of Sco X-l (sec Chodil et al. , 19GH). In this case the luminosity ratio 
Ixjcomcs approximately 10' » 

2 

The attenuation length in air for photons of 30 ReV energy is 3 g/cm , 

which makes it necessary to conduct measurements of cosmic X-radiation near 

the top of the atmosphere. At 3 kcV the attenuation length is reduced to 
2 

0. 005 g/cm* - so that cosmic X-rays of this energy cannot be measured at 
balloon altitudes. The first evidence for the arrival of X-radiation from beyond 
the solar system was obtained in June 19G2 (Giacconi et al. , 1962), with a set 
of rocket- borne uncollimatcd Geiger counters which were sensitive to X-rays 
from 2 to s A. The experimental objective was the detection of fluorescence 
radiation produced on the surface of the moon by solar X-radiation; instead, 
radiation from a direction near the galactic center was observed. The existence 
of this source in the constellation Scorpius was later verified (Gursky et al. , 
19G3; Bowycr et al. , 1964a)* In addition, these later flight records indicated 
f the arrival of X-rays from the constellations Cygnus and Taurus. 

hi July 19G4 a rocket measurement was performed at the moment when 
the moon occulted the central region of the Crab Nebula (Bowyer et al. , 19G4b). 
This experiment demonstrated that the Crab Nebula was the source of the 
X-radiation from Taurus. The angular width of the X-ray source was deter- 
mined to be about 1 minute of arc. The finite angular extent provided a strong 


argument against the neutron-star hypothesis for the source Tau XR-1. Another 
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resu/L >i these rocket experiments was the evidence for the oxistcnee.of a 
diffuse cosmic X-ray background. This cosmic background flux is roughly 
isotropic (Grader et al. , 1966), and has been measured at energies up to 
approximately 1 MeV (Metzger et al. 1964). 

The first X-ray astronomy experiments used small- area detectors with 
poor angular resolution. In subsequent observations the effective detector 
area was increased and the angular resolution and rocket pointing accuracy 
were improved. It became possible to resolve the X-radiation from the con- 
stellations Cygnus, Sagittarius, and Scorpius into discrete X-ray sources. 

Early in 1967, a list of 30 cosmic X-ray sources was published by Friedman 
et al. (1967b), but the typical uncertainty in the source position was still of the 
order of a few degrees. Independent surveys by other groups (Gursky et al. , 
1967; Fisher et al. , 1968; Bradt et al. , 1968a) showed only partial agreement 
with the source locations given by Friedman et al. (1967b). At present 22 
cosmic X-ray sources (Webber, 1968) have been confirmed by various rocket 
experiments. 

Very few of these 22 sources of X-radiation have been identified with 
celestial objects. Almost 2 years after the identification of Tau XR-1, an 
inge.-i^wS collimator made it possible to measure the location of Sco X-l to 
within 30 seconds of arc (Gursky et al. , 1966). The search for the correspond- 
ing optical object was guided by an extrapolation of the observed X-ray spectrum 
into the optical region. Sco X-l was then identified with a flickering blue object 
of magnitude V = 12.6 (Sandage et al. , 1966). Similar methods were applied to 
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sources in the constellation Cygnus. The X-ray source Cyg X-2 was tentatively 
identified with a starlike object (Giacconi et al. , 1967b), which was later shown 
to be a binary system (Kraft and Demoulin, 1967). Cyg XR-1 could not be 
identified with an optical object because of the heavy obscuration in the direction 
of that particular region of the galaxy (Giacconi et al. , 1967b). 

Low-energy X-ray measurements from rockets are restricted by the 
absorption of X-radiation in interstellar hydrogen. The lower energy limit for 
rocket observations of nearby galactic X-ray sources is at approximately 
0. 2 keV. The upper limit is imposed by the short exposure time in rocket 
experiments since the statistical significance of the observed high-energy 
photon flux is reduced for a steep photon spectrum. Consequently, measure- 
ments from rockets are usually carried out in the energy interval from 
approximately 1 to 20 keV. 

The attenuation of X-radiation in air restricts experiments at balloon 
altitudes to energies above approximately 15 keV. I u addition to the attenuated 
cosmic diffuse background, there is a diffuse background of atmospheric 
X-radiation at all altitudes. This background and the attenuation of X-radiation 
in the residual atmosphere limit the photon flux from cosmic sources that can 
be measured at a given altitude. Practical detectors also have a certain 
counting-rate background that is internally generated and further restricts the 
smallest photon flux that can be measured with a given field of view. The great 
advantage of balloon observations lies in the possibility of long exposure times 
to cosmic X-ray sources. Exposure times of a few hours are feasible, in 

• ' ■ • ■ ' lV . '• ■ ' 

contrast to rocket experiments which have durations of a few minutes. 
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Balloon observations of cosmic X-ray sources began in 1964, 2 years 
after the first rocket measurements. Clark (1965) observed X-radiation from 
Tau XK-l up to energies of 62 keV. This result furnished new evidence against 
the hypothesis that X-radiation from the Crab Nebula is produced by black-body 
emission from the surface of a neutron star. Subsequent experiments from 
balloon altitudes, for example by Peterson et al. (1968), resulted in refined 
measurements of the emission spectra of Tau XR-1 and of other X-ray sources. 
Long exposure times made it possible to extend the observed spectrum up to 
energies of a few hundred keV. Together with the results of rocket measure- 
ments, the X-ray spectrum for certain sources has been determined over more 
than two decades in energy. 

The present series of experiments sought to study cosmic X-radiation 
at balloon altitudes. Detectors were designed for photon measurements at 
energies between 10 and approximately 100 keV. The two-fold purpose of the 
survey was to find sources of hard X-radiation and to measure their spectra. 
This aim dictated the choice of the angular resolution and the number of 
channels for the pulse-height analysis. 

Various sectors of the northern and southern sky were scanned in a 
series of balloon flights. Some regions were selected because they were 
already known to contain sources of hard X-radiation. The vicinity of the 
galactic center was studied because the detection of soft X-rays from several 
sources suggested that a flux of high-energy radiation might be observable. 

In the absence of discrete sources, anisotrophy of the diffuse radiation in the 
direction of the galactic center seemed to be physically possible. 
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Some of the more prominent production mechanisms for X-radiation will 
be reviewed in the following chapter. Chapter 3, "Detectors, " contains a 
description of the proportional and scintillation detectors as well as the 
properties of the shield and the methods of data handling. In chapter IV, 

" Summary of Experiments, " we discuss relevant information for all balloon 
experiments and the detector orientation systems which were used in the 
observations. The numerical and statistical handling of the data is described 
in the chapter on the " Methods of Data Analysis." The chapter entitled 
"Results" contains details of the various X-ray source observations, the 
results of the final data analysis and a discussion of each source or source 
region. 

A review of topics related to the production of X-radiation at balloon 
altitudes is given in appendix B. 


spa 








CHAPTER II 


X-RAY PRODUCTION MECHANISMS 


Several reviews of cosmic X-ray production mechanisms have been 
written (Morrison, 1958; Hayakawa et al. , 1964a; Hayakawa and Matsuoka, 1964b; 
Ginzburg and Syrovatskii, 1965; Oda, 1965; Hayakawa et al. , 1966; Gould and 
Burbidge, 1967b; and Gould, 1967a). The basic processes are well understood 
and can be grouped into two classes: electron-photon processes (synchrotron 
radiation, bremsstrahlung, and Compton scattering) and nuclear and atomic 
processes (gamma-ray line emission by radioactive or fissionable isotopes 
produced in r-process synthesis, internal conversion radiation from excited 
nuclei following radioactive decay, atomic recombination lines, etc.). 

Line emission will not be discussed here because the energy resolution of 
the present series of experiments is insufficient to detect any of the lines that 
might be expected. Also, we do not expect emission lines of high intensity 
within the energy region that we are considering (sections III B and III C) . For 
example, lines from hot plasmas with universal abundance occur at energies 
below 10 keV (Tucker, 1967). The most important lines from the decay of 
radioactive isotopes which were synthesized in the r-process have photon 
energies above 60 keV (Jacobson, 1968). 

The identification of a source spectrum as an exponential or power-law 
function of energy does not unambiguously define the source mechanism. Just as 
in the case of the optical continuum spectra, which are not well described by a 

11 
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simple black-body function but, rather, by solutions of the radiative transfer 
equation, an interpretation of X-ray spectra is not straightforward. Sartori 
and Morrison (1967a) have proposed a model where a non-thermal radio source 
emits thermal bremsstrahlung from plasmas at two different temperatures. 
Klimas (1968, private communication) has studied the emission spectra of 
plasmas with given distributions of temperature and electron density. The type 
of spectrum emitted from such sources was found to depend strongly on the 
spatial distribution of temperature and density. Manley (1966) has shown that 
the synchrotron - emission spectrum of electrons with a flat energy distribution 
is nearly exponential. 

With these considerations in mind, we will proceed to review the basic 
X-ray production mechanisms. 

A. Synchrotron Radiation 

The synchrotron process appears to be the principal mechanism for 

photon production in some radio sources. For example, the polarization of the 

radio and optical radiation from the Crab Nebula can best be explained as due to 

2 

the synchrotron mechanism. An electron of energy E = Ymc , moving in a 

e 

magnetic field H emits a synchrotron emission spectrum which has a maximum at 

—8 o —90 O 

E = 1. 16 • 10 H, Y" = 4. 47 • 10 H E , (1) 

m i. j_ e 

where the electron and photon energies are measured in eV, H is measured in 

Gauss, and = H cos 6 is the strength of the magnetic field perpendicular to 

the electron orbit (Hayakawa et al. , 1964a). 
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The synchrotron spectrum is usually approximated by a delta function 

j: 

I 

spectrum at E^. The power radiated by an electron of energy is 

P <E, E 0 ) = 0.98 ■ 10" 3 H 3 (E e / me 2 ) 2 6 (E - EJ , (2) L 

L 

where P is in units of eV/sec. If the electron spectrum can be described by a j 

il 

power-law of the form 

I I 

N(E ) dE =KE“ r dE , (3) ! j 

' e' e e e * ' ' j-j 

then the observed synchrotron spectrum in the case of a homogeneous magnetic 

r 

field becomes il 

l 

Ig (E) “ II P(E, E e ) N(E e ) dE e dn = A-K^r) e'^ 1 !^ , (4) 

; t 

where A is a constant and K is a function of the spectral index r. The emitted 

1 : | 

photon spectrum is a power-law in energy with an index that is related to the jl 

1 

! | 

index of the electron spectrum. j 1 

hi the presence of a random magnetic field, equation 4 is still applicable, l| 

ii 

j 'J 

provided K (r) is replaced by a slightly different function of T, and H. is re- || 

placed by |h| (Ginzburg and Syrovatskii, 1964). 

For an electron which loses its energy E only by the synchrotron pro- 
cess, we can define a lifetime T such that the rate of energy loss is equal to 

E/T. The lifetime is then equal to 

o mc 2 

T = 0. 513 • 10 J“J [2 (sec) (5) 

e l 

13 

(Hayakawa et al. , 1964a). For example, an electron of energy E = 2. 7 • 10 

v 

-4 

eV will originally produce radiation at E^ = 30 keV in a field of = 3 • 10 

8 

Gauss. Such an electron has a lifetime of 10 sec 3 years. 
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A magnetic field of the order of 10 Gauss is assumed to be present in 
the Crab Nebula. If high-energy electrons have been generated in the original 
supernova explosion with a spectrum which produces the observed radio emis- 
sion, then they cannot generate the measured X-ray spectrum because of the 
high rate of synchrotron loss (Tucker, 1967). If the X-ray flux is produced by 
synchrotron emission, then it must be assumed that continuous production of 
high-energy electrons takes place in the Crab Nebula. The detection of polari- 
zation of X-ray emission would support the synchrotron hypothesis. 


B. Inverse Compton Scattering 


2 

We consider an electron of total energy E =7mc , in the rest frame S 

0 

of the observer (Felten and Morrison, 1966a). The electron collides with a 
photon of energy E which is incident at an angle ol. It is easiest to perform the 
calculations in the rest frame S’ of the electron, where the outgoing photon angle 
is Qi’ and the electron scattering angle is 0 T . In this rest frame the interaction 
is just the standard Compton scattering. After transformation back to the frame 

I 

S, one obtains the energy E^ of the scattered photon (Felten and Morrison, 1966a] 

Y 2 E (1 + cos a) (1 - /? cos a') 

i “ 1 + (yE/xhc2) (1 + /? cos a) (1 - cox 6 1 ) ' ' 

2 

As long as the photon energy is small, 7E«m c , in the rest frame of the 
electron, equation 6 can be approximated by 

E 1 » i Y 2 E . (7) 

In this case classical Thomson scattering prevails in the electron rest frame. 
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The power scattered by one electron moving through a photon distribution of 


3 V 


energy density p (in eV/cm ) is 
P 0 (V, p) = - Thorason ° P 


E 


-14 2 

= 2.10 7 9 , 


( 8 ) 


where P is measured in eV/sec. For an isotropic distribution of photons, 
c 

equation 8 can be integrated over all angles; this results in an additional factor 
of 4/3 for the scattered power. 

As in the discussion of synchrotron radiation, it is customary" to approxi- 
mate the spectral power a? /dE (E, 7, p) by a delta function at a characteristi 

c 

energy E 


dP 

_c.(E, 7, p)«P (7, p) <5(E-E ) . 
dE c c 


( 9 ) 


The average energy of a black-body distribution of photons is 

<E> = a K T , (10) 

where a « 2.7. The characteristic energy of the emitted photons becomes 

E « 4"7 2 a kT (11) 

c 3 

(Ginzburg and Syrovatskii, 1964). With a power-law energy spectrum for the 


electron, 


-r 


N(E ) dE = K E * dE 
' e e o e 


e 


(12) 


we obtain the power scattered per unit energy in the electron source volume in 

the presence of an ambient black-body distribution of photons at temperature T: 

I (E) = (E, E , p) N(E ) dE = P ° . (yE> P } , (13) 

c' ' o dE e’ ' e' e 2*£. Y„ a kT 

3 E 


where 7„ is given by (equation 11) 
E 


7 = 
E - 


E 


c 


4 ofkT 


(14) 
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We finally have a power-law photon spectrum with the same index as the 
synchrotron spectrum 


.4 /3 V 2 =T ££ -TH, 

I c (E) = 10 p N 0 (^) 2 T 2 E 2 (15) 

The Compton-synchrotron process. A magnetic field and high-energy 

electrons in the presence of photons having, for example, a black- body 

distribution, give rise to both inverse Compton and synchrotron emission. The 

peak in the synchrotron photon spectrum, E is related to that of the Compton- 

m 

process photon spectrum, E . When equations 1 and 11 are solved for V, 

c 

we obtain 


ft kT 

E = 1.15 . 10 a E , 
c H m ’ 


(16) 


where E , E , and kT are measured in units of eV, H is in Gauss, and ol- 2. 7, 
c m * 

For most astrophysical applications , the ratio E /E is very large compared 

o m 

to one. If a power-law spectrum of electrons is assumed as before (equation 3) 
then the ratio of the emitted intensities for the synchrotron and Compton 


processes is 


T s 3 H 2 /l, 8 . 10 4 T\ 
2 8 frp \ H J 


3-r 

2 


(17) 


where p is again the photon energy density, T is in °K, and H is in Gauss. 

Felten and, Morrison (1963) have suggested that a power-law energy 
distribution of electrons in the galactic halo produces both the synchrotron- 
radiation background observed at radio wavelengths and the X- and gamma- ray 
backgrouhd via the Compton process. However, detailed calculations (Felten 
and Morrison, 1966a) show that the predicted Compton spectrum in the presence 


of a cosmic 3 °K black-body distribution of photons (Dicke et al. , 1965; Penzias 
and Wilson, 1966) falls short of explaining the observed cosmic X-ray background 
spectrum by about two orders of magnitude. The spectral shape of the predicted 
X-radiation background is compatible with the observed one. An extragalactic 
electron component would be required to match the observed X- and gamma-ray 
spectrum. 

Similar models have been forwarded for discrete sources. Gould (1965b) 
attempted to explain the Crab Nebula X-radiation by a Compton spectrum caused 
by the electron distribution that produced the observed radio spectrum via the 
synchrotron process. A discrepancy of approximately five orders of magnitude 
between the observed and the calculated X-ray intensities was obtained. 


C^ Bremsstrahlung 


We assume that electrons of density n have a Maxwellian energy distri- 

6 

bution around an energy kT and interact with ions of charge Z and density n.. 
The equation for the bremsstrahlung spectrum as given by Alien (1963) can be 
modified to give the number of photons emitted per unit time, volume, and 
frequency: 


[tT e^Z 2 fm ,-hi\ 1 

Y? ~ykF Vt e(hu ’ kT) exp <Tt> hF ■ 


. 6 2 

dN .16 / 7r e Z 

— — = 47T-- 

dv 3 


(18) 


where dN/dj; is measured in units of photons/cm -sec-Hz, and the Gaunt 


factor g(hv, kT) will be defined below. In an exact treatment one must fold a 


Maxwellian electron distribution with the bremsstrahlung formula in the 


nonrelativistic Born approximation including the Sommerfeld factor. The 
Sommerfeld factor can be set equal to 1 if 
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27tZ _ , , 2rrZ _ „ 

137(3 1 311(1 137/3' <<: 1 ’ (19) 

where (3 and /3 f are equal to v/c for the initial and final electrons, respectively. 

This is a reasonable approximation for electron energies above 100 keV and for 


Z = 1. 


A comparison of the final result (Wharton, 1961) with equation»18 yields 


the Gaunt factor for any value of hv/kT: 

g<h„, kT) =*§ K q (gj) exp (gj) , 


where K (x) is the zero-order modified Bessel function of the second kind. With 
o 

the asymptotic expansions of K (Abramowitz and Stegun, 1964) , we obtain 

|f [ In f^-) -0. 577j exp (~j for hr « kT, 


g(hr, kT) = 


k 2kT 
hv 


= 0.977 


for h^ » kT. 


Equation 21 Shows that the limiting value of the Gaunt factor for Z = 1 and for 

large photon energies is not equal to one, but is approximately equal to VkT/hv. 

The observed spectra of several X-ray sources can be interpreted as due 

to bremsstrahlung of electrons with a Maxwellian distribution at a temperature 
7 8 

between 10 and 10 °K. If the electron gas contains a wide range of kinetic 

■m 

temperatures, then the main contribution to the overall spectrum of the brems- 
strahlung curve for a temperature T comes from energies around kT. Both 
arguments show that the full Gaunt factor must be used in the formula for the 


bremsstrahlung spectrum. 
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For photon energies above 10 keV whieh are comparable to or greater 
than kT, the contribution from electron-electron bremsstrahlung becomes sig- 
nificant. Electron-ion (Z = 1) bremsstrahlung is dipole in nature while 

electron-electron bremsstrahlung is quadrupole in nature; one expects their ratios 
2 

to be me /kT. Maxon and Corman (1967) have calculated electron-electron 
bremsstrahlung from a Maxwellian hydrogen plasma in the non- relativistic Born 
approximation. The ratio of quadrupole to dipole emission at a particular photon 
energy, B ( hv/kT ), can be approximated by 


B (lr) = ( 2 + 0 - 64 It) S* ’ for bv Z kT - 10 keV. 


( 22 ) 


For example, in a hot Z = 1 gas at kT = 10 keV, the electron-electron contribution 
is 7% of the electron-ion emission at a photon energy of hv - 25 keV. The dif- 
ferential bremsstrahlung spectrum for Z = 1 is therefore best represented in the 
energy region considered here by a function of the form 

P (hv, kT) = A exp (-^)k 0 (^)- 

[ 1+ ( 2+0 ‘ 64 s)^ 5-] d • (23) 

This equation will be used later in the data analysis. 
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CHAPTER III 


THE DETECTORS 
A . Introduction 

The X-ray sensing equipment consists of three elements: a gas propor- 
tional counter, an alkali-halide scintillation crystal, and a combination of active 
and passive shields which also acts as a collimator. X-radiation from a well- 
defined aperture is detected by this counter arrangement, so that we may speak 
of it as an X-ray telescope. 

Figure 1 shows a schematic diagram of the X-ray telescope. A krypton 
gas proportional counter is located in front of a CsI(Tl) scintillation crystal, 
and both are located inside a plastic scintillator cup. The light output pulses 
of the Csl crystal and the plastic guard counter are distinguished by their 
rise times. The proportional counter is also connected in anticoincidence 
with the scintillators. The plastic cup is the innermost part of a three-element 
graded shield which consists of tin on the outside and copper in the middle. 

B. The Proportional Counter 

2 

The X-ray transmission of a layer of air with d = px = 3 g/cm exceeds 
0. 01 at 15 keV and exceeds 0. 1 at 20 keV. Krypton gas at 1 atm. pressure was 
selected for a proportional counter of 5. 08 cm depth so that the absorption of 
this gas exceeds 0. 5 at all photon energies between the K absorption edge at 

2 

14. 323 keV and approximately 24 keV. The absorption curve of 0. 0188 g/cm 
krypton as a function of photon energy is drawn in Figure 2 as a solid line. 


20 


The broken line shows the product of this absorption function times the trans- 
mission functions of a film of aluminized mylar (0. 00254 cm polyester, 0. 00381 
cm Al), of the beryllium plate at the top of the pressurized detector housing 
(0. 3175 cm Be), and of the beryllium entrance window of the proportional 
counter (0. 0127 cm Be). 

The proportional counter has the shape of a right cylinder, its inner 
dimensions being 5. 08 cm diam. x 5. 08 cm length. The original version had a 
button-type anode; the supporting grids for the entrance and exit beryllium 
windows (0. 0127 cm Be) and the side wails were made of Monel. Later counters 
had a wire electrode, and the wall material was aluminum. The full-width at 
half maximum (FWHM) of the response to a 25 keV line (main peak of the 
response to Sn-119 m, see below) was 20% for the earlier counters and 32. 5% 
for the later models. 

The energies of emission lines from some radioactive sources which 
were usei in the energy calibration of the proportional and scintillation 
counters are listed in Table I. The relative line intensities have been adopted 
from the measurements of Holt (1967a). The "effective energy" is the peak 
energy for the superposition of several Gaussian response curves with in- 
tensities and mean energies as given in the table. The widths of the Gaussian 
curves were adjusted to correspond to the energy resolution of the proportional 
counter and the scintillation crystal, respectively. Also shown in Table I is 
the line broadening due to the composition of each calibration line. 

A typical energy calibration curve of the proportional counter is shown in 
Figure 3. The solid line connects the effective line energies of radioactive 


sources plotted as a function of the channel number. The curve is not linear, 
but it never deviates more than 1 keV from the linear approximation (broken 
line). This nonlinearity is caused by the detector itself; independent tests with 
calibrated pulses of the same shape as the signals from the proportional counter 
showed that the pulse-height analyzer was linear. The highest energy calibration 
point shown in Figure 3 (31. 5 keV, channel 54) would be expected to fall into 
channel 59, but the upper cutoff at channel 56 distorts the line response so that 
it appears to be peaked at a lower channel number. 

With the known energy calibration, one can convert the observed response 
to radioactive calibration sources into an energy spectrum. This is done in 
Figure 4 for Sn-119 m, using the nonlinear (solid line) and the linear (broken 
line) energy calibration curves. The half-width energy intervals at half 
maximum for the main and the escape peak are 4. 3 and 3. 63 keV, respectively. 

It is apparent from Figure 4 that the escape peak is not symmetrical; this 
fact is only partially explained by the influence of the lower threshold. As will 
be discussed in more detail later, escaping characteristic krypton K radiation 
can be absorbed in the Csl crystal. The K-line energy (12. 6 keV) is below the 
threshold of the crystal, but a voltage equivalent to two channels is added to the 
combined proportional counter/ scintillator output which therefore shifts the 
simultaneous escape signal for the proportional counter. According to the 
krypton counter calibration, this channel increment corresponds to approxi- 
mately 1 keV. Including the line broadening, we finally arrive at a relative 
resolution of 32, 5% FWHM at 25 keV and 54% FWHM at 12 keV. The full-width 
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energy resolution at half maximum can be represented by 

AE = >/l9 + 1- 9E , (24) 

where energies are measured in keV. 

In an ideal proportional counter, a finite resolution is due to: (a) the 
number of ion pairs formed by an incident photon, and (b) the subsequent multi- 
plication process, both of which are subject to statistics on finite numbers. 

4 

The multiplication factor of our counters is of the order of 10 for 2000 *V 
operating potential at the proportional counter anode. Depending on the choice 
of the Fano function (Fano, 1947; Curran and Wilson, 1965; Aitken, 1968), 
one arrives at a theo/retical resolution of typically 8% at 25 keV. The main 
contribution is caused by fluctuations in the number of ion pairs formed so that 
the variance is proportional to the energy. In other words, the full width at 
half maximum is proportional to the square root of the energy. 

The strong energy-independent contribution in addition to the linear term 
under the square root in equation 24, indicates that the second step (b), above, 
is responsible for the large energy dispersion. Since the electrode and 
detector configurations are unfavorable, charge amplification, charge col- 
lection, and dead regions introduce a constant term into the variance of the 
energy dispersion. 

When a photon of energy E is absorbed in a given material, fluorescence 
X-radiation fromK, L, and other shells can be emitted by the material, pro- 
vided that E is above the absorption edge for these shells (Appendix A). The 
escape factor, v (E), is defined as the probability of emission of a fluorescence 
X-ray for every photon that has been absorbed in the detector. There is no 
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simple analytical way of calculating r? for the proportional counter because the 
calculation would have to take full account of the detector configuration. We can 
estimate the value of the escape factor in the following manner. The K-shell 
fluorescence yield for krypton is /? K = 0. 66 (Bergstrdm and Nordling, 1965; 
Fink et al. , 1966). The probability of creating a K-shell vacancy is taken to be 
(Grodstein, 1957) 


a K 


lv3 


= 0.86 . 


1 + (* 2 “) + 


(25) 


The product a p must be multiplied by the average attenuation of the K-shell 
K Jv 

fluorescence photons. If the average path length of these photons is assumed to 
be about 1/3 of the characteristic detector length, s = 5. 08 cm, one arrives 
at an escape factor 


ri ~ 0.48. 


(26) 


The escape factor was also determined experimentally by a comparison of 
the source-response curves with calculated response curves. The energy reso- 
lution as a function of energy must be accurately known in order to obtain 
reasonable precision in the determination of V . Observations of the primary 
and escape peaks in the response to Sn-119 m in different proportional counters 
resulted in escape factors of 


V = 0. 50 - 0. 52. (27) 

C. The CsI(Tl) Scintillation Crystal 
When the present series of experiments began, two alkali halide scintil- 
lators were in general use: Nal(Tl) and CsI(Tl). CsI(Tl) has an efficiency for 
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the conversion of absorbed energy into light equal to approximately 30% of that 
of Nal(Tl) (Aitken, 1968), and a rise time of 1. 8 }j.sec compared with approxi- 
mately 1 }j.sec for Nal(Tl). CsI(Tl) is therefore more useful for the Phoswitch 

i 

technique of pulse-shape discrimination; in addition, the mechanical properties 
of the Csl crystal are more favorable, allowing easier handling of the crystal. 

The scintillation crystal has the shape of a right cylinder, 0. 2 cm in height 
and 5. 08 cm in diameter. The counter is optically thick in the energy region 
from 20 to approximately 100 keV. 

The "effective" line energies for groups of closely spaced lines from 
radioactive sources have been computed with a reasonable estimate for the 
scintillator resolution. They are listed in Table I under the heading "Effective 
Energy, Csl. " A typical energy calibration curve is drawn in Figure 5. The 
lowest-energy points (Sn-119 m at 24. 7 keV and Cd-109 at 22. 5 keV) do not 
fall on the calibration curve since the lower cutoff distorts the source-response 
curve and effectively shifts the observed peaks. 

With the known energy calibration, we can measure the energy resolution 
of the scintillation counter with all processing electronics in operation. The 
half- width at half maximum for the 86. 5 keV line of Eu-155 is 14. 5 keV; the 
upper haif of the response is distorted by the high-energy cutoff. This value 
corresponds to a relative FWHM of 33. 5%. The effective line energy of 
Ba-133 is 31. 5 keV, a spreading of 0. 6 keV is caused by the line composition 
(Table I). At this energy the resolution is 5 keV HWHM, corresponding to a 


value of 29% FWHM. 
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The Phoswitch rejection efficiency was reduced at energies above approxi- 
mately 60 keV. We have therefore restricted the data analysis to energies 
below 60 keV. 

The calculation of the escape factor rj is discussed in appendix A. A 
graphical representation of rj as a function of energy is given in Figure 6. 

D. The Electronic Circuitry 

In this section we summarize the detector pulse-selection criteria and 
the method of data recording (Doong and Noordzy, 1966). For simplicity we 
introduce the following labels for pulses from the three detectors: "A" events 
occur in the proportional counter; ”B" events are pulses originating in the 
Csl scintillation crystal; and "C" events originate in the plastic scintillator 
cup. A block diagram of the pulse logical analysis is given in Figure 7. 

1. Proportional counter events (A counts). When an X-ray photon 
interacts in the krypton proportional counter, a charge Q, proportional to the 
amount of energy deposited in the detector, is produced. This charge is con- 
verted inter a voltage signal by a charge-sensitive amplifier. If this signal 
corresponds to an energy between 10 and 30 keV, then it is analyzed and 
recorded as a mode 1 event, provided all other requirements (Table II) are met 
and the system is not busy. Whether the system is busy processing some other 
event or not, the event is stored as an A count in the prescaled A accumulator. 
This bit-rate signal is read out every 5 seconds (mode 2, see Table IE). 

2. Csl scintillation crystal events (B counts) . Signals from either or 
both the Csl crystal and the plastic scintillator can be received in the 
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photomultiplier. The procedure by which a signal from the plastic scintillator 
is identified as a C event will be discussed below. At present we assume that 
B signals have been extracted. They are first amplified and shaped. A valid 
B signal must fulfill a pulse-height criterion equivalent to 20 keV s E £ 100 keV. 
The subsequent analysis and recording as a mode 0 event proceeds like that for 

* 

an A event, provided the other requirements (Table II) are met and the system is 
not busy. All B events are stored in a prescaled accumulator. The B bit rate 
is read out together with the A bit rate every 5 sec (mode 2). 

During the 100 msec period immediately preceding the mode 3 command 
pulse, B events with E >100 keV are accepted and stored in a B overload 
accumulator. This information is read out during the mode 3 readout sequence. 

3. Plastic scintillator events (C counts) . C events are distinguished 
from B events by their shorter rise time. The signal from the preamplifier 
is simultaneously applied to a 40 nsec shaper with a threshold detector (equiva- 
lent to E > 10 keV) and to a peaking circuit. Plastic events of typically 30 nsec 
rise-time and E >10 keV will cause the two outputs to coincide and are 
accepted as valid C pulses. A stop gate is then activated which prevents the 
preamplifier signal from being processed in the pulse-height analyzer and which 
blocks the B bit-rate accumulator. 

Csl events with a rise-time between 0. 5 and 1. 8 fisec will be analyzed 
in both the B and the C systems. These B pulses are analyzed by four C 
threshold detectors which are set at 10, 20, 40, and 100 keV. An indication of 
the number of C threshold detectors that have fired is recorded along with the 
pulse-height analysis information for every B event (mode 0). B events with 
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simultaneous indications that one or more C threshold detectors have fired 
were not used in the subsequent data analysis. Valid C events with E > 100 keV 
are counted during a 100 msec interval preceding a mode 3 command pulse and 
recorded in the C overload accumulator. Both the B and C overload rates are 
read out every minute (mode 3). 

4. Mode 2. As previously mentioned, the A and B bit rates are recorded 
every 5 seconds. The rate was first scaled by a factor of 12 (Experiments 
1-5), and later by a factor of 8 (Experiments 6 - 11). In the latter case, a 
maximum recording capacity of 64 x 8 = 512 counts per 5 sec in each of the 

A and B channels was available. The mode 2 information serves as a check on 
the deadtime correction. 

5. Mode 3. Mode 3 occurs every 12 mode-2 readouts, i. e. every 
minute. For a duration of 100 msec, the number of B(E > 100 keV, 0 with 
E > 10 keV) and the number of C (E > 100 keV) events are recorded. 

6. Mode 7. Mode 7 is a time mark which is recorded on tape every 12 
occurrences of mode 3 (every 12 minutes). The count of the number of times 
mode 7 has occurred since the system was last turned on is recorded in 64 
channels. After 64 x 12 = 768 minutes, this count recycles, 

7. Modes 4, 5, and 6. For detector telescopes which were rotated in 

azimuth angle and programmed in elevation angle, magnetometer pulses were 

recorded on tape as housekeeping information. The data recorded in mode 4, 

5, and 6 are identical io those in mode 0, 1, and 2, respectively. The only 

difference is an additional mark which is added in the third of the three mode- 

0 12 

symbol tracks on the recording tape with binary values of 2 , 2 , and 2 . 
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A signal from the magnetometer will cause addition of 2 = 4 to the mode symool 
of the next data record, for example mode 0, 1, or 2. The addition is maintained 
until the next occurrence of mode 2, after which the addition is inhibited for at 
least 10 sec. 

E. Interaction Between the Proportional Counter and the Scintillation Crystal 
Because of the proximity of the proportional counter and the Csl crystal, 
and the criteria for event analysis (Table n, section HE D), some interaction 
between the counters is possible. For example, (a) K fluorescence radiation 
of the proportional counter can be absorbed in the crystal, and (b) K fluorescence 
radiation of the Csl crystal can be absorbed in the proportional counter. 

(a) The K fluorescence radiation (section m B) of krypton can escape 
from the proportional counter and may be absorbed in the scintillator. There 
the event is too low in energy to be accepted as a valid Csl event; therefore, 
the anti-coincidence requirements are not fulfilled. It is apparent from section 
in D and from Figure 7 that a voltage equivalent to 12. 6 keV is added to the 
voltage output for the escape peak in the proportional counter. Therefore, the 
escape peak is shifted in energy by an amount which equals two channels in the 
system used here, i. e. , approximately 1 keV equivalent photon energy in the 
proportional counter. 

In order to obtain an estimate of the relative probability e of this inter- 
action, we multiply the escape factor of the proportional counter by the proba- 
bility of escape into the Csl crystal. The latter will be of the order of the solid 
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angle of the Csl crystal seen from the center of the proportional counter, 
divided by 47 r: 

(28) 


r\ (E) • -l~ s i 

^K, Kr W 4 t r 


The escape peak will be shifted by an amount of energy AE « 1 keV, where AE 
is subject to resolution spreading. For E = 25 keV, the probability e is 


e =0.52 * i =0.13 . (29) 

(b) The K fluorescence radiation of Csl can escape from the entrance 
surface (appendix A) and may be absorbed in the proportional counter. If the 
primary photon energy was between 35 and 50 keV, then the energy deposited in 
the crystal is between 5 and 20 keV, and again the pulse will not be accepted as 
a valid Csl crystal event. The 30 keV escape photon can be accepted as a valid 
event in the proportional counter even though it stems from a primary photon with 
an energy between 35 and 50 keV. 

In order to estimate the magnitude of this effect, one must consider the 

response of both detectors to a common source. We assume a power-law 

spectrum incident at the top of the atmosphere. Then we compute the attenuation 
2 

in 3 g/cm of air as well as in the beryllium covers, etc. , and determine the 

response to the remaining spectrum in the proportional and scintillation counters. 

The total flux of escaping photons for escape through the entrance surface of the 

crystal and for primary energies between 35 and 50 keV is given by 
50 keV 

Escape = / ” < E > ’reap < E >' dE * < 30 > 

35 keV 


iiPjapawggBMBg^ggBBS 
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I (E) is the attenuated source photon flux which is absorbed in the Csl 
resp w 

crystal. As the next step we determine the absorption probability and the count- 
ing rate distribution in the proportional counter for the line at E - 30 keV. For 
exposure of both detectors to a power-law source spectrum of index 2, we find 
that the contribution of escape from the Csl crystal into the proportional counter 
amounts to 6. 2 percent at 17. 5 keV, 2. 5 percent at 24 keV, and 8. 8 percent at 
30 keV. This small contribution is shifted by an amount equivalent to 
0 SE $1 keV and is again subject to resolution spreading. 
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F. Geometric Factors 

Figure 1 shows that the Sn - Cu - Foam - CH shield arrangement 
also acts as the collimator for the proportional and scintillation counters. The 
geometric factor for the scintillation crystal is defined as 


GF = <Afi> = 27r 


r 

| A o “ 

•Ti 


(0) sin0 d0 , 


where A (0) is the detector area for a parallel beam of photons incident on the 

detector without modulation. The cut-off angle 0 is defined by the collimator 

o 

geometry. The factor a (9) accounts for other modulation effects, notably the 
ring and disk placed in front of the X-ray telescope. 

Figure 8 shows the effective exposure area A (0) of the scintillation 
crystal as a function of the angle 0 from the detector axis. This quantity was 
found experimentally by placing a point X-ray source at a large distance com- 
pared to the detector and collimator dimensions and at various angles 0 . The 


\ 
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observed counting rate was multiplied by a normalization factor to give the 
calculated value .at zero angle 

A (g * 0°) = Trr^ = 20 cm 2 . (32) 

The solid line in Figure 8 represents the effective exposure area for the un- 
obstructed detector. It is very close to the calculated curve for the same 
detector geometry. The broken line shows A (0) for the detector with modu- 
lation ring and disk. 

A characteristic of the response curve is the angle at which the effective 
area drops to half its maximum value. The FWHM angle is given in Table in. 
The geometric factor is then simply calculated by evaluating the integral 

GF - 2 7r Ja (0) sin0 d0 . (33) 

o 

The curves in Figure 8 were obtained wi th an emitter of X-rays at an energy 
of E « 25 keV (Sn-119 m). The effective exposure area for high photon energies 
is larger than that for lower energies since the collimator is less effective in 
the former case., Edge effects and an increased contribution from Compton 
scattering off the inside walls increase the effective exposure area. A summary 
of geometric factors at 25 keV and at approximately 75 keV (Bi-109) is pre- 
sented in Table III for the proportional and scintillation counters with and without 
modulation. The value for the scintillation crystal response at 7 5 keV in the 
presence of the modulation assembly has been inferred from the other three 
distributions for this detector. The apprqpriate exposure factor is therefore 
shown m parentheses. 
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The definition of the geometric factor for the proportional counter is 
analogous to that of the scintillation crystal, but now the detector is not thin 
compared to the collimator length. One has to consider several thin layers 
each with a different absorption probability for incoming X-rays, and then 
evaluate the geometric factor of each layer of krypton gas. The variation in 
the charge-collection efficiency at different points in the counter also plays 
a role. 

The effective area as a function of angle Q from the detector axis has 
been evaluated experimentally. The resulting curves are similar to those for 
the scintillation detector. The full-width angle at half-maximum effective area 
is 27.6°. The geometric factors for the unmodulated and modulated counters 
are given in Table M. 
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G. Response of the X-Ray Telescope to Gamma Radiation 

The telescope shield consists of one active layer and several passive 

layers. The active part of the shield is a plastic scintillator cup (Pilot B, 

equivalent chemical formula CH ) which is a part of a graded shield con- 

X* £ 

sisting of tin on the outside, copper, Ecco Foam, plastic, and aluminum on 
the inside. 


This shield loses some of its effectiveness at higher energies; it is 
therefore necessary to investigate the quality of the shield for electromagnetic 
radiation of energy within the range of interest/ *1$ - TOO keV) and above. A 


suitable parameter is the transfer function s i 4 ' (S, E 


E ) for the Csl crystal. 

6s 


It is defined as the effective geometric factor of the Csl crystal (with 
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dimensions cm ster) for photons of energy E >. E winch deposit an energy E’ 
in the crystal, provided that E’ lies between E and E . In other words, we 

X ^ 

consider the X-ray telescope immersed in an isotropic flux of photons of energy 
E greater than E and define the transfer function TE (E, E , E ) as the ef- 

u 1 u 

fective geometry factor for photon measurements at energies between E^ and 

E . The quantity consists of two components: a contribution from photons E 

which enter through the aperture, TF^, and one due to photons which first 

traverse the shield before interacting in the crystal, TF : 

s 


TF ap (E ’ E l’ E 2> = If F (a, E, E x , E 2 )dfi, 


Aperture 


TF g (E, E r E 2 ) = f 

^Shield 


F (a, E, E ± , E 2 ) 


exp [- (£x. p., (E)/ sin aj* 27rd (cos 6!) . 


The definition of F (ci!, E, E , E ) is similar to the definition of TF (E, E , E ), 

1 Z 1 z 

except that it refers to photons which are incident at an angle a from the 

detector axis. The depth of the i-th layer of the shield is x^, and (E) is the 

total absorption coefficient of this layer material. 

We will now consider the experimental evaluation of the transfer function 

TF. For this purpose we cannot make use of equations 34 and 35 in their 

present form. We define 17 (E, E , E ) as the probability that a photon of 

ap 1 2 

energy E, which interacts in the crystal after entering through the aperture, 

deposits an amount of energy E' between E and E 0 . The definition of 77 is 

1 2 s 

similar: a photon of energy E, that would interact in the crystal in the 
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absence of a shield, must traverse the shield and deposit energy between 
and E in the crystal. 

V ap and rj s can be found experimentally. We first determine the total 
counting rate at all energies for the bare crystal caused by a source of photons 
of energy E. Then, with the same relative source-detector configuration plus 
the shield in the case of T 7 g , we measure the counting rate at energies between 
E^ and E^. The ratio of the latter to the former rate is r? a p and rj g in 
absence or presence of the shield, respectively. The transfer functions are 
given by 


TF ap < E > E i’ E ol ; = (E, E i , E„) <Fn> 


1 2 J 'ap 


V 2' 


ap 


[l “ exp (- fo.(E) X q ) 


(36) 


and TF g (E, E^, E^ = rj g (E, E^ E g ) <FS2> g 1 1 - exp (-fJ-(E) <x> g ) ] (37) 

2 

where x^ is the crystal thickness and < Ffi > gp =3 cm ster (section in F). For 

an isotropic photon distribution we have 

1 


< Ff2 > = 


i F - 

) 

-1 


and < x > = 


cos a • 2 7 r d cos a= 2ir F »<FI2> 

o s 


x 


(38) 


A cos a. — 2 7 T d cos a 
cos a. 


x 

f 


2 x « <x> 
o s 


(39) 


A cos ol * 2tt d cos oi 


Experiments verified that 77 (E, E , E ) is a weak function in a, and rj was 

S X a S 

therefore assumed to be constant. 

Figure 9 shows the total transfer function TF 


TF (E, E 2 ) = TF gp (E, E^ E g ) + TF g (E, E g ) , 


(40) 
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for = 20 keV and E^ = 60 keV. Four values of the energy E were used in the 
experiments: 122 keV (Co-57), 357 keV (Ba-133), 622 keV (Cs-137), and 1.3 
MeV (1. 17 and 1. 33 MeV of Co-60). 

We can use the transfer function TF and the geometric factor GF 
(previous section) to estimate the counting rate due to the atmospheric back- 
ground flux at balloon altitudes. We assume the atmospheric X- and gamma- 

ray flux to be isotropic and use the spectrum measured by Peterson et al. 

2 

(1967a) at d = 3. 5 g/cm and at a local vertical cutoff rigidity R = 5. 1 GV. 

The transfer function and the geometric factor as a function of energy are then 

folded with the atmospheric photon spectrum, and the resulting counting rate 

2 

is adjusted to a residual atmospheric depth of d = 3. 0 g/cm . The counting 
rate at this depth caused by the atmospheric and detector- generated background 
is shown in Figure 10 (section III J), For the Csl crystal of an X-ray telescope 
without modulation ring and disk, we find the following counting rates: 

1. 8 counts/sec is the total counting rate caused by the atmospheric and 
the detector-generated background (the contribution from the cosmic 
X-ray background, approximately 0. 8 counts/sec, is not included). 

Of this rate: 

1. 1 counts/sec are caused by the atmospheric background photons 
with energies between 22 and 60 keV, and 
0. 5 counts/ sec are caused by the atmospheric background gamma 
rays with energies above 60 keV. Therefore, the remainder, 

0.2 counts/sec, are generated in the telescope by other effects, e.g. , 

by charged particles because of incomplete rejection by the Phoswitch. 
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It is difficult to estimate the possible error in each of the numbers quoted 
above. Therefore, a large error may be attached to the last-mentioned counting 
rate which cannot be accounted for. Nevertheless, these numbers show that a 
reduction of the solid angle of the field of view would lower the counting rate 
produced by atmospheric and cosmic X-rays and would therefore reduce the total 
counting rate. These considerations also show that it is not useful to reduce the 
aperture angle below a certain value without simultaneously reducing the, count- 
ing rate background due to gamma rays and other effects. 

H. The Atmospheric Transition Curve 
Figure 11 shows typical atmospheric transition profiles for the. scintil- 
lation and proportional counters. The Csl crystal data in the 22 - 60 keV 
interval are shown as double bars, and the krypton proportional counter data 
between 12 and 28 keV are represented by single bars. The measurements 
were taken by detector E during Experiment 10 (chapter IV). 

In order to find the contribution of the attenuated cosmic diffuse X-ray 
background to the counting rate at ceiling altitude, it is necessary to extrapolate 
the transition curve to that altitude. Such an extrapolation predicts the dependence 
of the sum of the detector- generated and the true atmospheric backgrounds. 

We will study the atmospheric transition profile for altitudes above those 
corresponding to the maximum in the counting rate. This maximum is the 
equivalent to the "Pfotzer maximum" and will be referred to as the "peak" in the 
transition curve. Several analytical functions have been employed to approxi- 
mate this section of the atmospheric transition curve, for example 
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I (p) « a - b exp (-p/c) (41) 

(Anderson, 1961; Brini et al. , 1965a), 

j 

] 

I (p) = a + bp (42) 

j 

(Peterson, 1967b), and 

I (p) = a p b (43) 

(Bleeker et al. , 1968). 

The sum of the primary and secondary cosmic-ray intensities as a 
function of pressure can be described in lowest-order approximation by a 

F 

function of the type f 

1 

£ 

X (p) -- a exp (-p/b) - c exp (-p/d), (44) j 

(Ghielmetti et al. , 1964b), 

2 i 

where b represents the absorption length, of the order of 160 g/cm , and d is j 

2 ■ 

the mean interaction length, of the order of 86 g/cm . The observed atmos- 
pheric transition curves for X-radiation can be fitted by functions of the same 
type. By means of computer calculations, least-square fits for the parameters j 

a, b, c, and d were found for the transition curves of the various data sets. 

,2 ! 

For example, the (6. 7 - 68. 0) g/cm section of the proportional counter data ' I 

of Experiment 10, detector E, can be best described by j 

m t 

I (p) a 3976 exp (-p/143) - 3617 exp (-p/75) . (45) j 

The appropriate curve is given in Figure 11 as a solid line. j 

I 

1 

It is apparent that the calculated lines present a good fit to the atmos- | 

pherie transition profiles, but it is not suggested that the observed counting rate 

i 

f 

is simply due to the primary and secondary atmospheric cosmic ray intensity. j 

ft 

... /.‘Vi . 'I 

A brief discussion will illuminate this point. , 1 
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(a) Atmospheric X-radiation at high altitudes is assumed to be primarily 
caused by the soft component of secondary cosmic radiation in the atmosphere. 
Some contributions to the atmospheric production of X-radiation and their 
relation to the cosmic-ray flux will be discussed in appendix B. 

(b) A study of the Phoswitch properties showed that the rejection threshold 
shifts at high counting rates, so that an increased number of plastic pulses are 
accepted as valid Csl crystal events. If we define a peak-to-ceiling ratio (PCR) 

as the ratio of counting rates at the peak of the atmospheric transition curve, 

2 2 
e. g. , at 90 g/cm , to that extrapolated to ceiling altitude, e. g. , at 3 g/cm , 

then the effect just described will increase the PCR for the scintillation crystal,. 

The PCR for the proportional counter will, of course, not be affected. 

(c) Observations have shown that the PCR increases between 20 and approx- 
imately 180 keV, and then gradually decreases again up to a few MeV (Peterson 

et al. , 1966a; Peterson et al. , 1967a; Peterson, 1967b). If this trend can be 
extrapolated to lower energies, a lower PCR will be expected for the proportional 
counter data than for the scintillator data. 

(d) Some Csl crystal events are due to higher-energy photons that have 
undergone Compton scattering (see section III G). The considerations of (c) 
apply in this case. 

We can now extrapolate the atmospheric transition curves to ceiling altitude 
and then use the calculated contribution of the diffuse cosmic X-ray background 
to infer the expected counting rate at ceiling altitude. Considering the effects 
discussed above, no great emphasis should be placed on the extrapolation or on 
the results derived from this procedure. 
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Table IV lists the extrapolated and the observed counting rates at ceiling 

2 

altitude (in counts/sec-cm -ster) for various detectors and several experiments. 
For each case, Table IV also shows the calculated contribution from the diffuse 
cosmic X-ray background. According to a summary given by Gould (1967), an 
average of several observations was used to define a two-component photon 
spectrum of the diffuse cosmic X-ray background: 

J (E) = 10 E -1, 5 for E < 23 keV 

(46) 

-2 22 

and J (E) = 100 E for E > 23 keV , 

2 

where J is the photon flux in units of photons/cm -sec-ster-keV and E is 
measured in keV. Appropriate adjustments for absorption above the detector 
as well as the detector efficiency and the geometry factor have been included. 

The counting rate at ceiling altitude is approximately 10% higher than 
the rate predicted from the extrapolation plus the computed contribution from 
the cosmic diffuse X-ray background. This difference can be understood as 
caused by the shape of the atmospheric transition curve. The "peak" section 
is too high compared with the higher- altitude section, thereby producing an 
extrapolated ceiling rate which is too low. 

J. Latitude Dependence of the Detector-Generated and the 
Atmospheric X-ray Backgrounds 

The following procedure was carried out for several experiments: The 
contribution of the diffuse cosmic X-ray background was subtracted from the 
counting rate at ceiling which was observed when no discrete cosmic X-ray 
sources were within the field of view. The resulting net counting rate was 

_ -^]LH jW Jj r i' lx i fic ii tr i jE ifntnniTOf ^i^'wfa ir 1 ^' T 
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adjusted to a common atmospheric depth of 3. 0 g/cm“' and was divided by the 
geometry factor GF (section III F) . Figure 10 shows the net adjusted event rates 
(in counts/cm 2 ~sec~ster) for the scintillation crystal, 

r = (R obs " R d> 1 GF = a + D / GF - (47) 


caused by the atmospheric (a) and the detector background (D/GF), as a function 
of the local vertical cutoff rigidity (in GF). Data from several experiments 
(Experiments 4, 9, 10 and 11) and several detectors (A, E, G and W) have been 
used. The points 4 A, 10 G, 4 A, 11 G, and 11 G define a curve of decreasing 
rate r with increasing local cutoff rigidity. The points 4 W, 9 E, 10 E, 4 W, 

11 E, and 11 E show more scattering around another curve, not shown in Figure 
10, which is roughly equidistant from the first one. The detectors E and W were 
partially occulted by a modulation disk and ring, while A and G were not so 
occulted. The geometry factors GF (section III F) were applied to the net event 
rate r * GF = R ^ - R^, which is derived from the true atmospheric (a • GF) 
and detector-produced events (D) (a and D can be assumed to be independent of 
the geometric factor of the detector). The net rate r is therefore larger for 
the occulted detectors which have a smaller geometry factor. This explains 
the shift between the two sets of data points shown in Figure 10 and the resulting 
background- versus-rigidity curves . 

The corresponding curve for the proportional counter could not be 
derived because some of the data points contradict each other. In general, the 
curve of proportional counter event rate versus rigidity is expected to reflect 
only the change of the atmospheric X-ray background with rigidity. 
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The shape of the curve of background versus rigidity shown in Figure 10 
caii be expected for any observation which is related to the charged-particle 
flux in the earth’s magnetic field. A comparison of the latitude dependence of 
the X-ray flux at ceiling altitude with the latitude dependences of other particle 
fluxes might give a clue as to the production mechanisms for the observed flux. 
Clearly, a much wider range in latitude would be required to make any definitive 
statements. Some aspects of this problem are discussed in appendix B. 
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CHAPTER IV 


SUMMARY OF THE BALLOON FLIGHTS 
A. Summary of Flight Data 

Between December 1965 and July 1967, 11 balloon-borne X-ray 
experiments were conducted. The flights have been labelled in chronological 
order. 

Table V lists the relevant information for eight of the flights. Launch data 
and the times for beginning and end of ceiling coverage are given, as well as 
the range of residual atmospheric depths and the range in local vertical cutoff 
rigidities that were encountered during the main portion of each experiment. 

The K indices are also presented to indicate geomagnetic activity during the 

P 

period of ceiling coverage. The names of the detectors that were flown in each 
experiment together with the scanning method that was used for each detector 
are also listed. The following rows in Table V summarize information on 
celestial objects that are relevant to each experiment. Solar flare information 
is given in the first row: the universal time of the maximum phase and the 
importance of each flare are quoted (ESSA, 1966, 1967). Finally, the names 
of known X-ray sources or of prominent directions in the celestial sphere are 
stated along with the times of their meridian transits. Ail abbreviations used 
are explained at the bottom of Table V. 
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B. The Balloon Gondola 

Three scanning mechanisms were used in the present series of experi- 
ments: the vertical, the oriented, and ; the rotating schemes (Table V). 

A vertical detector experiment, with or without modulation assembly (see 
section in F), performs a simple meridian-transit observation of an X-ray 
source. A necessary restriction is that the source must be within a band of 
width « 25° centered at the declination of the detector axis. A misalignment 
of the detector axis or of the modulation assembly is averaged out for free- 
rotating gondolas or for gondolas rotating with constant azimuthal velocity. 

Optical observations during Experiment 1 showed that a free balloon 
gondola rotates with changing direction and angular velocities of between 
0. 06 and 0. 2 rpm. The upper limit to the pendulum motion of a free balloon 
gondola was found to be approximately 0. 5°. 

The oriented detection scheme changes the effective latitude of the 
detector by tilting it to a certain fixed zenith angle. A band of the sky can be 
scanned by the meridian transit method; this band is then centered at the 
declination given by the new direction of the detector axis. Of course, the 
azimuth angle must be kept constant. We have aligned the detector axis in the 
plane of the meridian and have used the local direction of the magnetic field as 
a reference direction. Upon a deviation from the correct alignment with respect 
to the magnetic field, a magnetometer generated an error signal which controlled 
the release of an appropriate amount of gas through a jet nozzle and the resulting 
torque brought the gondola back to its proper azimuthal attitude. The correct 
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alignment of the detector axis and the modulation assembly is very important 
in this detection scheme. 

The rotating scanning method uses the magnetic field as a reference to 
adjust the otherwise continuous rate of change of the azimuth angle. The gondola 
rotates through 360° over a period of 10 minutes; a reaction wheel serves as a 
reservoir of angular momentum. A signal is recorded on magnetic tape 
(section in D) whenever the magnetometer senses the direction of magnetic 
east or magnetic west; another pulse is generated when the detector points 
toward north. These three marks provide the directional information necessary 
for the data analysis. 

In addition to the rotational motion, the zenith angle of the detector can 
be programmed so that the axis scans through a celestial point or tangentially 
along a celestial curve once during every rotation. Thus the projection of the 
detector axis executes a spiralling motion on the celestial sphere. Inaccurate 
alignment of the detector axis is not critical with this scanning' method. 


CHAPTER V 


METHODS OF DATA ANALYSIS 

Two different procedures have been used in the analy sis of the data and 
the extraction of the most probable source spectra. Procedure 1 was used when 
a channel-by-channel analysis of the data was warranted. Procedure 2 was 
used for all observations of cosmic X-ray sources. 

A. Data Analysis Procedure 1 

In each of the present series of measurements, an observed spectrum is 
divided into a number of discrete energy intervals of finite width, corresponding 
to the 64‘ channels of the pulse-height analyzer. The number of counts ob- 
served in each channel versus the energy- equivalent of the channel number may 
be represented by a vector V^. The source function can be approximated by a 

g 

similar vector V . Various media and devices act as operators on the source 
function. In our representation we can describe these operators by matrices M. 

The matrix which describes the influence of a layer of air of thickness 
d above the detector is simply given by 

< r = exp [- W E i> ' d ] < 48) 

and Mf * r =0 if i * j , 
ij 

where p . is the scattering coefficient of air without coherent scattering 
air 

(Grodstein, 1957), Similar matrices describe the attenuation in the various 
layers of aluminized mylar, beryllium, etc. , and in the krypton proportional 
counter gas in the case of the Csl crystal data analysis. All matrices are sym- 
metric and the off-diagonal elements are equal to zero. 
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The matrix M EsC describes the escape effect in the proportional counter 
and in the scintillation crystal. M Esc has the same shape in each case, but the 
numerical values are different. The probability of escape is assumed to be 
described by 17 (E) (sections HI B and III C and appendix A). For simplicity we 
combine all the characteristic line energies into one line energy E^, i. e. , 

- 12. 6 keV for the proportional counter and E = 30. 0 keV for the Csl 

rv 

crystal. If the subscript i denotes the rows of a matrix and j denotes the columns, 
then M Esc is described by 


..Esc * . 

M.. = 1-?)(E.), 

TTcri 

M.. = f) (E.) if E. = (E. + E, ) , 

ij v 3 .11 k 


(49) 


and M ESC - 0 if E. * (E. + E.) . 


13 


3 


The resolution matrix is not symmetric, but each row is symmetric 

with respect to the element of this row which lies along the diagonal of the 

matrix. The resolution is assumed to follow a Gaussian distribution with 

Res 

standard deviation o(E), The matrix M is defined as 


M* 68 

ij 


-W ex * ■%) • 


(50) 


The detector efficiency is described by a diagonal matrix M 
Eff 

M = 1 - exp [ - 0 . 0188 /v (E.)] 

11 Kr 1 

Eff 

and M. . =0 if i ^ 3 , 

13 

for the proportional counter , and by 


Eff 


(51) 
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..Eff . 

M.. - e (E.) 

n 1 


Eff 

and M, . 

iJ 


0 if i # j , 


for the scintillation counter. Tests on some scintillation detectors showed that 

the efficiency was equal to e = 0. 7 and was independent of energy. 

The relation between the source vector and that for the expected counting- 
E 

rate spectral function V is 

E Eff , T Res ..Esc ...Kt. _ Be ,Air 
V = M M M (M ) K . . . . M • V = ^ 

_ .total __S 
M * V . 

g 

For an assumed source spectral function V , we compute the expected counting- 
rate distribution as a function of the channel number i: 

f v fl I vE < e i>1 


v E = v E 


= V E (E 2> 


L V n J [_V E (E n )J 

This spectrum is now compared with the experimentally observed spectrum 
which is the distribution of observed counting rate versus the channel number i: 

'v?l [vV ' 

v° = V° = V°(E 2 ) (55) 

V° V°(E ) 

n L a 

J 

We might also take a different approach and invert the matrices , 
a , etc. , in order to find an expected source vector. A numerical 
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application of this method was investigated with the help of fast, large-core 
computers, but the results were disappointing, and the method was dropped. 
The main reason for this failure seems to be that a larger interval in energy is 
required for the matrix-inversion analysis. In this study, only 0. 5 decade of 
energy for the scintillation crystal and 0. 37 decade for the proportional counter 
were used. 

E O 

In order to compare the sets of numbers V and V , a meaningful 
quantity must be introduced as a figure of merit. Following Pearson’s sug- 
gestion (1900), this quantity is defined as the sum of the squares of the differ- 
ences between the observed set and the expected set of numbers, weighted by 

the inverse variance of the expected set of numbers, 

r O El 2 
V. - V. 

ft - 2 — - ( 56 ) 

1 <V E > 
i 

In our experiments , the observed spectrum V is always found by sub- 

B B+S 

tracting a background spectrum S from another spectrum S , which is as- 
sumed to be the sum of the background spectrum plus that of the X-ray source. 
The variance of the net observed counting rate in a particular channel i is 
vO o B u. e B+S 

<V. >-S. + S. . (57) 

ill 

E 

The ’’expected” variance < VV >in equation 56 must be derived by a similar 

procedure in order to be applicable in our data analysis. It is therefore a 

E O 

reasonable approximation to replace <V > in 56 by < V >, giving 


f = 2. 


ro E "1 2 

L v i - v i J . 


<v°> 
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As a next step we determine the probability P that, upon repetition of the 

}j 

i? 

experiment, a larger value of f would be encountered. For this purpose one lj 

h 

II 

h 

uses a probability distribution function for the number of degrees of freedom N jj 


which exists in the problem considered. In particular, for N ~ 30 the "chi- 


square distribution" is used: 



t 

e£ 

N 

2^ 



_ . N . 

r <T> 


(59) 


Pearson’s chi-square test then consists of evaluating a normalized least-squares 
sum according to equation 58 and of computing the probability Q of exceeding 
this value upon repetition of the experiment; 


Q(f|N) = 


°p dP 
f dt 


(N) dt, 


(60) 


where N is again the number of degrees of freedom. This function Q has been 
extensively tabulated (Abramowitz and Stegun, 1964). An approximation is 
available as a subroutine in the IBM system/360 Scientific Subroutine Package. 

In the practical analysis a given set of data was compared with 

E E 

several "expected" data sets V . Each set V , for a given spectral function and 
a particular set of two parameters, was obtained in the manner described above. 
One parameter, A, described the amplitude; the other, n, stood for the spec- 
tral index. The index n could be the exponent in a power- law spectrum or the 
value of kT for a thermal spectrum. Both parameters were varied until the 
corresponding value of f (eq. 58) reached an absolute minimum f^. The cor- 
responding probability Q therefore has the maximum possible value for the 
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given data set and the functional form fitted to it. The appropriate parameters 
will be called A Q and n Q . From the computation of f for varying A and n we 
obtain two functions, 

A ~ A(f,n) 

( 61 ) 

and n - n(f,A) . 

We must now define the statistical errors of the ,r best" parameters A^ 
and & . The mean value of a variable t which follows the chi-square distribu- 
tion (59) is t = N. The variance for this distribution is equal to 2N (Abramo- 
witz and Stegun, 1964, see below) . We adopt 

o- * n /Sn (62) 

as the "standard deviation" of a chi-square distribution with N degrees of 

freedom. The equivalent standard deviation in the parameters A and n is then 

the variation in these parameters that is generated when f- is varied by 

cr = v/2N. This is in contrast with the definition of Gorenstein et al. (1968) in 

which fg is varied by 1, independent of the number of degrees of freedom. 

In general, the variation will not be symmetric so that two values, e.g. , a 

and <T n2 , must be stated. For example, for a fixed A = A Q , we have 

n 0 =n(f 0 , A q ), (63) 

and the two values of <r are obtained from 

n 

n 0 +<T n = n ( f o +<T » V ' ( 64 ) 

The function n(f, A rt ) is double-valued for f>f_ and single-valued for f ~ f n . 

The analogous equations apply to the calculation of for a fixed n ~nQ. 

This procedure for the determination of a A and cr is lent support by a 

A n 

comparison with the more familiar situation in which a variable follows a 
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normal distribution. If the variable of the chi-square dis tribution is /it 
rather than t, then the mean of the distribution is (Abramowitz and Stegun, 
1964) 

l i+ i 6 -ife)3 + ° < n " 7/2 > - < 65 > 

and the variance becomes 

1 -dr + 0(tr3 >- < 66 > 

The following table summarizes the mean and variance of the chi-square 
distributions for two variables. 


Variable 

— ■ 

Mean 

Variance 

t 

v/iT 

N- 1 

N 

“ 11 
L 1 + 16N (N-1)J 

+ 

2N 

, 1 1 j. 

1_ 4N'iP 


For N £ 30 the chi-square distribution as a function oiy/2~t resembles a normal 
distribution of a variable s = yjl T, with unit variance and a mean s m =/2N -1 . 
In particular, for large N we can approximate the integral Q(fc|N) under the 
chi-square distribution by the integral under a unit normal distribution of the 
variable 

s' ~>/2% - v/2N - 1 . (67) 

Indeed, one-half of the distance in t-space between the points which correspond 
to s’ = +1 and s’ = -1 is equal to */2 N - 1 , in close agreement with the value 
or = y/2 N , used above. 

The reproductive law of the chi-square distribution (Alexander, 1961) 
states that if x and y are independent values from the chi-square distributes 
for M and N degrees of freedom, then x + y is distributed according to the 
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chi-square distribution for M + N degrees of freedom. This theorem can be 
used to combine independent sets of f~ values for certain parameter combinations , 
For example, the f- values for the proportional and scintillation counter expo- 
sure to a source were added for identical parameter sets and compared with 
the chi-square distribution for the total number of degrees of freedom. 

B. Data Analysis Procedure 2 

In this section we describe the method of data analysis which was used 
for all cosmic X-ray source observations. 

The observed counting-rate distribution as a function of channel number 
was grouped into a small number of bins j; for example, 2 bins for the pro- 
portional counter and 3 or 4 bins for the scintillation crystal. Each bin com- 
bined m channels. With a reasonable first assumption of the source spectrum 
S F 

V , we computed W (equation 54) by the procedure described in section A. 

E E 

Then a reduced vector V was found by averaging V . over all m channels 
which make up the larger bin j. Thus the extrapolation factors: 


F. = “g— , (68) 

J V. 

J 

account for the influence of the atmosphere and other absorbing materials and 
the physical processes within the detector. The factors F. were used to con- 
vert the observed counting rate in the bin j to a photon-flux spectrum at the top 
of the atmosphere. With a method to be described below, a best fit to these 
"extrapolated" data points was determined. As has been pointed out by other 


observers (Stein and Lewin, 1967), the extrapolation factors F. are relatively 
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insensitive to the detailed form of the input spectrum because of the approxi- 


mations made in their calculation. 


The result of this procedure was a set of photon-flux values y'(E^) with 

their statistical uncertainties <5y’(E.). The flux values are averages over energy 

3 

intervals AE^, centered at E^. Although the width of the energy intervals was 
not small in this analysis, its effect on the best-fit spectrum has nevertheless 
been ignored. 

A transformation was applied to the data points in order to be able to fit 

a straight line to the transformed points [ (y. + <5y.), x. ]. For a power-law 

3 3 3 

spectrum; 


y. = In [y'(E.) ] , 

x. = In (E.) , 

3 3 

ay, 

and 5y j = wt) = 

J 


111 

y'(E.) 


where 5y. denotes the statistical error of the transformed photon flux y.. 
3 3 

a thin- source thermal spectrum; 

y j = 1x1 [E j * y ’ (E j ) ] ’ 

x. = E. , 

3 3 


j y'(E.) 

The parameters a and b for the best straight line y = a + bx were calculated by 


means of a weighted least-squares method discussed by Beacon (1953). The 
relative weights have been defined with respect to the uncertainty 6 y in the 
first point of the spectrum: 


2 2 2 2 

v =(<5y 1 ) ~ w 2 <<5y 2 ) = w 3 («5y 3 ) 


w m^ y m) * 


The average x is equal to 


— _ Sx 
X Sw 


(72) 


where the summation indices have been suppressed for the sake of brevity. 
The values a and b for the best straight lines are then given by 


a = Swy Swx - Swx Swxy 
Sw Sw (x - x) 2 


(73) 


and b , ?w. assr - 

Sw Sw (x - x) 


The variances of a and b are 


2 2 

, 2 Swx cr 

<6 a ) = ~2 

S wSw (x - x) 

_2 


and ( 6 h ) 


2 = 


^ / -, 2 
Sw (x - x) 


(74) 
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CHAPTER VI 


RESULTS OF THE OBSERVATIONS 

During the present seri.es of eight balloon flights, four of the five cosmic 
X-ray sources which have been unambiguously measured by other workers at 
energies above approximately 15 keV have been observed. Ihe fifth source, 

Cen XR-2, is located in a region of the celestial sphere which was not surveyed 
in this study. 

Two sources, the Crab Nebula (Tau XR-1) and Cyg XR-1, have been 
observed in this series as well as by several other experimenters. Ihe Crab 
Nebula has become a calibration standard for X-ray experiments, while the 
absolute intensity of Cyg XR-1 was reported to be varying in time. Sco X-l 
was also observed during one of our experiments, and in a small number of 
balloon observations by other experimenters; again, the absolute intensity as 
well as the spectral shape of this source appears to change with time. 

We have detected hard X-radiation from a source or source distribution in 
the vicinity of the galactic center. Two other balloon observations have con- 
firmed these findings, but this source region is not yet clearly resolved at high 
X-ray energies. Finally, we have conducted three. observations of the region 
near the north galactic pole. Ihe detector response during the first of these 
observations can be interpreted as being due to a source near the direction of the 
north galactic pole, but subsequent measurements have failed to detect a signifi- 
cant flux of hard X-radiation from this region of the celestial sphere. 
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A. Hie Crab Nebula 

1. Experiment 4 . Two vertical detectors, designated W and A, were 
used in Experiment 4 on 27 April 1966. A modulation disk and ring' were 
mounted in front of telescope W; the detector A was not modulated. 

Figure 12 shews the counting-rate profile for the scintillation crystal 

(22 - 60 keV) of detector W. The modulated exposure curves to Cyg XR-1 and 

to the Crab Nebula are shown in the lower part of the figure. It is evident that 

the Csl crystal counting rate increased prior to the exposure to the Crab 

Nebula. Between 1420 UT and 2020 UT, no known X-ray source was within 

the field of view, and the available information indicates that the floating altitude 

2 2 

remained constant at 2. 8 g/cm , to within 0.2 g/cm . Between 1630 UT and 
2230 UT, the balloon drifted in geomagnetic latitude from a local vertical cutoff 
rigidity of about 4. 7 GV (Shea et al. , 1968) to a rigidity of 3. 3 GV. We inter- 
pret the change in the background counting rate in the Csl crystal as a geo- 
magnetic latitude effect. Data points for 4. 7 and 4. 0 GV, together with data 
points from other experiments at 4. 3, 4.8, and 5.3 GV (Figure 10), were used 
to find a curve of Csl crystal background rate versus the cutoff rigidity 
(section III J). In the absence of a proven model, a simple straight line 
through the low-rigidity data points was used to extrapolate the background 
flux to lower values of the cutoff rigidity. The background flux during the time 
of exposure to the Crab Nebula was inferred from this extrapolation. The 
source contribution was obtained as the difference between the observed 
counting rate and the computed background rate. The resulting net spectrum 
for exposure of the detector A to the Crab Nebula is shown in Figure 13 

~ mu m . , ggg I UlllHI Jil.mi.1 
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together with the best-fitting lines to the data of Peterson et al. (1968), and the 
summary given by Gould (1967). Although the statistical accuracy of the scin- 
tillation crystal points of Figure 13 is relatively poor, the results demonstrate 
that the background-determination method is, in principle, correct. 

The geomagnetic latitude effect for the proportional counter is small and 
somewhat uncertain over the rigidity range considered here. Without any 
correction for this effect, a net spectrum was obtained which is also shown in 
Figure 13. 

2. Experiment 11 . Experiment 11 had a three-fold purpose: First, a 
measurement of the Crab Nebula with higher total closure than in Experiment 
4 was desired; second, a calibration of the detectors and the modulation mecha- 
nisms was needed in view of the new results concerning hard X-radiation from 
the vicinity of the galactic center (section VI D, Exp. 10); and finally, upper 
limits to the X-ray flux from the sun were sought in order to strengthen argu- 
ments against the sun as the source of the observed X-radiation in Experiment 
10 . 

Two detectors were flown in Experiment 11 on 15 July 1967. The detector 
E was vertical at all times and had a modulation ring and disk mounted in front 
of the telescope. The balloon gondola rotated continuously at the rate of one 
revolution per 10 minutes . The zenith angle of detector G was programmed so 
that the telescope axis traversed the Crab Nebula once in every rotation. In an 
effort to calibrate our Scorpius and Sagittarius measurements on a well-known 
source, the detectors G and E were flown without any changes in the pulse- 
height analyzer or the data recording subsystem. 


Hie detector E performed a meridian transit of the Crab Nebula with 
additional modulation during this experiment. A clear signal was observed by 
the Csl scintillation crystal (Figure 14), while the observation by the proportional 
counter had a lower statistical significance. A plot of the observed counting 
rate in the crystal over the energy interval from 22 to 60 keV followed the 
predicted exposure curve (center curve in Figure 14) quite well. 

Hie occultation disk and ring of detector E were designed to provide a 
structured response pattern in the observation of discrete sources. Hie accu- 
racy with which a source position may be determined by this procedure was 
tested with the data from this experiment. 

During the meridian- transit of a source, the occultation disk and ring 
produce a certain counting- rate curve, the shape of which depends on the 
difference in latitude (declination) of the telescope axis and the source. Hie 
position of the pattern is determined by the right ascension of the source and 
the instantaneous longitude of the detector. Hie observed counting rate 
dN/dt(t) of the proportional counter, can be folded with the response function 
E(t,a;,( 5 ) of the detector. Hie function F at any time t is equal to the effective 
exposure area of the counter to a source at right ascension a and declination <5. 

We must now construct a correlation function R(a, t), which has a maximum, for 
the counting- rate profile dN/dt(t) caused by a source at (a,<5), if this profile is 
identical to that expected from a source at (a, 6 ). A suitable normalized folded 
counting rate is 
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4f(t) F (t, o;, 6) dt 


’£()+ At 


R (a, a) 


#(*) dt 


ST 

f F (t, a 


,6) dt 


Using data from the present observation, we scan in declination for 
sources at a fixed right ascension a = 5^3i m . The folded rate, shown in the 
upper part of Figure 15, is necessarily symmetrical about the declination of the 
detector axis (approximately 31°). Two maxima of the curve are seen at decli- 
nations <5 » 21. 5° and <5 «40. 5°, but because of the known positions of cosmic 
X-ray sources in this region of the sky, the second peak may be ignored. As a 
test of the method, an artificial counting- rate pattern for a point source at the 
position of Tau XR-1 was generated, and the above-mentioned procedure was then 
applied to this idealized pattern. The resultant curve is shown in the lower part 
of Figure 15 and indicates close agreement with the correlation function of the 
observed pattern. 

Errors in the determination of the detector exposure area as a function of 
the zenith angle and the geographical position of the balloon gondola limit the ac- 
curacy of this method to approximately ±2° in declination and right ascension. 
Errors caused by a possible misalignment of the detector axis and the modulation 
gear are not important in this experiment because of the continuous rotation of 
the gondola. 

The net counting rate due to the Crab Nebula was analyzed with the help of 
procedures 1 and 2 (chapter V). The results of an extrapolation according to the 
second method are shown in Figure 16. Diamonds represent the proportional 


counter data, and double-bar crosses show the scintillation crystal data of 


61 

detector E. Single-bar crosses represent the scintillation crystal data of 
detector G. 

The scanning method of the detector G was described in chapter IV. 

Because of some minor difficulties in the data handling for certain rotations, 
only 21 rotations with unambiguous azimuth information were used for the 
analysis. Thn zenith angle during each of these 21 rotations was such that the 
detector aiternatingly viewed the Crab Nebula and regions of the sky which are 
void of any known emitters of X-radiation. During certain periods the sun was 
within the field of view. However, sun exposure and exposure to the Crab 
Nebula were clearly separated in azimuth angle and were therefore easily 
distinguished in viewing time. 

The data from each rotation were grouped into 10 1-minute bins, each 
corresponding to an azimuth angle interval of 36°, and superposed so that bins 
from successive rotations with maximum exposure to the Crab Nebula were 
added. Then the bins with 18° -54° azimuth angle with respect to the Crab Nebula 
were added, then the 54° -90° bins, etc. The resulting curves of observed count- 
ing rate versus azimuth angle with respect to the X-ray source are given in Fig- 
ure 17 for the scintillation crystal and in Figure 18 for the proportional counter. 
In each figure the top two curves (A and B) refer to Experiment 10, and the next 
4 curves (C, D, E, and F) are relevant to the Crab Nebula observations of 
Experiment 11. The exposure to the source was calculated from the instanta- 
neous detector positions; the resulting curves (B, D, and F) are drawn below 
the counting- rate superpositions (A, C, and E). Crab Nebula data from a group 
of seven rotations (C) are shown for this experiment in order to allow a 


comparison with the seven-rotation exposure to the vicinity of the galactic 
center (Experiment 10). The full Crab Nebula data set for 21 rotations is 
represented by curves E in Figures 17 and 18. 

The proportional counter ?’esults indicate that, although the background 
level was higher, there was a clear signal in Experiment 10 (curve A in 
Figure 18) and no significant signal in the present experiment (curves C and E 
in Figure 18). The absence of a signal is statistically significant, i.e. , the 
observed number of counts in seven rotations within 54° azimuth angle from the 
direction to the Crab Nebula is approximately 4cr below the value which is 
expected from an extrapolation of the Csl scintillation crystal spectrum. 

Data from the scintillation crystal of detector G were again analyzed 
according to both procedures 1 and 2. The extrapolated spectrum at the top 
of the atmosphere is given in Figure 16, together with the results from 
detector E. 

The observed counting-rate distribution is due to one specific source 
spectrum. However, different methods of data analysis result in different 
numerical values for the spectral parameters, for the statistical uncertainty, 
and for the relative probability of a particular spectral shape (chapter n). An 
intrinsic property of the analysis procedure 2 (section V B) is that the statistical 
•significance is lower than in procedure 1 because of the smaller number of 
degrees of freedom. Analytic methods of procedure 2 can be used for power- 
law and exponential fits to the data, in contrast to those of procedure 1, where 
a figure of merit is used to determine the best fit. On the other hand, the 
tatistical error for transformed straight-line fits is calculable with standard 


procedures, and in the absence of a known method for the determination of the 
statistical error for chi-square tests, we have used the method described in 
section V A. 

We have analyzed data from three independent detectors: the proportional 
counter in E and the scintillation crystals in E and G. The f-values for the 
same parameter combination were added according to the reproductive law of 
the chi-square distribution. The parameters of the best fits, their statistical 
uncertainties and relative probabilities (where applicable)* are summarized in 
Table VI. Whenever procedure 1 resulted in equally high values of Q for several 
combinations of the parameters A and n, then a value of A close to the "best" 
value for the amplitude A from procedure 2 was adopted. We have tabulated 
the results for one data set, from the Csl crystal of detector G, and the results 
for the sum of all three data sets. The table illustrates that the resulting 
parameters and their statistical uncertainties are relatively stable when two 
distinctly different analysis procedures are applied. 

None of the three model spectra can be recognized as the correct one on 
the basis of this analysis alone. However, a thermal spectrum with the Gaunt 
factor for hydrogen-like atoms and with a correction for electron-electron 
collisions (chapter II), seems to be favored by a few percentage points. The 
same observation was made by Jacobson (1968) in his analysis of a high-energy 
measurement of the Crab Nebula. We have come to a similar conclusion in the 
analysis of our observations of the vicinity of the galactic center (section VI D) . 

The present measurements and observations by other experimenters at 
energies between 1 and approximately 500 keV indicate that the spectrum can be 
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explained by a single power-law with an index n 2. 0, In particular, a recent 
rocket observation of the spectrum of the Crab Nebula between 2 and 20 keV by 
Boldt et al. (1968b) resulted in a spectral index n = 1. 93 ± 0. 05, The fact that 
the Crab Nebula spectrum between 2 and 60 keV is consistent with a power-law 
of index n = 1. 93 appears to rule out a two- component thermal model as proposed 
by Sartori and Morrison (1967a). In a recent study (Holt et al. , 1968), it has 
been shown that the expected line emission from such a thermal source should 
be evident. Data from the above-mentioned rocket experiment were used to 
exhibit that, at a statistical level exceeding 1 <t, there is no evidence of such 
line emission for any relevant model of elemental abundances. The sum of all 
available information on X-radiation from the Crab Nebula indicates that the 
observed power-law spectrum is caused by synchrotron radiation. 

B. CygXR-i 

1. Experiment 3 . Experiment 3 was launched on 13 Jan. 1966. One 
X-ray telescope with modulation gear was suspended vertically so that it exe- 
cuted a meridian- transit scan of Cyg XR-1. The observed counting-rate 
distribution has been folded with the pattern expected from a source at the 
position of Cyg XR-1 (compare section VI A). The results showed that the most 
probable source location of the observed X-radiation is indeed consistent with 
the reported position of Cyg XR-1 (Bowyer et al. , 1965; Giacconi et al. , 1967a), 

The net spectrum of Cyg XR-1 was determined from the proportional and 
scintillation counter data. The net flux at 25 keV, where the spectra of the two 
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detectors overlap, was extrapolated to the top of the atmosphere. The resulting 
primary photon flux is J(E = 25 keV) « 0. 01 photons/cm -sec-keV. 

2. Experiment 4 . Altitude information was not available during the 
exposure to Cyg XR-1 in Experiment 4, but the pressure profile prior to 
12:30 UT, indicates that floating altitude was reached near the beginning of 
this period. 

The response of the scintillation detector of the modulated telescope W 
to Cyg XR-1 is shown in Figure 12. The observed counting-rate pattern (upper 
curve) follows the predicted response (lower curve) quite well. The decrease 
in the background rate between the end of the exposure to Cyg XR-1 and approx- 
imately 1630 UT cannot be accounted for. A reduction of the event rate would 
be expected to be caused by the change in the geomagnetic latitude which 
occurred during that period, but the observed decrease is too strong to be ex- 
plained in this way. Altitude changes may have caused the observed pattern; 
however, the lack of altitude information precludes confirmation of any explana- 
tion of this type. 

Data analysis procedure 2 was employed for the evaluation of the Cyg XR-1 
data. Because of the lower exposure factor of detector W for the Cygnus ob- 
servation, the results from this detector have a lower statistical significance. 
The extrapolated spectrum from detector A is shown in Figure 19. It is 
apparent that the proportional counter information below 20 keV is not con- 
sistent with the remaining points of the extrapolated spectrum, and it has 
therefore been omitted in the following analysis. 
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A weighted least-squares fit of a power-law spectrum, 

J(E) = A-E‘ n , (76 ) 

2 

(J in photons/cm -sec-keV, E in keV) yields A ~ 7.47 and n = 2. 0 ± 0. 3. The 
best fit of an exponential spectrum, 

J(E) = A'*exp (-E/kT) / E , (77) 

was obtained for A' = 0. 53 and kT = 34. 2 +12 ‘ 4 keV. 

- 7.2 

These parameters may be compared with the results of other observers: 
for example A = 3. 58, n = 1. 93 ± 0. 20 (Peterson et al. , 1968, experiment on 
13 Sept. 1966); A = 16, n = 2. 2 ± 0. 8, kT = 35 (Bleeker et al. , 1967, experi- 
ment on 5 April 1966); n = 1. 6 ± 0. 2, kT = 60 ± 20 (Clark et al. , 1968, experi- 
ment on 19 July 1966); and A = 3. 4, n = 1. 8, kT = 94 (Haymes et al. , 1968, 
experiment on 29 Aug. 1967). Haymes et al. (1968) have observed a steepen- 
ing of the spectrum of Cyg XR-1 above 140 keV; below this energy a power- law 
spectrum is consistent with the available observations. Thus, the X-ray 
source Cyg XR-1 may be of the same type as Tau XR-1. 

The intensity of Cyg XR-1 is reported to vary in time (Overbeck and 
Tananbaum, 1968). While the statistical accuracy of our first observations of 
Cyg XR-1 is insufficient to make a clear comparison, we can say that the 
intensity at 30 keV observed in the present experiment (27, April 1966) is 
consistent with the spectrum measured by Bleeker et al. (1967) on 5 April 1966. 

C. Sco X-l 

The purpose of Experiment 9 was to measure X- 
The experiment was col ^ted over Mildura, 
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Australia, on 15 Dec. 1966. The detector E, with modulation disk and ring, 
was mounted on the balloon gondola so that the telescope axis formed an angle 



0 of 20° with the vertical. The gas-jet orientation system was used to align 
and maintain the horizontal projection of the detector axis at a fixed angle (3 
with respect to the local direction of the horizontal projection of the geo- 
magnetic field lines. According to the information given on aeronautical maps 
of the region, the angle (3 was set at +8.7°, i. e. , the local magnetic field line 
is shifted by (3 with respect to the meridian line, and magnetic north appears 
to be 8. 7° east of geographical north. Subsequent studies of satellite data 
indicated that f3 is actually equal to +7. 6°, thus producing an error, A/3 ? = 1. 1° , 
in the magnetometer setting. However, the effective error A (3 at zenith angle 

©is considerably smaller than A/3 T : 

2 2 

cos A /? = cos 0 + sin 0 cos A/3'. (78) 

A small compass, which was photographed together with the barograph at 
intervals of 2. 5 minutes during the flight, indicated that the gonauia was 
oriented properly throughout the experiment. 

This particular balloon flight had an unusually slow ascent rate with the 
result that floating altitude was reached after the meridian transit of Sco X-l, 

and the changing altitude during the exposure to the source consequently made 

* 

the background subtraction more uncertain than in Other experiments. We have 
used the altitude profiles of other experiments to infer the background counting 
rate during the last part of the balloon ascent. 

/" Figure 20 shows the difference between the observed counting rate and 

| ‘~N ' 

V the assumed background. The histograms represent the net rate for 
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10-minute intervals without any adjustment for the varying altitude. The relation 
between the residual atmospheric depth and Universal Time is shown on the 
scale at the bottom of Figure 20. Data from the Csl crystal (22 - 60 keV) are 
shown by the top curve, and data from the proportional counter are represented 
by the second curve. The statistical uncertainty for each 10-minute interval 
is typically 55 counts/10 min for the Csl crystal data and 30 counts/10 min for 
the proportional counter data. The continuous line (third curve) represents the 
exposure area (cm 2 ) for Sco X-l. It was derived from the instantaneous 
detector- source configuration under the assumption that the detector axis was 
in the plane of the meridian under a zenith angle of 20°. This curve is valid 
for the proportional counter at 25 keV, and the appropriate curves for the 
scintillation crystal at other energies are similar. 

There is apparently a shift of about 15 minutes between the predicted 
profile and the observed histograms. This difference must be interpreted as an 
error of approximately 4° in the detector alignment with respect to the magne- 
tometer axis. With this interpretation we can now shift the curve of exposure 
area by the proper amount and compute the net source spectrum from Sco X-l. 
Data analysis procedure 2 was used to correct the net counting rate for the 


detector efficiency, and to extrapolate it to the top of the atmosphere. The 

correction factors (chapter V) were determined under the assumption of 

7 

thermal emission from an optically thin gas at T = 5 • 10 ° K. the resulting 
energy spectrum is plotted in Figure 21. The proportional counter results 
are represented by diamonds, and the Csl crystal data are drawn as crosses. 
It must again be pointed out that the absolute values of the spectrum points 
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in Figure 21 are somewhat uncertain because of the uncertainty in the in- 
stantaneous background flux during the exposure to Sco X-l. The shape 
of the resulting spectrum, however, as well as the relative positions of pro- 
portional and scintillation counter results in the spectrum are more positive. 

The broken line in Figure 21 represents the energy flux due to thin-source 
bremsstrahlung of a gas at T = 5* 10 °K. Thermal spectra of this type can 
explain several of the lower-energy observations of Sco X-l (Chodil et al. , 

1965; Grader et al. , 1966; Gorenstein et al. , 1967; Chodil et al. , 1968; and 
Gorenstein et al. , 1968). Figure 21 shows that our measurement is also con- 
sistent with a thermal spectrum at energies below approximately 35 keV. At 
higher energies, the observed flux is definitely above the extrapolation of the 
low-energy spectrum. 

A similar high-energy feature was measured by Peterson and Jacobson 
(1966b) in a balloon experiment on 18 June 1965. The reported spectrum is 
drawn schematically in Figure 21. From 20 to approximately 35 keV, the data 
fit an exponential law corresponding to T = 5* 10 °K. Above 35 keV, in 
particular around 45 keV, the energy flux appears to be significantly higher 
than the extrapolation of the thermal spectrum would indicate. Our results are 
in substantial agreement with the observations of Peterson and Jacobson (1966b). 
Whether this feature is a flat high-energy tail or a single line at approximately 
45 keV cannot be decided upon from these two observations. Both spectra are, 
however, not consistent with power-law spectra, so that we can rule out 
synchrotron emission from Sco X-l. This statement is supported by the fact 


that the optical polarization of Sco X-l (Hiltner et al. , 1967) can be explained 

as due to the interstellar material alone. 

Experimenters at the University of Adelaide also have observed an 

apparent high-energy tail in the Scorpius spectrum during an experiment which 

was conducted on 29 Feb. 1968 (Buselli et al. , 1968). The best fit to 

their data is provided by the power-law spectrum: 

-2 79 

1(E) = 925- E * , (79) 

2 

where I is the energy flux in keV/cm -sec-keV and E is measured in keV. This 
function is drawn in Figure 21 as a dashed-dotted line. The data may also be 

g 

fitted to a thermal spectrum for T «1. 7- 10 °K. In either case, the observed 

high-energy flux is definitely above that corresponding to the 5* 10 °K thermal 

spectrum which can be inferred from measurements at lower energies (between 

,2 

1 and 30 keV). In fact, the average energy fluxes (in keV/cm -sec-keV) in 
the energy interval between 40 and 50 keV were 0. 06 as observed by Peterson 
and Jacobson (1966b); 0. 07 in the present experiment; and 0. 03 as observed 

by Buselli et al. (1968). In contrast, the measurement by Lewin et al. (1967) 

.* , 

on 13 Feb. 1967 yielded a lor-upper limit of 0. 011, while the observation by 

o 

Overbeck and Tananbaum (1968) yielded an upper limit of 0. 008 keV/cm -sec- 
keV. The upper limit to the energy flux prior to the X-ray flare on 15 Oct. 

1967 was given by Lewin et al. (1968c) as approximately 0. 01 keV/cm -sec- 
keV. 

2. The stellar accretion model for Sco X-l . Several authors (Shklovsky, 
1967; Cameron and Mock, 1967; Prendergast and Burbidge, 1968; and Westphal 
et al. , 1968) have built models for the optical and X-ray emission from Sco X-l. 


; 
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; 

From energy-balance arguments they sho;y that energy must be supplied con- I 

tinuously, for example by infalling matter through a conversion of gravitational f 

potential energy into kinetic energy. | 

Shklovsky (1967) has proposed a model for Sco X-l in which gas streams 
from the secondary component of a close binary system toward a compact 
(primary) star. The latter might be a neutron star, a white dwarf, or any 

$ 

other massive star with suitable parameters (Cameron and Mock ? 1967; [ 

Prendergast and Burbidge, 1968). During the mass accretion of the primary 
star, gravitational potential energy of the incoming matter is released and 
converted into kinetic energy. The coronal plasma is thus heated to a high 
temperature, and thin-source bremsstrahlung can occur. Some of the 
resultant X-ray emission can leave without significant absorption (the overlaying 
material is optically thin) while another part may be absorbed in the photo- 
sphere and reradiated in a black-body spectrum corresponding to a lower 
temperature. 

In Shklovsky J s model, a high-energy tail in the Sco X-l spectrum, as j 

observed by Peterson and Jacobson (1966b) and in the present experiment, may 
be explained readily as the thermal emission of a still hotter plasma, e. g. , at 
a temperature of the order of 5* 10 °K. 

The maximum coronal temperature that can be reached through the 
release of gravitational energy by the infalling gas is given by 

3/2 Nk T = G M \i/ R , (80) 

where the mean molecular weight is 2/3. With the stellar mass M and 
the radius R measured in solar units, we obtain 
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7 M /BA 

T = 1, 02* 10 (81) 

As Cameron and Mock (1967) and Prendergast and Burbidge (1968) have 
pointed out, temperatures of T » 10 - 10 °K will be obtained from binary 
systems where the normalized value of M/R is between 1 and 10, such as in 
white dwarfs. For a neutron star, however, M/R (R 0 /M 0 ) « io^, an d t can 
be as high as 10 °K. The observation of temperatures of the order of 5* 10 °K 

would then extend the range of possible values for the normalized ratio M/R of 
the accreting star. 

As we have stated above, no high-energy tail was observed on 13 Feb. 

1967 by Lewin et al. (1967). If all the observations are correct, then we must 
assume that the influx of mass to the primary star decreased between the time 
of our measurement, 15 Dec. 1966, and the next observation, 13 Feb. 1967. 

This change over a period of two months must have been sufficient to decrease 
the flux of high-energy X-radiation by a factor of 6 or more. Between 15 Oct. 
1967 and 29 Feb. 1968, the high-energy X-ray flux apparently increased again 
by a factor of 3 or more. In this connection, it may be interesting to note that 
the energy flux at 25 keV on 13 Feb. 1967 (Lewin et al. , 1967) was approxi- 
mately one-half of the flux that was measured on 18 June 1965 by Peterson 
and Jacobson (1966b). 

Such fluctuations in the observed spectrum of Sco X-l are not surprising 
in view of the fact that an intensity change by a factor of 2 over a period of one 
month was measured by Overbeck and Tananbaum (1968). Changes in the 
absolute flux level of low-energy X-radiation as well as in the spectral shape 
have been reported by other observers (Chodil et al. , 1968; Gorenstein et al. , 



1968; and Hill et al. , 1968). Several observations of fluctuations in the optical 
spectrum of Sco X-l have also been reported (Chodil et al. , 1968), 

3. Bremsstrahlung from an isovelocity plasma. Bless et al. (1968) and 
Wallerstein (1968) have suggested that Wolf-Raj'et and Of stars may be sources 
of X-radiation. 

Wolf-Rayet stars are an interesting type of objects whose line-emission 
spectra show some features which are also observed in O stars. However, 
they deviate markedly from typical O spectra in that the He lines are stronger. 
Some objects, called ”WC stars,” show C, O and strong He lines, while the 
”WN stars” show strong N and He lines and weak C and O lines. 

UV observations of Wolf-Rayet and Of stars have shown that these objects 
lose mass with velocities up to 2000 km/sec. Although it is not clear whether 
these velocities are associated with gas streams In a binary system or with the 
symmetrical expansion of a hot corona; coronal temperatures of the order of a 
few • 10 °K have been inferred. The high temperatures may also be generated 
through interaction between the outflowing gas and the interstellar medium. 

It has been attempted to compare the distribution of X-ray sources with 
that of associations (O’Dell, 196.7), but such a correlation does not account for 
the absence of X-ray sources in the region 1^ = 120° - 260° (we exclude the 
Crab Nebula in this argument). Wolf-Rayet stars show a better correlation 
since they are very scarce in the above-mentioned interval of galactic longi- 
tude; however, a detailed statistical analysis has not been carried out. 

Two X-ray sources have been tentatively correlated with Wolf-Rayet 
stars. The optical spectrum of GX 3+1 (see section VI D 3) is similar to that 
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2 8 
of a WC star, while the Wolf-Rayet star y Vel was recently observed in soft |; 

It 

X-radiation (Grader et al. 1968), jj 

Jl 

Fischel and Stecher (1908) have suggested that the observed X-radiation || 

1 ! 

Jj 

may be due to bremsstrahlung from an isovelocity plasma which interacts with 

* ! 

a stationary medium. The heavier ions have greater kinetic energy and there- 
fore produce more energetic photons during their thermaiization. In such 

j. 

a model, the abundance distribution as well as the degree of ionization play a 

crucial role. For any reasonable abundance model, helium is the main con- l 

f; 

f; 

tributor to the spectrum at typically 3 keV while heavier ions at velocities of a 

few hundred km/sec produce the bulk of the observed radiation at energies around f 

j 

50 keV. j 

The mathematical formalism of Fischel and Stecher (1968) was used to \ 

find the total X-ray flux emitted during the thermaiization of an isovelocity 
plasma with a universal abundance distribution (Aller, 1961) for the 10 most 
abundant elements up to iron. Fischel and Stecher (1968) found that this 
plasma, with degrees of ionization as derived from optical and UV observa- 
tions of hot stars, provides the best agreement with the spectral shape of 
Sco X-l as determined by Grader et al. (1968), and by Buselli et al. (1968), 

for a velocity of 380 km/sec. On the other hand, we find that an isovelocity j 

£ 

plasma at v = 200 km/sec provides good agreement with the data of Goren- 

• : ' • • ' ■ : t 

stein et al. (1968), Chodil et al. (1968), and the results of the present ex- 
periment (Figure 21a). Sample calculations for other model plasmas, for j 

example with enhanced helium and heavier-element abundances, result in j 

t 

; . .. [ 

spectral shapes which are not consistent with the experimental spectrum 

• | 


shown in Figure 21a. 


We have attempted to understand the observed fluctuation in the optical 
and X-ray spectrum of Sco X-l within the framework of this model. An 
increase in the available energy may increase the degree of ionization and 
therefore increase the intensity of high-energy bremsstrahlung by shifting 
the spectral contribution of each element to higher energies. Calculations for 
the model plasma which. was discussed above (curve 3 in Figure 21a) show that 
upon an increase of 2 in the degree of ionization for the more abundant ele- 
ments from magnesium to iron, the photon flux at 3 keV increases by 14% 
while the flux at 50 keV increases by a factor of 2. 25 (curve 4 in Figure 21a). 

Such a difference in the photon flux at high X-ray energies is apparent from a 
comparison of the results of the present experiment with the data of 

Lewin et al. (1967) . We are not aware of any low-energy X-ray measure- 
ments which were taken during the times when the two sets of high-energy 
X-ray observations were carried out. We have, however, plotted in 
Figure 21a the results of observations made on 11 Oct. 1966 (Gorenstein 
et al. , 1968, curve 2), and on 29 Sept. 1967 (Chodil et al. , 1968, curve 1). 

The photon flux at optical wavelengths is not affected by the partial 
change of the state of ionization in the particular model plasma used here, 
but may increase somewhat if the relative abundances of helium and the 
heavier elements is increased. 

It is apparent that the presently available data are not sufficiently 
accurate to confirm or rule out the validity of the present model. Complete 
spectra in the optical, low- and high-energy X-ray region are needed in 
order to build consistent models for the chemical composition, the state of 
ionization and the velocity of the streaming plasma. 
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D. The VicinJty of the Galactic Center 

1. Experiment 10. The purpose of Experiment 10 was to scan the 
vicinity of the galactic center in hard X-radiation, This region of the sky 
was chosen for balloon observations because several sources of low-energy 
X-radiation had previously been detected there in rocket experiments. 

The experiment was conducted over Australia on 18 Dec. 1966. The 
orientation system for the detector G, and the modulation disk and ring for 
detector E were described in the discussion of Experiment 11 (section VI A), 
The zenith angle of detector G was programmed so that the detector axis passed 
approximately tangentially through the galactic plane once upon each rotation of 
the balloon gondola. 

Six passes from galactic north as well as a group of six passes from 

galactic south are drawn schematically in Figure 22. A few relevant sources 

of low-energy X-radiation are also given there. In the manner described in 

section VI A 2, the raw counting rate was grouped and superposed into 10 bins 

of equal azimuth angle with respect to a fixed point in the celestial sphere. 

This point was assumed to be the X-ray source Sgr XR-1, which was detected 

by Friedman et al, *(1967) (a= 17 h 55 m , 6 = -29. 2°). The results Of the 
superpositions for two groups of 6 rotations are shown in Figure 23. Both 

the passes from the galactic north and south exhibit a minimum in the counting 

rate for azimuthal angles near 180°. The dominant feature of the superposition 

of passes from the galactic north was the signal detected from Sco X-l by the 

krypton proportional counter. In contrast, the dominant feature of the super- 
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position of passes from the galactic south was the signal from the vicinity of 

Sgr XR-1, detected by the scintillation crystal. This is a strong qualitative 

indication of the relative hardness of the two spectra. In both cases, we 
derived the positional information from the azimuth angles with respect to the 

direction to Sgr XR-1. 

A total of seven passes from the galactic south was used for the evaluation 

of the X-ray flux from the vicinity of the source Sgr XR-1. The appropriate 

data superpositions and a comparison with measurements of the Crab Nebula 

are presented in Figures 17 and 18 (section VI A 2). Two different procedures 

(chapter V) were used for the data analysis. The results of procedure 2 are 

given in Figure 24, in which diamonds represent the proportional counter data, 

and crosses symbolize the scintillator crystal data. A power-law of index 

n = 2. 69 ± 0. 92 (dotted line) provides the best of such a spectrum. The best 

fit for a bremsstrahlung spectrum from an optically thin, hot plasma is 

+4 7 

obtained for a temperature corresponding to kT = 16. 0 * keV (solid line in 

"u« U 

Figure 24). The numerical results of both analysis procedures are summarized 
in Table VII. The thermal spectrum for Z = 1, including electron-electron 
collisions (chapter n) , appears to provide the best overall fit to the data from 
the vicinity of the galactic center. 
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2. The sun as source of the observed X-radiation (Experiment 10) . It is 
apparent from Figure 22 that the sun ( a~ 17 4i m , 6 = -23° 22’) was within 
the field of view during exposure to the vicinity of the galactic center. How- 
ever, it is believed that the sun’s contribution to the observed counting rate 
was negligible. The following arguments support this belief: 

(a) The response of detector E was used to determine the location of 
the source, or the center of the source distribution, which gave rise to the 

observed radiation. The shape of the modulated response depends strongly 
on the difference in declination between the detector axis and the X-ray source. 

We have therefore computed the folded rates F (a, 6 ) for a fixed right ascension 
of oi = 17^55 m . A scan in declinatior with the data from this experiment is 
shown in the upper part of Figure 25. Counting-rate distributions were gene- 
rated for hypothetical sources at the positions of Sgr XR-1 and the sun. The 
folded rates for the expected profiles are given in the lower part of Figure 25 

i 

for a source at the position of Sgr XR-1 (solid line), and for one at the location 
of the sun (broken line). The left peak of the solid line can be ignored because 
of the absence of reported sources in the vicinity of -39° declination and at 
this right ascension. Figure 25 indicates that the declination of the source or 
the center of the source distribution is closer to -29° than to -23° . A single 
source or k source distribution centered at the position of the sun seems to 
be ruled out. 

(b) At one time during Experiment 11, the detector G had full exposure 
to the sun. The response in the 22 to 60 keV band deviated from background 
by less than Icr during exposure to the sun. 1 This deviation would account for 

v : 
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approximately 10% of the observed response from the vicinity of the galactic 
center over the same energy interval. Indirect evidence supporting the 
validity of this correlation is supplied by the fact that the daily average of the 
9400 MHz solar radio emission (Nagoya University, 1967), which usually 
shows strong correlation with long-time changes in the soft X-ray emission, 
was constant to within 5% on the days on which Experiments 10 and 11 were 
conducted. 

(c) An importance 1-flare. (GRP 2160) was reported (ESSA, 1967) towards 

the end of the one-hour data-collection interval of Experiment 10 which was 

used in the analysis. No change in the counting-rate pattern followed this flare. 

The planetary magnetic Index, K , was 2- during exposure to the vicinity of the 

P 

galactic center. The experiment was conducted on a geomagnetically quiet day. 

3. Discussion of the data from the vicinity of the galactic center . The 
vicinity of the galactic center has been scanned from rocket altitudes by 
several experimenters. The most recent publications of four research groups 
(Bradt et al. , 1968a; Gursky et al. , 1967; Fisher et al. , 1968; and Friedman 
et al. , 1967b), as well as the results of a recent balloon observation by Buselli 
et al. (1968) have been used in the analysis. Figure 26 shows the relevant 
reported positions of X-ray sources which emit at photon energies of typically 
2-6 keV. Except for the thin- slit observations of the ASE group (Gursky et 
al. , 1967), the uncertainty in the source position is not shown in this graph. 

The broken lines represent the trajectories of the axis of the detector G during 
the first and the last of the seven rotations which were used in this analysis. 


80 


In Figure 26, considerable confusion as to the true positions of X-ray 

sources exists. The situation is relatively clear, however, at galactic lati- 
II II 

tudes between 1 =0° and 1 « 20°, where we have adopted the positions of 
5 sources as reported by the MIT group (Bradt et al. , 1968a). For galactic 
longitudes between 340° and 360° , three source positions were, perhaps some- 
what arbitrarily, selected so that they show maximum correlation with the 
reported positions. These sources lie on the "position lines" reported by the 
ASE group and will possibly have to be revised as more reliable observations 
become available. For source number 2 (Table VIII), we have adopted the coordi- 
nates of GX 354-5 (Buselli et al. , 1968). Table VIII lists the eight sources thus 

selected. The source designations of the various experimental groups are given, 

o 2 

along with an estimate for the source flux density, J , in photons/cm -sec- 

s 

keV at a photon energy of E « 4 keV. For the source GX 3+1, we have adopted 
the position of an intense ultraviolet object of magnitude V = 13.4 (Blanco et al. , 
1968a, and Blanco et al. , 1968b). The optical spectrum of this star is consistent 
with a Wolf-Rayet type spectrum (Freeman et al. , 1968). (Compare section 
VI C 3). 

We can now attempt to compare the observations reported here with the 
results from rocket-borne instruments. Because of the relatively large angle 
of view of our X-ray detectors, direct comparison of the results with the 
spectra or flux-values for individual sources is not possible. 

Let us consider a hypothetical detector which is suitable for low-energy 
X-ray measurements and has the same angular response as our detector. We 
can compute the photon flux which this detector would have observed if it had 
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followed the same trajectory as our detector. Obviously, the hypothetical 

low-energy detector would have to be located above the atmosphere. 

We compute the angles 0 . between the source s and the direction of 

s, 1 

the detector axis during the half-minute interval i for M = 8 sources and N = 56 

data accumulation points. With the relative exposure factors F (0 .), we 

s, 1 

determine the flux 


_o 


J . = J . F ( 0 .), 

s, l s ■ s, r 


(82) 


where the estimated source flux density at 4 keV, J , is taken from Table VIII. 

s 


The net flux due to each source, 
1 N 

J = r T 2 J 


s N ,• 


= 1 s, 1 


(83) 


is listed in Table VIII. The total net flux from all sources is then 

M 2 

J= 2 J = 0. 78 photons/cm -sec-keV. 
s = 1 s 


(84) 


This photon flux is shown in Figure 24, together with the results of Ex- 
periment 10. A thin-source bremsstrahlung spectrum with kT = 11 keV pro- 
vides the best fit to the data measured by Lewin et al. (1968b), in a balloon-borne 
experiment. The appropriate curve is also shown in Figure 24 as a broken 
line. Experimenters at the University of Adelaide (Buselli et al. , 1968) 
measured the spectra of GX 3+1 and GX 354-5 on 29 Feb. 1968, and their data 
from GX 3+1 can be fitted either to a power-law spectrum of index 2. 0 (dash- 
dotted line in Figure 24) or to a thermal spectrum with kT = 55 keV. The 
photon flux of GX 3+1 is lower than that observed in this experiment, for 
example, by a factor of 2 at 25 keV. The spectrum of GX 354-5 is also lower 
than that measured here, except for a narrow interval of energy around 20 keV 


where the flux from GX 354-5 would be sufficient to explain the flux observed 
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here. Buselli et al. (1968) have not speculated on the spectrum of this source; 

|; 

if the 20 keV-data point is real, then the X-ray spectrum might be almost flat 

between 2 and 20 keV. f 

i 

Comparisons between the observed spectra or their line representations 

* 

in Figure 24 suggest that a solid-angle effect exists in the various measure- 
ments. Our detector telescope had the largest opening angle and observed the 
highest photon flux; the detector of the University of Adelaide group was most 
strongly collimated and thus observed single- source spectra of lower intensity; 

j 

the observations of Lewin et al. (1968b) are intermediate in both respects. j 

j 

This seems to support the suggestion that we observed hard X-radiation from j 

I' 

• . ¥. 

several sources in the vicinity of the galactic center having, perhaps, different [ 

degrees of spectral hardness. Of course, we cannot rule out the hypothesis that 

■ ■ ' I 

only one source with variable hard X-ray emission was observed. Freeman 
et al. (1968) have given evidence for the short- and long-term stability of 
the optical counterpart of GX 3 + 1, which was the most intense X-ray emit- 
ter during the survey of the vicinity of the galactic center by Buselli et al. 

(1968). 

•r • 

If it is assumed that the observed hard X-ray flux is due to the set of | 

| 

sources whose soft X-ray contribution has just been evaluated* then Figure 24 

f 

demonstrates that it is not possible to decide between synchrotron radiation j 

■■ - -•■■■ ... jr 

• ... '■ ■ .. • I 

and bremsstrahlung as the dominant emission process for the total source 

. ’ . • ■ . . t. 

! 

spectrum. j 

I 

: j 

The ASE and Lockheed groups have reported apparent temperatures for I 

some sources in the vicinity of the galactic center which emit bremsstrahlung 
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from an optically thin, hot plasma. The apparent temperatures range from 
kT - 2. 1 keV to kT = 4. 3 keV (Gursky et al. , 1967) ancl from kT = 6.0 keV to 
kT = 13, 8 keV (Fisher et al. , 1968), for the sources 1, 5, 6, 7 and 8 (Table Vin). 
A wide spread of apparent temperatures can be expected when the emitting gas 
has a temperature distribution rather than one definite temperature. 


E. The Vicinity of the North Galactic Pole 

1. Experiment 1 . Experiment 1 was conducted on 6 Dec. 1965 with one 
vertical X-ray telescope, equipped with a modulation ring and disk. The 

a 

detector axis scanned a band of the celestial sphere at declinations between 
31° and 28° and at right ascensions between 10^12 m and 18^26 m , i. e. clearly 
above the galactic plane. The pulse-height analysis system for the scintillation 
detector malfunctioned during this flight so that the analysis had to be restricted 
to data from the proportional counter. 

An increase in the counting rate was observed between 1400 and approxi- 
mately 1500 UT, corresponding to a celestial position in the vicinity of the north 

galactic pole. Possible altitude changes as the cause for this feature can be 

2 

ruled out, since the floating altitude was stable at 2. 9 g/cm to within less than 
2 

0. 1 g/cm , during the main part of the experiment. 

According to radar information, the balloon drifted in a northeastern 
direction between 1340 and 1410 UT, corresponding to a change in the local 
geomagnetic rigidity cutoff from 4. 95 GV to 4. 5 GV (Shea et al. , 1968). Before 
1340 UT and after approximately 1600 UT, the local geomagnetic rigidity cutoff 
was approximately 4. 9 GV. The shape of the observed counting rate pattern 


\ 

j! 

I 


i 

£ 

t 


? 


I 

l 

l 


! 



r, 

i; 


rt 11 


84 

differs from the background rate pattern as deduced from the change of the 
local cutoff rigidity. A definite cutoff- dependence of the atmospheric and 
detector background was established for the scintillation crystal from several 
experiments. The corresponding data for the proportional counter contradict 
each other somewhat, and the dependence on the rigidity cutoff is unclear 
(section III J). 

The identification of the actual source position makes use of the particular 
geometrical features of the X-ray telescope and the modulation disk and ring. 

This method was tested in Experiment 11 and found to be quite accurate in 
determining the right ascension and declination of an X-ray source. 

The definition of the folded rate R (a, 6 ) was discussed in section VI A. 
Figure 27 shows the functions R (o>, 6 ) computed with the counting- rate patterns 
observed in Experiment 1. The upper part shows the curve R (a - 13^00 m , <5 ) 
as a function of declination. The lower section of Figure 27 represents the 
curves R {a, 6 = +28°) (solid line) and R ( a , 6 = +31° ) (broken line) as a function 
of right ascension. 

An intrinsic feature of this type of analysis is that the function R ( a, 6 ) 
for a fixed a is symmetric in declination with respect to the declination of the 
detector axis. The broad structured peak in the upper part of Figure 27 indi- 
cates that the source of the observed radiation is located at a declination of 
either 26° or 36°. For a point source we would then expect a dip in the cor- 
relation function R at the symmetry axis, i.e. at 6 - 31°. However, Figure 27 
indicates a broad central peak rather than a dip at 6 - 31° . i 
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The scan in right ascension yields the curves shown in the lower half of 
Figure 27. The correlation function for <5= +28 shows a broad central peak 
at a - 13^00 m . The appropriate curve for 6 - +31° also exhibits a narrow 
peak at the same right ascension, indicating that some X-radiation may come 
from this direction of the celestial sphere. 

We have previously suggested (Boldt et al. , 1966) that X-radiation was 
observed in Experiment 1 from an extended region centered at a right ascension 
of approximately 13^00 m and a declination of approximately 28°. The Coma 
cluster of galaxies is an extended object at right ascension a= 12^56 m ±8 m and 
at declination <5 = +28° ±2° (Abell, 1965). The direction to the center of the 
Coma cluster of galaxies is approximately 3° away from the direction to the 
north galactic pole. 

The positional evidence presented above and in Figure 28 allows several 
interpretations: (a) The observed feature in the counting- rate profile of 
Experiment 1 is due to an extended source of X-radiation in the vicinity of the 
north galactic pole, (b) We have observed radiation from a point source with 
varying intensity of X-ray emission, (c) The observed feature is due to a 
slowly varying point X-ray source together with a changing background counting 
rate, due to the variation in the local rigidity cutoff. 

The available information is insufficient for a discussion of the third 
suggestion. The fact that subsequent observations (see below) failed to detect 
a significant flux of X-radiation from the vicinity of the north galactic pole is 
relevant to the first two models. It is improbable that a strong flux of X- 
radiation could have been generated in an extended region of the sky with 


86 



dimensions of a few degrees, over a period of time of the order of one year or 
less (see Table V). This argument would rule out model (a) and leave us with 
the second interpretation of our observations. However, a correlation of our 
observations with possible models for the time dependence of the X-ray emission 
from a point source does not seem to be warranted because of the lack of suf- 
ficient statistics for such a detailed hypothesis. 

From the proportional counter data during the 30-minute period of 
maximum exposure to the north galactic pole we determine the net spectrum 

of events in several energy intervals. Figure 28 shows the observed spec- 

2 

trum for a total exposure of 9000 cm -sec; the width of each energy inter- 
val is 5. 5 keV. An extrapolation to the top of the atmosphere according to 
the data analysis procedure 2 (chapter V) results in a net photon flux 

J(E = 25 keV) « 0. 014 photons/cm 2 -sec-keV. This result, together with 

those of other experiments, will be discussed in section VI E 5. 

2. Experiment 6 . Experiment 6 was launched on 2 Sept. 1966. Six 

importance 1-flares (ESSA, 1967) occurred while the balloon floated at a 

. 2 

residual atmospheric depth of 2. 9 g/cm . This experiment used the "rotating” 
scanning method to examine the vicinity of the north galactic pole. The X-ray 
telescope was programmed to change its zenith angle continuously so that it 
scanned the north galactic pole during every rotation. 

Figure 29 shows a superposition of epochs phased in azimuth angle 
relative to the direction to the north galactic pole. The upper curve shows a 
superposition of the Csl scintillation crystal data from 22 to 60 keV for 30 
rotations of the gondola. The center curve represents the proportional counter 
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data in the energy interval from 12 to 28 keV during the same exposure. The 
lowest curve is proportional to the exposure to a point source along the 
direction of the north galactic pole. 

The mean counting rate for the half rotations opposite iii azimuth to the 
north galactic pole was used as the background level in the determination of 
the standard deviation a . The 3cr level was assumed to be an indication of the 
limiting signal for exposure to the vicinity of the galactic pole. Figure 30 
shows the results of the determination of the 3u limiting signals for several 
energy bands in both counters. In each of the 7 energy intervals, the observed 
rate of events at azimuth angles corresponding to the vicinity of the north 

galactic pole deviated from the background rate by less than 2 a. The 
extrapolation to the top of the atmosphere yields an upper limit to the photon 

flux of J(E “ 25 keV) & 0. 0072 photons/cm -sec-keV. This upper limit applies 
to point sources at celestial positions near the direction to the north galactic 
pole, for example, to the Coma cluster of galaxies. 

3. Experiment 7 . Experiment 7 was intended to duplicate Experiment 1. 
Due to a shift of the Phos witch threshold an excessively high number of events 
were accepted as valid scintillation crystal counts. Because of this, the 
dead-time correction for the proportional counter amounted to 60 to 77% during 
this observation, in contrast to approximately 15% during most of the other 
experiments described here, thus reducing the statistical significance of the 
observations from the proportional counter. The Csl crystal data were 
not analyzed. 

Figure 31 shows the counting-rate profile of the krypton proportional 
counter as a function of time, including the dead- time correction. Events in 
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the energy region from 12 to 28 keV are represented by the upper curve. The 
center curve shows events from a sub-sample of energy for the band 22-28 keV. 
The lower curve is proportional to the exposure to a point source at the 
direction of the north galactic pole. 

The balloon reached ceiling altitude at approximately 1230 UT. The 
counting-rate profile was essentially flat from 1300 to 1900 UT. The observed 
rate was somewhat lower after the meridian transit of the galactic pole, but 
this feature can be accounted for by a loss of altitude which occurred at that 
time. The average background before the transit of the galactic pole was 

used to construct the 3 a lines shown in Figure 31. The total counting rate for 
the 12-28 keV band shows no statistically significant deviation from the back- 
ground. The largest deviation for a sub-sample occurs in the 22-28 keV band in 
which one half-hour interval exhibited a deviation from the background of almost 
3 awhile the galactic pole was still approximately half an hour away. Since this 
feature was not followed by further periods of increased counting rate, it must 
be assumed that it was not the response to a source in the vicinity of the north 
galactic pole. 

During the period from 1600 to 1630 UT, the counting rate could not be 
determined due to difficulties experienced in the retrieval of data from the 
magnetic recording tape. 

Figure 32 shows the 3a upper limits for three energy bands. The energy 
intervals were chosen to be identical to those used in the published analysis of 
Experiment 1 (Boldt et at* , 1966). We assume that the 3 a level determines 
the upper limit to the photon flux from point sources in the vicinity of the north* 



galactic pole at the time of maximum exposure. Except for a single energy 

interval which is centered at 20 keV the upper limits are below the net values 

found in Experiment 1. The extrapolation to the top of the atmosohere results 

2 

in a photon flux of J(E = 25 keV) < 0. 005 photons/cm -sec-keV. 

Towards the end of this experiment the detector’s field of view included 
the constellation Cygnus. The extraction of a net spectrum for the X-ray 
source Cyg XR-1 was made difficult by the large dead-time correction and by 
a change in the background counting rate. However, the resulting spectrum of 
Cyg XR-1 showed sufficient agreement with the spectra of Experiment 4 and 
with the results of other observers to lend support to the validity of the upper 
limits established in this experiment for the vicinity of the galactic 
pole. 

4. Summary of experimental information . In addition to the results of 
Experiments 1, 6 and 7, upper limits to the photon flux from the vicinity of 
the north galactic pole have been given by other observers. During a rocket 
experiment on 25 Nov. 1964, Friedman and Byram (1967a) found that no sig- 
nal with an intensity as great as 7.4% of th*it from the Crab Nebula was ob- 
served in the 1 - 10A wavelength interval from point objects near the north 
galactic pole, e. g. at the position of the Coma cluster of galaxies. If the 
spectra of the objects are assumed to follow power-laws with the same index 
as that of the Crab Nebula, 

c , J(E)«E" 2 , (85) 

then the extrapolated upper limit to the photon flux at the top of the atmosphere 

-4 2 

is J(E = 25 keV) < 9 *10 photons/cm -sec-keV. In an observation on 13 Feb. 

- ■■■ ^ , . . 
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1967, Lewin et al. (i968a) placed an upper limit to the flux from the same 

region of the celestial sphere. From the reported information we deduce that 

J(E = 25 keV) < 15 • 10“ 4 photons/cm 2 -sec-keV. 

The available experimental information on the photon flux from the 

vicinity of the north galactic pole, for example from the direction of the 

Coma cluster of galaxies, is summarized in Table IX. 

5. Discussion. The observation of the Coma cluster as an extended 

scarce of X-rays could be relevant to the interpretation of the measured diffuse 

cosmic X-ray background dj/dJ2. Several models have been proposed for this 

background; one of them assumes that the observed flux is due to the sum 

of many unresolved galaxies. If<J > is the mean emission from a galaxy, 

g 

then the observed flux due to all unresolved galaxies up to a radius 

R « R TT , . . can be estimated from 
Hubble 


dj n < J > R 

g g , (86) 

dfZ 4 tt 

"75 -3 

where n is the mean density of galaxies (n « 10 cm ). When <J > is 
S g g 

set equal to the observed emission of our galaxy, then the resulting value 
of dj/dfi is about two orders of magnitude lower than the measured diffuse 
X-ray flux (Oda, 1965; and Gould and Burbidge, 1967). However, the assump- 
tion that our galaxy is typical with respect to its X-ray emission is somewhat 


arbitrary. If we assume that X-radiation from the Coma cluster of galaxies 

was observed in experiment 1, then we obtain a better estimate for < J > 

g 

since the Coma cluster contains about 1500 galaxies at a mean distance of 
90 Mpc {Abell, 1965). Felten et al. (1966b) have, however, pointed out that 
the resulting diffuse X-ray flux would be too high by a factor of 30 to account 
for the observed background spectrum. 
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The fact that subsequent observations failed to detect comparable X- 
radiation from the Coma cluster of galaxies has somewhat removed this dis- 
crepancy, but has not solved the problem. The upper limits of Friedman and 
Byram (1967a) and Lewin et al. (1968a) reduce the discrepancy factor from 
30 to approximately 3. Considering the uncertainty in the number of galaxies 

in the Coma cluster , the various parameters may actually combine to yield 
the observed diffuse X-ray flux. In this case the "unresolved galaxies 

model for the cosmic X-ray background is not yet ruled out by the limits 
on the X-ray emission from galaxies of the Coma cluster. The conclusion 
would then be that the mean emission per galaxy is approximately two 
orders of magnitude higher than that of our own galaxy. In this connection 
it is interesting to note that, according to recent measurements (Bradt 
et al. , 1967),, the X-ray luminosity of the radio galaxy M 87 is two to 
three orders of magnitude higher than that of our galaxy (Friedman et al. , 
1967b>. 

In view of the fact that other experiments failed to observe a 
significant photon flux from the vicinity of the north galactic pole, we may 
tentatively attribute the observations of 6 December 1965 to a transient 
phenomenon associated with a discrete object. 

No information from optical observations during the time of our first 
experiment is available. Any unusual activity in the vicinity of the north 
galactic pole would probably have remained undetected since this region of the 
celestial sphere was in a very unfavorable position for reliable optical 
observations. 
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Four flare stars were observed during 1964 and 1965 in a region in Coma 
Berenices, approximately 6° away from the direction of the north galactic pole 
(Haro, 1966). The number of flares that were observed per hour in that region 
was 0. 05. We are not aware of flare observations during the time when Ex- 
periment 1 was conducted. 

Pencil-beam radio observations of the Coma cluster of galaxies (Bozyan, 

1968) have recently resolved a small-diameter source, Coma A, which is 

centered at a = 12^51. 8 m , 6 = 27 ° 53. 8 f . The energy spectrum of Coma A 

has an index of T~ 0. 6. An extrapolation of the power -law radio spectrum 

can account for the optical luminosity (Matthews and Sandage, 1963; Schmidt, 

-4 2 

1965; Wyndham, 1966) and yields a photon flux of 2. 5 . 10 photons/cm -sec- 
keV at E = 25 keV. Table IX shows that this flux is approximately two orders of 
magnitude below the previously reported value (Boldt et al. , 1966) and about 
one order of magnitude below the upper limits found in other experiments. 

It is of interest to discuss the hypothesis that the observed feature in 
Experiment 1 was due to a supernova explosion in the Coma cluster of galaxies. 
The duration of X-ray emission during a supernova outburst has not yet 

been established. Colgate's hydrodynamic model (Colgate, 1968) suggests a 

-3 

time-span for X-ray emission of the order of 10 sec. The fluorescence 
theory of optical emission from a supernova (Sartori and Morrison, 1967b) 
assumes that the explosion takes place within a few seconds. If the 
above hypothesis were correct, the observations on 6 Dec. 1965 would in- 
dicate a period of at least one hour. The photographic light curve for a 
type I supernova suggests a time- scale of about 20 days for the main phase. 
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The absence of a significant flux of soft X-ray emission within about 30 days 
after a supernova explosion (Bradt et al. , 1968b) indicates that an upper 
limit of about 30 days can be set on the duration of X-ray emission. The 
nova-like behavior of the low-energy X-ray emission from the variable 
source Cen XR-2 (Harries et al. , 1967; Cooke et al. , 1967; and Chodil et 

al. , 1967) indicates an e-folding time of about 30 days. 

With the assumption that the Coma cluster contains 1500 galaxies, about 

4 supernova explosions per year are expected when a rate of 1 event per 400 

years per galaxy is adopted. Considering the uncertainty of the observations, 

this estimate is supported by the fact that 12 supernovae have been observed 

in the Coma cluster during a recent 10 year period (Zwicky, 1967). The 

probability of observing the main phase of a supernova explosion during a 

one-hour observation is therefore 20% if the main phane lasts 20 days, and 

is 0. 04% if it is assumed to last 1 hour. 

We can also estimate the energy output in order to evaluate the plausibility 

of the supernova hypothesis. We assume a power-law photon spectrum of index 

2 and compute the energy flux in the X-ray interval from 1 to 100 keV. If 

the distance to the source is 90 Mpc, then we arrive at a total X-ray energy 

51 

output of q = 8 • 10 ergs emitted over a period of 20 days. This is not 

unreasonably high. In fact, the fluorescence model (Sartori and Morrison, 

52 

1967b) requires an output of 10 ergs in the ultraviolet band. If the 

emission is assumed to last only one hour, however, then we find an output 
49 

of q = 2 • 10 ergs. The X-ray energy output increases by an order of 

7 

magnitude if thin- source bremsstrahlung at T = 5 • 10 °K is assumed. 
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Colgate’s model (Colgate, 1968) predicts an X-ray pulse with a total energy 

47 

output of only q ^ 5 • 10 ergs within a fraction of a second and can there- 


fore not be invoked to explain the observation of experiment 1. 



CHAPTER VH 


SUMMARY 


A series of balloon-borne X-ray experiments was carried out between 
December 1965 and July 1967. The sensing elements of the X-ray telescope con- 
sisted of a Krypton proportional counter and a CsI(Tl) scintillation crystal with 
an effective energy range from 12 to 60 keV. The full- width angle of the detector 
aperture at half maximum was approximately 24 In the course of the balloon 
observations, a broad band of the celestial sphere about a declination of +30°, 
the vicinity of the galactic center, and the Constellation Scorpius were surveyed. 

A study of the data from several balloon flights, showed that there exists 
a definite relation between the background counting rate and the geomagnetic 
latitude. After subtracting the contribution from the diffuse cosmic X-raj^ back- 

i 

ground, we are left with a counting rate which is largely caused by the atmos- 
pheric X- and gamma- ray background. It is assumed that most of this photon 
flux is a component of the atmospheric electromagnetic cascade and is therefore { 

linked to the flux of charged particles incident at the top of the atmosphere; how- 
ever, a complete calculation of this process has not been carried out to date. A 
correct prediction of the geomagnetic-latitude dependence of this flux rate con- 

- . . . • t- 

] 

stitutes a stringent test for any model of the production of atmospheric X- \ 

radiation. J 

■ - ■ . . . ■ ■ i 
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Two observations of the Crab Nebula (Tau XR-1) and three observations of 

• .. ' ■ ■ ....... ■ ■ . • • . - i 

the X-ray source Cyg XR-1 resulted in net source spectra in substantial 
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i 
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agreement with the results of other experimenters. Our results and the observa- 
tions of other experimenters at lower and higher energies are consistent with a 
power-law spectrum of index n « 2. 0 for the photon flux from Tau XR-1. One 
exposure to Tau XR-1 with high statistical accuracy was used to demonstrate two 
distinctly different methods of data analysis and to show that the resulting sets 
of spectral parameters are mutually consistent, although one of the methods is 
superior because it involves fewer approximations. 

The measurement of Cyg XR-1 with the highest statistical significance 
yields a spectrum that is similar to that of Tau XR-1 in shape and in the absolute 
flux level. With regard to the reported variability of Cyg XR-1, we find that our 
results are consistent with the spectrum derived by Bleeker et al. (1967) from an 
experiment that was conducted three weeks prior to the measurement reported 
here. 

A meridian-transit scan of Sco X-l was carried out during a balloon experi- 
ment over Australia. The data analysis for this flight was somewhat difficult 
because of the variation of the residual atmospheric depth during the exposure to 
Sco X-l. At energies between 15 and 35 keV, the resulting source spectrum 
was consistent with thin-source bremsstrahlung at a temperature of approxi- 

7 

mately 5 • 10 °K. However, the photon flux at higher energies was found to 
be distinctly above the extrapolation of the low-energy spectrum. A similar 
spectrum was observed in June 1965 by Peterson and Jacobson (1966) . The 
comparison with the results of other balloon experiments shows that both the 
shape of the spectrum and the absolute flux above approximate! y 15 keV vary in 


time. The high-energy feature that was observed in this experiment can be 
understood qualitatively in terms of both Shklovsky* s (1967) stellar accretion 
hypothesis and the bremsstrahlung from an iso velocity plasma as suggested by 
Fischel and Stecher (1968). Any detailed model of Sco X-l must take into 
account the existence of a variable high-energy component. 

An observation of the vicinity of the galactic center detected a significant 
flux of hard X-radiation. The resulting spectrum appears to be somewhat 
softer than that of the Crab Nebula. Subsequent experiments confirmed our 
finding, but they cannot be directly compared with our observation because of 
the difference in the solid angles of the various detector apertures. We conclude 
from these observations and from a calculation of the expected equivalent flux 
for rocket-borne measurements that we have observed X-radiation from more 
than one source in the vicinity of the galactic center , perhaps with different 
degrees of spectral, hardness. Further detailed balloon-borne observations will 
be required in order to resolve the hard X-radiation from the many sources in 
the vicinity of the galactic center. 

Three observations of the vicinity of the north galactic pole brought con- 
tradictory results. An increase in the counting rate during the first experiment 
was interpreted as the X-ray flux from an extended source in the direction of the 
Coma cluster of galaxies, which is a few degrees away from the direction of the 
north galactic pole. However, subsequent observations failed to detect a simi- 
lar feature and set upper limits to the photon flux from celestial objects near 
the direction of the north galactic pole, including the Coma cluster of galaxies. 
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It is concluded that if the observed response was caused by a cosmic source of 
X-radiation, then the source intensity must have varied in time. The hypothesis 
that a supernova explosion at the distance of the Coma cluster of galaxies caused 
the signal detected in the first experiment is shown to lead to an energy output in 
the X-ray region which is not excessively large for such a phenomenon. 

In the present series of experiments, we have measured four of the 
five cosmic X-ray sources or source regions which have to date been clearly 
observed from balloons (the fifth source is located in Centaurus, a region 
which was not scanned). The detector that was used for the measurements 
was designed for an exploratory survey of the celestial sphere. The telescope 
aperture had a rather large opening angle so that a considerable counting rate 
was produced by the atmospheric and cosmic diffuse X-ray background. This 
contribution amounted to approximately two-thirds of the ceiling counting rate, 
so that only about one-third of this rate represented the inherent detector 
background. The major portion of the detector background is caused by Compton 
interactions of gamma radiation in the detector and the surrounding material. 

The remaining event rate, a small fraction of the total counting rate, is proba- 
bly caused by an incomplete rejection of charged particles and other effects. 
Closing the detector aperture without increasing the inherent background counting 
rate would have reduced the total event rate at ceiling by a factor of 3 to approxi- 
mately 13 * 10 counts/cm -sec-keV. This is in contrast to the best presently 
available scintillation crystal detector telescope which has an inherent detector 

_4 2 

background level of 3 • 10 counts/cm -sec-keV (Peterson et al. , 1968a) at a 
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residual atmospheric depth of 3. 4 g/cm and a geomagnetic latitude of 40° . The 

gamma-ray induced inherent background can be reduced by a factor of 2 to 3 if 

the experiment is conducted near the geomagnetic equator. 

It appears to be feasible to build a proportional chamber with the same 

-4 2 

background event rate (3 ♦ 10 counts/cm -sec-keV), an efficiency of typically 

2 

GO percent, and an exposure area of the order of 1000 cm . With net exposure 
times during balloon experiments of approximately 3 hours, it will then be pos- 
sible to detect hard X-radiation from several other sources which have been ob- 

* 

served in rocket-borne experiments, provided the low-energy spectra can be 
extrapolated to higher energies. The available pointing accuracy for balloon- 
borne experiments is already sufficient to allow the reduction of the field of 
view to a value where the resulting counting rate caused by the diffuse atmos- 
pheric and cosmic background is comparable to the lowest attainable inherent 
detector background. The inherent and diffuse-background counting rates of the 

above-mentioned scintillation detector are equal for a solid angle of the field of 

-2 

view of approximately 2 * 10 ster, which corresponds to a circular aperture 
with 6° full width opening angle at half maximum. This demonstrates that the 
total counting rate cannot be significantly reduced by making the aperture angle 
smaller than about 4 ° . 

It is expected that X-ra y astronomy measurements will evolve along two 
somewhat overlapping lines of development. In the first, one will try to obtain 
information on the spatial distribution and the types of continuum spectra of 
X-ray sources and their correlation with optical and radio objects. Surveys of 



the celestial sphere with improved sensitivity and precise position determina- 
tions are best carried out at low X-ray energies at the top of the atmosphere. 
Balloon or higher-energy satellite observations with high angular resolution and 
good pointing accuracy during extended exposure times can then be used to 
determine the high-energy continuum spectra of these X-ray sources. 

In a second line of development, one may attempt to resolve lines in the 
X-ray spectrum. The most prominent low-energy lines are generated in 
stellar atmospheres and may expand our knowledge of the physical processes 
occurring there. However, a unique aspect of X-ray astronomy lies in the fact 
that some intense emission lines in this energy range can be produced in nu- 
clear processes and could provide first-hand information on energy-releasing 
mechanisms. For example, the 60 keV line of Am 241 is expected to be gener- 
ated during the decay of isotopes that were produced via the r-process in super- 
nova explosions. The intensity of this line has been estimated for the case of 
the Crab Nebula (Jacobson, 1968). The above-mentioned proportional chamber 
could detect a So- line signal above the continuum in a 5 keV-bin if its energy 
resolution is 10 percent at 60 keV, and if we assume a 30-hour exposure to the 
Crab Nebula from the top of the atmosphere. This would require that satellites 
be pointed accurately for extended periods of time instead of used for sky 
patrols at X-ray energies. Such a shift in the mode of operation of X-ray 
satellite experiments would aid both wide-band continuum and line measurements. 

Balloon-borne experiments will continue to play an important role in the 
field of X-ray astronomy. While rocket-borne experiments are restricted both 
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in exposure time and maximum available detector area, and satellite experi- 
ments are limited in weight and availability, balloon observations are more 
favorable in all of these respects. 

The continuing development of new detectors is expected to reduce the 
inherent detector background and improve the energy resolution. However, for 
a given detector system, optimum results can be obtained through a detailed 
study of the behavior of the detector in the radiation environment near the top of 
the atmosphere* Knowledge of the detector event rate as a function of atmos- 
pheric depth, geomagnetic latitude, and energy will reduce systematic errors 
in the extraction of net source counting rates, particularly for exposures to 
variable sources. It will also make it possible to use a large fraction of the 
flight time at ceiling for the exp Dsure to a source or source region. 

Variable sources require simultaneous observations over a large energy 
region. As the recent studies of Sco X-l have indicated, programs for optical, 
ultraviolet, and low- and high-energy X-ray measurements must be set up 
with the emphasis on both line and continuum observations. It is hoped that the 
resulting data will supply the Information necessary for the construction of con- 
sistent models for discrete X-ray sources. 


APPENDIX A 


THE ESCAPE FACTOR 

A photon of energy E incident on a detector material can be absorbed 
through the photoelectric effect. The emitted electron will come from the K 
shell with probability a , from the L shell with probability a r , etc. If . he 

& Jj 

relative probabilities for K, L^, M^, etc. shells or subshells are related to 

3 3 3 

the azimuthal quantum numbers 1 as l/l , i.e. (1, 1/2 , 1/3 , etc.), then 
a = 0.86, a = 0.106, a = 0.031 , etc. 

1 X 

(Grodstein, 1957; Lay, 1934). The remaining atom can be de-excited via the 

Auger effect or by fluorescence. The latter occurs with probability /3(Z) 

which is a strong function of Z. The fluorescence radiation (in practice only 

the characteristic K-, L-, etc. radiations with energies E , E , etc. need 

xv L 

to be considered) can escape from the detector; in this case the amount of 

energy deposited in the detector is equal to E - E_ _ , and forms the 

K, (L. . . 

"escape peak. " 

The escape factor 77 (E) is defined as the fraction of the absorbed number’ 
of photons which escapes in the form of fluorescence photons. For an infinite 
plane-parallel slab of thickness x q , and escape of characteristic radiation E^ 
from both surfaces, we have 
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V (E, x ) 


27 r 1 


S exp(- /z(E)x) n(E)a/3 


x. 


4 7 r / exp(-^(E)x) p. (E) dx 
o 


(87) 


I f [ exp • ) + exp ( — ] d cos 0 d <p dx 

% a 1 * ' COS v / ' nno A T 


o o 


cos 0 


4 7r / exp(- p(E)x) ^z (E) dx 
o 

where we have assumed that the incident photon arrives normal to the 
entrance surface. 

With the definition for the exponential integral (Abramowitz and Stegun, 1964) 


E X W = 7 dz , 


( 88 ) 


the escape factor can be re-written as 

x. 


>1 (E, x Q ) = 1/2 a/3 


/ exp(-p(E)x) [ E 1 (n(E c )x) + Ej (M E c ) (x q -x) ] dx 


1 - exp (- p(E)x 0 ) 


(89) 


The equation for escape from the entrance surface only is 

x o 

H (E, xj = 1/2 a /3 o ex P(- M(E)x) E i( m(E c )x) dx 

1 - exp<-p(E)x 0 ) 


(90) 


Both equations require numerical evaluation. For a semi-infinite medium, 
x o -*°o, equation 87 can be integrated directly: 


^< E c> 


v < E ’°°> ■ 1/2 1 *~ m ~ ln (1 * m ) } 1 * 


(91) 


We apply these equations to the mm of Csl(Tl). Table X shows the 
characteristic line energies (BIoki& mzd relative intensities 




-■ ^StSF"**?** 


(interpolated and extrapolated from data given by Williams, 1933) for Cs, I, 
md Kr. Only the K shell is considered in table X. 


104 


hi order to simplify the calculations, the CsI(Tl) crystal is assumed to 

emit only one fluorescence line at E.^ = 30 keV; the mean line energy, weighted 

with the relative intensities, is 30. 5 keV. The ,K- fluorescence yield is taken 

to be /3 = 0. 88 (Bergstrom and Nordling, 1965; Fink et al. , 1966). The 

probability of creating an L-shell vacancy and the fluorescence yield for the 

L shell ( (3 =0.1) are considerably smaller than those for the K shell 
L 

OLi /?, 

(- 7 - j—n — = 0. 014); this justifies the neglect of L shell fluorescence. 

°K P K 

2 2 

Figure 6 shows 77 (E, x Q ) as a function of E for x Q = 0. 1 g/cm , x q = 0.4 g/cm 
and x q -* 00 . The escape factors for both surfaces are drawn as solid lines 
(compare Stein and Lewin (1967) for Nal). A broken line represents the func- 
tion 77 (E) for escape from the entrance surface only for a crystal of thickness 
x q 0 . 1 g/cm . The latter type of function for a crystal depth of x q =0.4 g/cm 2 
is almost identical with 77 for the semi- infinite detector, where both cases ob- 
viously result in the same escape factor. Escape from the entrance surface is 
important since the principle K X-radiation of Cs and I can contribute to the 
energy region of interest in the proportional counter (section III E). 

The curves for escape through one or both surfaces of the Csl crystals 

2 

which were used in our experiments (2 mm Csl or 0. 92 g/cm Csl) are prac- 


tically identical to the curve for a semi -infinite medium ( x q — 00 ). The cor- 
rection for the finite lateral extent of the crystals is negligible. 


APPENDIX B 


v.. . r^r-nu... 


PRODUCTION OF X-RADIATION NEAR THE TOP OF THE ATMOSPHERE 

Figure 10 shows the sum of the detector-generated, plus the atmospheric, 
X-ray background as a function of the local vertical cutoff rigidity. Most of the 
events shown here are due to the atmospheric X- and gamma radiation; a small 
fraction is due to corpuscular radiation and other effects (section IH G). 

It is tempting to draw a curve showing the "pure atmospheric X-ray 
background" versus the cutoff rigidity, and to compare it with other particle 
intensity curves, for example, for the proton, neutral pion and neutron fluxes. 
However, in the rigidity interval from 4 GV to approximately 5. 3 GV the slopes 
of the appropriate curves for energetic neutrons, protons, and pions are roughly 
equal to each other and to such a "background X-radiation" curve, so that no 
definite statement as to the origin of atmospheric X-radiation could be made 
on the basis of such a correlation alone. 

We will now present a short discussion of some mechanisms for the 
production of atmospheric X-radiation. It is assumed that the bulk of the 
observed photon flux is caused by the electromagnetic cascade in the atmos- 
phere. Other mechanisms, discussed below, can perhaps account for as much 
as 20% of the observed X-ray flux at balloon altitudes and at A «40°. Gamma- 
ray emission from thermal-neutron capture in 0 16 and N 14 is assumed to make 
a negligible contribution to the atmospheric X-, and gamma-ray flux 
(Peterson, 1963). Bremsstrahlung from re-entrant electrons has not been 
considered here because oi. the lack of relevant experimental information. 
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Internal or "inverse" bremsstrahlung (Hayakawa and Matsuoka, 1964b; 
Boldt and Serlemitsos, 1968a), during collisions of primary and secondary 
cosmic rays with atomic electrons, makes a small contribution to the at- 
mospheric X-radiation. At approximately 30 keV the inverse bremsstrahlung 

o —4: 

flux at balloon altitudes, A «40 , is less than 5 * 10 of the observed atmos- 
pheric X-ray flux (Peterson et al. , 1967a). 

Unless otherwise stated, the spectral flux J will, in this appendix, be 
measured in particles/m 2 -sec-MeV; the differential flux per unit solid angle, 


N = , 
d£2 


(92) 


will be measured in particles/m -sec-st^r-MeV, and energies are measured 


in MOV. The production flux 
dN d 2 J 


Q 


dz dz d£2 


(93) 


2 2 

has the dimension particles/m -sec-ster-MeV-g/cm , unless otherwise noted. 

I. Compton- Scattered Photons from Primary X-rays 
The spectrum of the primary cosmic X-ray background has been measured 
by several observers (Klarman, 1962; Metzger et al. , 1964; Garmire and 
Kraushaar, 1965; Brini et al. , 1965b; Rothenflug et aL. , 1966; Kohn et al. , 

1966; Peterson et al. , 1967a; Peterson, 1967b; Bleeker et al. , 1968) in balloon, 
rocket, and satellite experiments. Balloon-borne observations in the 10 - 100 
keV energy range must be extrapolated to the top of the atmosphere; the re- 
sulting spectra fit reasonably Well with other observations at lower and higher 
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energies, where this extrapolation is not necessary (cf. the summary given 
by Gould, 1967a). Thus, we may assume the spectrum of the primary cosmic 
X-ray background to be known. 

We are interested in the spectrum of secondary photons created through 
Compton scattering of the primary cosmic X-rays by atmospheric atomic elec- 
trons. If we limit ourselves to one interaction into an angle within 1 ster 
of the forward direction, then the differential production spectrum Q(E , z), 

X 


deriving from a primary spectrum N(E , z) is 

X 


Q(B, z) = N q \ / 


E (E , © ) 

■w» nv !| v ' r% v' 


max’ x max 


d0 


x 


E 


N(E x , z) ~ (E, E x )dE. (94) 


x 


For small-angle scattering and E< 100 keV, the interval |e , E (E ,'©)! 

I x max x I 

becomes very small, so that we can approximate 


Q(E x , z> ~N 0 | N(E x , z) J 


E (E , © • 

max x' max do „ v 

(E, E ) dE 
dE x 

E 

x 


(95) 


N (E , z) n (E , © ) . 

' x’ ' r c ' x max' 


The differential equation for the depth-dependence of the secondary radiation, 


s, in a one-dimensional model is 


~ = Q - ( P T - ix J s = m c N (z = 0) exp (-p T z) - (-P T - P c )s, (96) 
where fx^ is the total (photoelectric and Compton) absorption coefficient. The 
solution 

s(z) = N(z = 0) jexp (-( fXj- n c )z) - exp (- n T z)J (97) 

has a maximum at 


z 

max 







T “ 



(98) 
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For example, at E - 30 keV; z = 3. 27 g/cm we obtain the maximum 

' max ° 

secondary flux 

s(Zmax) " 0, 0013 photons/cm 2 -sec-ster-keV. (99) 

This is to be compared with an experimental observation of the atmospheric 
X-ray background flux at z = 3.5 g/cm , A = 40° (Peterson et al. , 1967a): 

N(E = 30 keV) »0. 01 photons/cm 2 -sec-ster-keV. (100) 

This production mechanism is independent of latitude, provided the primary 
cosmic X-ray background is isotropic. An upper limit to the anisotropy of 
the diffuse X-ray background was set by Seward et al. (1967), at approximately 
10 %. 

n. Bremsstrahlung from Precipitating Electrons in the Atmosphere 
Energetic electrons caa be trapped in the earth’s magnetic field. At- 
mospheric absorption of electrons (Cladis and Dessler, 1961; Walt, 1964) is 
an important loss mechanism; we are here interested in the bremsstrahlung 
produced in this process. 

Cladis and Dessler (1961) have evaluated the precipitation flux of elec- 
trons over the Capetown region. With an experimentally determined pitch-angle 
and energy distribution of radiation-belt electrons, they determined the elec- 
tron flux entering the atmosphere over the Capetown anomaly, where the mirror 
altitude of electrons is very low. We can compute the bremsstrahlung spectrum 
in a production layer which is assumed to be thin compared to the atmospheric 

depths at which balloon-borne detectors float. The resulting photon flux which 

2 

would arrive at a residual depth of 3. 5 g/cm oyer the Capetown anomaly was 
compared with one-half of the integral flux into 47r which was measured by 


Peterson et al. (1967a) at this atmospheric depth and A = 40°. The ratio 
of these values is 0. 02 at 30 keV and has a maximum of 0. 034 at approxi- 
mately 100 keV, 
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During magnetic disturbances the field lines can be altered so that the 
mirrow points are lowered deeper into the atmosphere; an increased flux of 
bremsstrahlung radiation would then be expected. Ghielmetti et al, (1964a), 
reported the observation of an enhancement of the X-ray intensity at balloon 
altitudes over the South American anomaly. However, as Hudson et al. (1966), 
have pointed out, the X-ray intensity observed in this experiment can be ex- 
plained as being due to the diffuse cosmic X-radiation. Therefore, ex- 
istence of X-radiation from precipitating trapped electrons to a degree 
comparable with the observed atmospheric X-radiation is not confirmed. 

in, Bremsstrahlung from Electrons 
The electron bremsstrahlung cross-section can be approximated by the 
function 


(E , E ) = AE a • E" b 
dE e x' e x 

x 


( 101 ) 


where E is the electron energy and E x is the photon energy. If we assume a 


power-law spectrum for the electrons, 

J (E ) = BE" C , 
e v e e 


(102) 


then the resulting photon spectrum is 


i ( g«r 


no 


W - I J e < E > W < E ’ E x )dE 
E x 

E. 


(103) 


- / /2 J(E)H-(E, E x )dE+ / J (E)|| (E, E x )dE , 

E x x Bi/a x 

where the energy E.^ that raakes the two integrals on the right 

side of the last equation equal to each other. In other words , one-half of 

the photon flux at a particular energy E is due to electrons with E < E ]/£, 

and one-half is due to electrons with E g > E^yg, The ratio 


E 


E 


1/2 .In 2 t 

= exp ( r- ) , 

' c-a-1 ' 


x 


(104) 


is then a good indication as to which part of the electron spectrum is relevant 
for the bremsstrahlung flux. The equation holds only for c-a-1 > 0. 

A study of the bremsstrahlung cross-section as given by Bethe and Ashkin 

j-vw *' 111 flirt** env»pAVtirt^ dimuc fliof fAV pvqmnlp 1? = 0 0^- Qnrl 

iJ UiJ i} uiviliuliig puwirw vj.iw.1 Ivi C AUjuiyiC j ava ■ — vi v j. ItAv V 

X 

0. 16 < E <40 MeV, the exponent, a, becomes «0.2. For C - 1. 5, E., /o / 
e ~ x/ & 

E - 10, and for c = 2. 5, E„ E =1.7. This demonstrates that electrons of 
x 1/2 x 

energies below a few MeV are of greatest importance in the production of X- 
radiation, provided the election spectrum is of the steepness indicated above. 

There exist only a few measurements of the atmospheric electron compo- 
nent below 10 MeV; they are summarized in Figure 33. The data points of 

Brunstein (1967) and of Cline and Brunstein (1966) have been obtained at at~ 

2 

mospheric depths of 5, 5 and 4 g/cm , respectively. The data points of 

2 

Beedle and Webber (1968) are the result of an extrapolation to A g/cm . 

These measurements were made at the latitudes of Minnesota and further north. 


** ii 1 ^/ " 1 



The atmospheric electron flux is related to the cosmic-ray flux, and the lati- 
tude dependence is assumed to depend roughly on the integral cosmic-ray flux. 
It is therefore necessary to compute the atmospheric electron spectrum at the 
latitudes where measurements are available, and to compare the calculations 
with the experimental results. The spectrum can then be extrapolated to or 
re-calculated for other latitudes. 

The calculation can be divided into four sections: 

1. The electron production spectrum 

2. The electron equilibrium spectrum 

3. Bremsstrahlung of electrons (with the calculated or experimental 
spectrum) 

4. The equilibrium photon spectrum (e.g. , at balloon altitudes), 

1. The electron production spectrum. Low-energy electrons are pro- 


duced in the atmosphere by several processes. 


We will discuss the four most 


important mechanisms (compare Beedle, 1966). 

(a) Compton collisions of primary and secondary photons with atmospheric 
atomic electrons. For this and for the next two processes we need to know the 
spectrum of primary and secondary photons as a function of the residual at- 
mospheric depth, d, and the geomagnetic latitude. From 0. 05 MeV to approxi- 
mately 10 MeV we use the measurements of Peterson, et al. (1967a) at 

2 

A - 40* and a depth d = 3. 5 g/cm to define the secondary photon spectrum 

-1 7^5 

N(E X ) - 600 E * ' (105) 

The cosmic X-ray background between 0. 03 and 1 MeV at the top of the 
atmosphere is assumed to be represented by (Gould, 1967; Bleeker et al. , 1968) 
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N(E ) * 150 E 

' X X 


*2,5 
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The electron production spectrum is 


«< E e> = l 5b < E e’ E x> N < E x> dE x' 
E . e 
mm 


E , =Ke 

min 2 p5 e\j E 0 J 


(107) 


where Q has the dimensions particles/m -sec-ster-MeV-g/em , and 

(E e> E^) is the cross-section per g of air per unit electron energy for the 
i e 

production of an electron of energy between E and E + dE by an incident 

e e e 

photon of energy E , The relevant curves for electron production due to the 

X 

secondary and primary photon spectra are shown in Figure 34 as curves 
number 2 and 3, respectively. It is apparent that the secondary photons 
(at 3. 5 g/cm ) are more effective than the primary photons (at 0 g/cm ). The 
altitude dependence of the photon spectra and the resulting electron spectra is 
rather complicated; however, we can make the qualitative statement that the 
relative contributions of the two components remain the same over the first 


few g/cm near the top of the atmosphere. 

(b) Production of photoelectrons due to the primary and secondary photon 


£ 

f | 


spectrum in air. The cross-section for the production of electrons via the 
photoelectric effect has been given by Evans (1955). We adopt the approxima- 
tion 


0 <E e , E x ) = 1. 9 • l(f 6 E x " 3 (cm 2 / g) 

= 3.3 ■ 10" 6 2 ~ 1, 5 (cm 2 /g) 

X. 


E < 0. 7 MeV 
x 

E > 0. 7 MeV 
x 


(108) 
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With the approximation that a photon of energy E produces an electron of 
energy E - E we obtain the production spectrum 

6 X 


Q(E e ) = 0 <E e , E x ) N <E x ) . 

The results are shown in Figure 34 as curves 6 and 7. 

(c) Pair production of electrons and positrons in air. 

for pair production was adopted from Evans (1955) 

0. 53 E - 0. 6 

5 

E x - 2 m 

E -2m 
x 

The resulting electron production spectrum is 


(E , E ) = 0.0012 
dE ' e x' 
e 


= 0 . 0012 


(109) 


The cross-section 


£ E <5 MeV 
x 

( 110 ) 

E > 5 MeV 
x 


Q(E e> = l , 51 < E e’ E *> N < E x> dE x < 1U > 

E + 2m e 
e 

The production spectra due to the primary and secondary photon flux are drawn 
schematically in Figure 34 as the curves 5 and 4, respectively. 

(d) Production of electrons by collisions of primary and secondary 
cosmic rays with atmospheric atoms. 

The primary cosmic-ray spectrum was assumed to follow the analytical 
expression given by Balasubrahmanyan et al. , (1967), adjusted to fit the sum- 
mary spectrum given by Gloeckler and Jokipii, (1967), for measurements near 
the solar minimum in 1965: 

N(E p ) = 2 • 10 8 Ep 1 * 5 (E p + 500)” 4 (112) 

(N in particles/m -sec-ster-MeV, E p in MeV). 


114 


Measurements of the secondary cosmic ray spectrum have been carried out 

2 

by Ormes (1967; at 3 g/cm , Ft. Churchill, 1965) and by Teegarden (1967; at 
2 

4 g/cm , Sioux Fails, S. Dak., 1961 and Fairibault, Minn., 1964). In order 
to simplify the calculation we have approximated the secondary cosmic-ray 
spectrum with the following function: 


N(E p ) = 0.65 

N(E p ) = 52/Ep 

N(E ) = 10400/E 
P P 


E < 80 MeV 
P 

80 MeV < E < 200 MeV 
P 

200 MeV < E . 

“ P 


(113) 


At the geomagnetic latitudes of Palestine, Texas the secondary cosmic ray 
spectrum was assumed to be half as intense as at the latitudes of Sioux Falls 
and further north (Teegarden, 1967). The dependence of the secondary proton 
intensity on atmospheric depth was taken from the measurements of Teegarden 


(1967). 

The proton-electron knock-on cross-section is well known (Evans, 1955; 
Rossi, 1952): 


I 00 ' < E e’ V = 2 ’ N f r e “ 
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1 - 


M 


(E + W 


1 “ { 1 “ ( 
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E + M 
P 


_ E {2m E + (M +m) } 
x 2 \ e P 
} J 2m E (E + 2M) 

P P 


E 


(114) 


2 (E + M)' 
P 


where M is the proton mass and m represents the electron mass. Kinetical 

considerations show that a minimum proton energy E . (E ) is required 

* p, mm e 


it- 
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to produce an electron of energy E : 

e 


E . (E)=iE m + 7E 2 +~E **~-E (115) 

p, mm e' 2 e V 2e 4e2em 4m e ' ' 

The production spectra of electrons at different geomagnetic latitudes A 

and atmospheric depths d are obtained from 


On 00 y co l 

Q <E e , A , d) = ^ J - (E e> E ) N (E , A , d) dE , (116) 


e E_ e 
1 


e p' 


where denotes the effective solid angle of the primary and secondary cosmic 
rays and represents the solid angle into which knock-on electrons are 
emitted. We have assumed a sharp geomagnetic cutoff at E (A); E 1 is the 
larger of the two energies E q (A) and E^ Primary particles heavier 

v>an protons were included under the usual assumption (Abraham et al. , 1966) 
that these nuclei behave like ’.’A" nucleons in an interaction, where A is the 
number of nucleons in the primary particle. 

At this point, some comments are in place regarding the relative contri- 
bution of secondary cosmic rays and the effect of multiple collisions. 

A study of equation 114 shows that the knock-on cross-section for fixed 


electron energies of typically 0. 1 - 1 MeV does not vary significantly between 

E . and proton energies of a few GeV. This means that primary cosmic 
p , min 

rays are the dominant producers of knock-on electrons. The primary cosmie- 

2 

ray spectrum is not attenuated significantly within the first 10 g/cm near the 

2 

top of the atmosphere (the collision length for protons in air is 86 g/cm ), so 

that the total production flux of electrons near the top of the atmosphere is, to 
a first approximation, independent of the atmospheric depth. 
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The effect of multiple collisions can be evaluated by considering the next 

step, electron-electron collisions. An estimate of the production spectrum Q f 

due to electron-electron collisions showed that Q f « 0. 4 Q at E =0.1 MeV and 

© 

is less at higher energies. 

The east-west effect has not been evaluated, but it is not assumed to 
change the results significantly. 

In Figure 34 we show the calculated production spectrum (curve 1) at an 

. 2 

atmospheric depth d = 1 g/cm for a proton kinetic -energy cutoff E = 3600 MeV, 

c 

corresponding to the geomagnetic latitude of Palestine, Texas. It is apparent 

that the cosmic-ray collisions provide the largest contribution to the electron 

flux at energies below a few MeV. In Figure 34 we have also plotted (curve 8) 

2 

the electron spectrum which is produced by pion decay at d = 1 g/cm and at a 

cutoff energy ’E = 3600 MeV (Verma, 1966). 
c 

2, The electron equilibrium spectrum . We must now calculate the 
equilibrium spectrum since the range of low-energy electrons is very small 
compared with the atmospheric depths considered. In a one-dimensional model 
one solves a transport equation 

aj , E max (E) d0 col 

(E, z) = Q(E, z) - J(E, z) J Ri52Ji d “, (E, E’) dE' 

d0 < 117 > 

op y Pn l 

+ J J (E + E\ Z) — U1 (E + E', E') dE’ , 
and the solution leads to a straggling equation 

oo d 

J(E , d) = / dE / d2 Q <E, z) G (E, z), 

6 E o 

e 


( 118 ) 


where the Green’s function (Rossi, 1952) is 


G(E, z) 


2 7T p Z 


(e - e y 

' e a 

7 2 

2 p z 


E = E - k (d - z) 

a 


k (E) = / E’ ~ 0 (E, E') dE’ 

o 

pV)= / E ' 2 i 0 eo! (E + E', E') dE- 

° dE 

The resulting equilibrium spectra are shown in Figure 35 for a cutoff energy 

E =20 MeV (curve 1) and for E = 3600 MeV (curve 2). We also give the data 
c c 

points of Cline and Brunstein (1966) and of Brunstein (1967). The straight 
line (curve 3), at energies between 0. 05 MeV and 5 MeV, represents a fit to 
the measurements of Beedle (1966), in an experiment at Ft. Churchill. There 
appears to be a discrepancy of about one order of magnitude between the 
curves 1 and 3 and the measurements of Cline and Brunstein. In the following 
considerations we will use the latitude dependence as previously calculated 
and normalize to the measurements of Cline and Brunstein (1966). 

3. Bremsstrahlung of electrons. The energy range considered here is 
between the energy intervals where the non- relativistic or the extreme 
relativistic limits of the bremsstrahlung cross-section are valid. We have 
used an approximation (Holt and Cline, 1967b) which assumes the appropriate 
values in the high- and low-energy limits. 

4. The observable X-ray flux. The observable X-ray spectrum was 


calculated according to the equation 




1X8 


d © 

max 

J (E x , A , d) = 2 7T / / Q x (E x , A , z) 

o o 

exp [l» air (E x ) i ^-] F (6) sine de dz , (m) 

where f(0) represents the detector angular response function. Figure 38 shows 

the resulting atmospheric photon spectra for E = 400 MeV and E = 3600 MeV, 

c c 

using the electron spectrum of Cline and Brunstein (1966) as the correct 
spectrum at a proton cutoff energy of 400 MeV. We have also shown the result 
of a measurement of the atmospheric photon spectrum by Peterson et al. 

(1967a), at E = 3600 MeV. 

V 

It is apparent that there is a difference of more than an order of magnitude 
between the observed and calculated values of the photon flux at energies above 
30 keV. The difference in the slopes of the observed and predicted spectra may 
actually be more significant as an indication that bremsstrahlung from knock-on 
and Compton-accelerated electrons in the atmosphere does not account for the 
observed atmospheric X-radiation. 


IV. The Photon Component of the Electromagnetic Cascade 
Before considering the electromagnetic cascade itself we will briefly 
summarize the available information on the angular distribution of atmospheric 
X- and gamma radiation near the top of the atmosphere. 

Above 70 MeV the downward photon flux is approximately twice as intense 
as the upward flux (Svensson, 1958; Cline, 1961). A "horizon brightening" was 
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observed (Svensson, 1958), where the horizontal photon flux is approximately 
6 times stronger than the downward flux. X-ray observations {Frost et al. , 
1966; Frost, 1968, private communication) indicate that the upward flux is 
larger than the downward flux by about 50% at energies of approximately 
300 keV. If the contribution of the cosmic diffuse background is subtracted, 
then tho upward flux at approximately 30 keV is still larger than the downward 
flux by 20 - 30%. The same measurements indicate that a horizon brightening 
is not observed at X-ray energies; the photon flux at 90° zenith angle is close 
to that of the upward going photons. 

Measurements by Kniffen (1967) showed that the atmospheric gamma 
spectrum above approximately 60 MeV is comparable to the spectrum which is 
observed from the decay of neutral pions (Perola and Scarsi, 1966), but appears 
to be higher in intensity by a factor of 2. 

Several authors (Anderson, 1961* Peterson, 1963; Brini et al. , 1967) 
have suggested that the large flux of low-energy photons is in equilibrium 
with high-energy photons in the electromagnetic cascade. The great abundance 
of soft photons is explained by the fact that a photon, once it has dropped below 
the pair threshold, will interact by the Compton process, which scatters 
photons many times before their energy is significantly reduced. 

Figure 37 shows the transition curve for photons of energy E = 1 MeV. 

This curve was adapted frbm 1-3 MeV measurements (Peterson et al. , 1967a; 

Peterson, 1967b) with an omnidirectional detector and normalized to the 

/ 2 

1 MeV-photon flux at d = 3. 5 g/cm . 
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The second curve in Figure 37 represents a numerical solution to a one- 
dimensional transport equation for the photon component of the electromagnetic 
cascade: 




The source function for photons from pion decay, Q n , was adapted from the 
calculations of Perola and Scars! (1966). The altitude dependence of this 
function was scaled according to the integral flux of particles with P - 0. 83, 
corresponding to proton kinetic energies above 750 MeV (McDonald and Webber, 
1959). In the calculation we have assumed the spectrum of atmospheric electrons 
to be known as a function of depth z (Beedle and Webber, 1968). The limits 
E . (E ) and E (E ) are defined by the kinematics of Compton interactions. 

Figure 37 shows that a one -dimensional photon transport equation does not 
reproduce the observed transition curve for 1 MeV photons. The observed 
anisotropy of X-radiation near the top of the atmosphere indicates that the con- 
tribution from cascades at all zenith angles must be considered. Also, a one- 
dimensional model clearly breaks down at energies of the order of a few 
MeV and lower. 
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It is evident that a three-dimensional calculation of the electron and photon 

components of the electromagnetic cascade will be required. 

Alsmiller (1968) has followed electromagnetic cascades due to an isotropic 

flux of 100 MeV gamma radiation incident at the top of the atmosphere. The shape 

2 

of the resulting photon spectrum at d ~ 5 g/cm between 0, 01 MeV and a few 
MeV showed qualitative agreement with the observed spectrum (Chapman et al. , 
1968), However, a discrepancy was found between the relative intensities of 
the 511 keV annihilation line in the observed and the predicted photon spectra. 








REFERENCES 


Abell, G. O. , Ann, Rev, Astron. Astrophys. , 3, 1 (196,5). 

Abraham, P. B. , Brunstein, K. A. , and Cline, T. L, , Phys. Rev, 150 , 

1088 (1966), 

Abramowitz, M. and Stegun, I, E. , Handbook of Mathematical Functions 
(National Bureau of Standards, 1964), 

Aitken, D. W, , IEEE Trans , Nucl , Science NS15, 10 (1968). 

Alexander, H. W. , Elements of Mathematical Statistics (Wiley, New York, 

1961). 

Allen, C. W. , Astrophysical Quantities (Athlone Press, London, 1963). 

Aller, L. H. , Handbuch der Physik, S. Fliigge, Ed. (Springer Verlag, 

Berlin, 1958), Vol. 51, p. 324. 

Aller, L. H. , The Abundance of the Elements (Interscience Publishers, 

New York, 1961). 

Alsmiller, R. G. , Working Conference on Upper Atmosphere Radiation and 
Shielding, Oak Ridge (1968), 

Anderson, K. A„ , Phys. Rev. 123 , 1435 (1961), 

B alasubr ahm any an , V. K. , Boldt, E. , and Palmeira, R. A. R. , J. Geophys . 
Res. 72 , » 27 (1967). 

Beacon, R. H. , Am, J, Phys. 2L, 428 (1953). 

Beadle, R, E. , A Balloon Measurement of the Intensity of Low Energy Electrons 
&& a Function of Altitude at Fort Churchill, Canada , M. S. Thesis, Univ. 

•'ft*-! ' 11 1 111 ■■ wiiii — — — n — — mmrnm * — . - - * 

of Minnesota (1966). 

■t 

Beedle, R. E. and Webber, TV. R. , Can. J. Phys., 46, 1014 (1968). 

Bergstrom, I. and Nordling, C, i Alpha-, Beta- and Gamma-Ray Spectroscopy, 
K. Siegbahn, Ed. (North-Holland Publishing Co. , Amsterdam, 1965), 
p. 1523. 

Bethe, H. A. and Ashkin, J. , Experimental Nuclear Physics, E. Segre, Ed. 
(John Wiley & Sons, Inc., New York, 1953), Vol. I, p. 166. 


Blanco, V., Kunkel, W. , Hiltner, W, A. , Lynga, G. , Bradt, II, , Clark, G 
Naranan, S. , Rappaport, S. , and Spada, G. , Ap, J, 152, 1015 (1968a), 


• » 


Blanco, V. , Kunkel, W. , and Hiltner, W. A. , Ap, J. 152, L137 (1968b). 

Bleeker, J. A. M. , Burger, J. J. , Deerenberg, A. J. M. , Scheepmaker, A. , 
Swanenburg, B. N. , and Tanaka, Y. , Ap. J. 147, 391 (1967). 

Bleeker, J. A. M. , Burger, J. J. , Deerenberg, A. J. M. , Scheepmaker, A. , 
Swanenburg, B. N. , Tanaka, Y. , Hayakawa, S, , Maki.no, F. , and Ogawa, 
H. > Can. J. Phys, 46, 461 (1968). 

Bless, R. C., Fischel, D. , and Stecher, T. P. , Ap. J. 151, 1117 (1968). 

Blokhin, M. A. , The Physics of X-rays (State Publishing house of Technical- 
Theoretical Literature, Moscow, 1957). 

Boldt, E., McDonald, F. B. , Riegler, G. , and Serlemitsos, P. , Phys. Rev. 
Letters 17, 447 (1966). 

Boldt. E. and Serlemitsos, P., Suprathermal Proton X-Ray Bremsstrahlung, 
NASA document X-611-68-385 (1968a). 

Boldt, E, A., Desai, U. D. , and Holt, S. S. , 2-20 keV Spectrum of X-Rays 
from the Crab Nebula and the Diffuse Background near Galactic Anticenter , 
NASA document X-611-68-353 (1968b). 

Bowyer, S. , Byram, E. T. , Chubb, T. A. , and Friedman, H. , Nature 201 , 
1307 (1964a). 

Bowyer, S. , Byram, E. T. , Chubb, T. A., and Friedman, H. , Science 146, 
912 (1964b). 

Bowyer, S. , Byram, E. T. , Chubb, T. A. , and Friedman, H. , Science 147, 
394 (1965). 

Bozyan, E. P. , Ap. J . 152 , L155 (1968). 

Bradt, H. , Mayer, W. , Naranan, S. , Rappaport, S. , and Spada, G. , 

Ap. J . 150 , L199 (1967). 


Bradt, H. , Naranan, S. , Rappaport, S. , and Spada, G. , Ap. J. 152, 1005 
(1968a). 

Bradt, H. Naranan, S. , Rappaport, S. , Zwicky, F. , Ogelman, H. , and 
Boldt, E. , Nature 218, 856 (1968b). 

Brini, D. , Ciriegi, U. , Fuligni, F. , Gandolfi, A., and More tti, E., 

J. Geophys, Res, 70, 5460 (1965a). 



124 


Brini, D. , Ciriegi, U. , Fuligni, F. , Gandolfi, A. , and Moretti, E. , Nuovo 
Cimento 38, 130 (1965b). 

Brini, D. , Ciriegi, U. , Fuligni, F. , and Moretti, E. , J. Geophys, Res. 72, 

903 (1967). 

Brunstein, K. A. , J. Geophys. Res. 72, 5510 (1967). 

Burbidge, E. M. , Burbidge, G. R* , Fowler, W. A. , and Hoyle, F. , Rev. Mod. 
Phys. 29, 547 (1957). 

Buselli, G. , Clancy, M. C,, Davison, P. J. N. , Edwards, P. J. , McCracken, 
K. G. , and Thomas, R. M. , Nature 219 , 1124 (1968). 

Cameron, A. G. W. and Mock, M. , Nature 215, 464 (1967). 

Chapman, G. T. , Macklin, J. H. , Gibbons, J. H. , and Nutt, R. , Bull. Am. 
Phys. So c. 13, 695 (1968). 

Chodil, G. , Jopson, R. C. , Mark, H. , Seward, F. D. , and Swift, C. D. , 

Phys. Rev. Letters 15^ 605 (1965). 

Chodil, G. , Mark, H. , Rodrigues, R. , Seward, F. , Swift, C. D. , 

Hiltner, W. A. , Wallerstein, G. , and Mannery, E. J. , Phys. Rev . 

Letters 19, 681 (1967). 

Chodil, G. , Mark, H. , Rodrigues, R. , Seward, F. D. , Swift, C. D. , Turiel, 

I. , Hiltner, W. A. , Wallerstein, G. , and Mannery, E. J. , Ap. J, (to be 
published) (1968). 

Cladis, J. B. and Dess ler, A. J. , J. Geophys. Res. 66, 343 (1961). 

Clark, G. W. , Phys. Rev. Letters 14, 91 (1965). 

Clark, G. W. , Lewin, W. H. G. , and Smith, W. B. , Ap. J. 151, 21 (1968). 

Cline, T. L. , Phys. Rev. Letters 7_, 109 (1961). 

Cline, T. L. and Brunstein, K. A., Proposal to Study the Energy Spectra , 

Time Variations and Anisotropies of Cosmic-Ray, Solar, and Distant 
Magnetos pheric Electrons of Energy 50 keV to 2 MeV with the IMP-H and -J 
and the IMP-I Satellites, NASA document X-611- 66-466 (1966). 


Colgate, S. A. , Can. J. Phys. 46 ? 476 (1968). 

Cooke, B. A., Pounds, K. A., Stewardson, E. A., and Adams, D. J. , 
Ap. J. 150 , L189 (1967). 




125 


Curran, S. C. and Wilson, H. W. , Alpha-, Beta- and Gamma-Ray Spectroscopy, 
K. Siegbahn, Ed. , p. 303 (North-Holland Publishing Co. , Amsterdam, 1965). 

Dicke, R. H. , Peebles, P. J. E. , Roll, P. G. , and Wilkinson, D. T. , Ap. J. 
142,414 (1965). , 

Doong, H. and Noordzy, M. , The Electronics Brain on Board the X-Ray 

Telescope Spectrometer Experiment, NASA document X-611-66-70 (1966). 

ESSA Institute for Environmental Research. Reports of Solar-Geophysical 
Data (U. S. Dept, of Commerce, 1966, 1967). 

Evans, R. D. , The Atomic Nucleus (McGraw-Hill Book Company, Inc. , New 
York, 1955). 

Fano, U. , Phys. Rev. 72, 26 (1947). 

Felten, J. E. and Morrison, P. , Phys. Rev. Letters 10, 453 (1963). 

Felten, J. E. and Morrison, P. , Ap. J. 146, 686 (1966a). 

Felten, J. E. , Gould, R. J. , Stein, W. A. , and Woolf, N. J. , Ap. J. 146, 

955 (1966b). 

Fink, R. W. , Jopson, R. C. , Mark, H. , and Swift, C. D. , Rev. Mod. Phys. 

38, 513 (1966). 

Fischel, D. and Stecher, T. P. , Ap. J. (to be published) (1968). 

Fisher, P. C. , Jordan, W. C. , Meyerott, A. J. , Acton, L. W. , and Roethig, 

D. T. , Ap. J. 151, 1 (1968). 

Freeman, K. C. , Rodgers, A. W. , and Lvnga, G. , Natu re 219, 251 (1968). 

Friedman, H. and Byram, E. T. , Ap. J. 147, 399 (1967a). 

Friedman, H. , Byram, E. T. , and Chubb, T. A. , Scien ce 156, 374 (1967b). 

Frost, K. J. , Rothe, E. D. , and Peterson, L. E. , J. Geophys. Res. 71 , 

4079 (1966). 

Garmire, G. and Kraushaar, W. L., Space Sci. Rev. 4, 123 (1965). 

Ghielmetti, H. S. , Becerra, N. , Godel, A. M. , Heredia, H. , and Roederer, 

G. , Phys. Rev. Letters 12, 388 (1964ai). 

Ghielmetti., H. S. , Becerra, N. , Godel, A. M. , Heredia, H. , Marzulli, 

L* C. and Roederer, J. G. , J. Geophys. Res. 69, 3959 (1964b). 

Giacconi, R. , Gursky, H. , Paolini, F. R. , and Rossi, S. R. , Phys. Rev. 
Letters 9^ 439 (1962). 






126 


Giacconi, R. , Gorenstein, P. , Gursky, H. , and Waters, J. R. , Ap. J. 148 , 
L126 (1967a). 

Giacconi, R. , Gorenstein, P. , Gursky, H. , Usher, P. D. , Waters, J. R. , 
Sandage, A. , Osmer, P. , and Peach, J. V. , Ap. J. 148, L129 (1967b). 

Ginzburg, V, L. and Syrovatskii,, S. I. , The Origin of Cosmic Rays (The 
Macmillan Company, New York, 1964). 


Ginzburg, V. L. and Syro vats Idi, S. I., Soviet Phys.-Usp. 7 ^ 696 (1965). 
Gloeckler, G, and Jokipii, J. R. , Ap. J. 148, L41 (1967). 

Gorenstein, P. , Giacconi, R. , and Gursky, H. , Ap. J. 150 , L85 (T967). 
Gorenstein, P. , Gursky, H. , and Garmire, G. , Ap, J, 153, 885 (1968). 
Gould, R. J. , P hys. Rev. Letters 15, 511 (1965a). 

Gould, R. J., Phys. Rev. Letters 15, 577 (1965b). 

Gould, R. J. , Am. J. Phy s. 35, 376 (1967a). 

Gould, R. J. and Burbidge, G. R. , Handbuch aer Physik, S. Fliigge, Ed. 
(Springer Verlag, Berlin, 1967) Vol. 46-11, p. 265. 


Grader, R. J. , Hill, R. W. , Seward, F. D. , and Toor, A. , Science 152, 
1499 (1966). 

Grader, R. J. , Hill, R. W. , and Seward, F. D. , presented at the 127th 
Meeting of the American Astronomical Society (1968). 

Grodstein, G. W. , X-Ray Attenuation Coefficients from 10 keV to 100 MeV 
(National Bureau of Standards, 1957). 

Gursky, H. , Giacconi, R. , Paolini, F. R. , and Rossi, B. B. , Phys. Rev. 
Letters 11, 530 (1963). ” 

Gursky, H. , Giacconi, R. , Gorenstein, P. , Waters, J. R. , Oda, M. , Bradt, 
H. , Garmire, G, , and Sreekantan, B. V. , Ap. J. 146, 310 (1966). 


Gursky, H. , Gorenstein, P. , and Giacconi, R. , Ap. J. 150, L75 (1967). 

Haro, G. , Nebulae and Interstellar Matter, B. M. Middlehurst and L. H. 
Aller, Ed. , (Univ. of Chicago Press, Chicago, 1968), p. 141. 










f 


i 



I 



127 


Hayakawa, S. , Okuda, H. , Tanaka, Y. , and Yamamoto, Y. , Suppl. Prog. 

Theoret. Phys. (Kyoto) 30, 154 (1964a). 

Hayakawa, S. and Matsuo ka, M. , Suppl. Prog. Theoret. Phys. (Kyoto) 30 , 204 
(1964b). 

Hayakawa, S. , Matsuoka, M. , and Sugimoto, D. , Space Sci, Rev. J5, 109 (1966). 

Haymes, R. C. , Ellis, D. V. , Fishman, G. J. , Glenn, S. W. , and Kurfess, 

J. D. , Ap. J. 151, L125 (1968). 

t ' 

i 

Hill, R. W. , Grader, R. J. , and Seward, F. VV. , Ap. J. (to be published),, 

(19'68). 

Hiltner, W. A. , Mook, D. E., Ludden, D. J. , and Graham, D. , Ap. J. 148, 

L47 (1967). 

Holt, S. S. , Calibration Standards for Photon Spectroscopy between 20-200 keV , 
NASA document X-611-67-433 (1967a). 

Holt, S. S. , and Cline, T. L. , On the Generation of Synchrotron and X-Ray 
Emission from Electrons with Energy below mc^ in Solar Flares , NASA 
document X-611-67-525 (1967b). 

Holt, S. S. , Boldt, E. A., and Serlemitsos, P., Line Emission from X-Ray 
Sour c es, NASA document X-611-68-343 (1968). 

Hoyle, F. , Phys. Rev. Letters 15, 131 (1965). 

Hudson, H. S. , Primbsch, J. H. , and Anderson, K. A. , J. Geophys. Res. 71, 
5665 (1966). 

Jacobson, A. S. , A Search for Gamma-Ray Line Emissi o n from the Crab Nebula, 
Ph. D. Thesis , Univ. of California (1968). 

Johnson, H. M. and Golson, J. G. , Ap. J. 153, 307 (1968). 

Klarmann, J. , Nuovo Cimento 24, 540 (1962). 

Kniffen, D. A., A Search for Gamma-Rays from the Crab Nebula Using a Digitized 
Spark Chamber, P h. D. Thesis , Cath. Univ. of America (1967). 

Kohn, D. , Pinkay, K. , and Wibberenz, G. , Ein Gammaquantenspektrometer fur 
Messungen in der Atmcsphare (Chris tian-Albrechts -Uni versitat Kiel, 1966). 

Kraft, R. P. and Demoulin, M. H. , Ap. J. 150, L183 (196"). 




128 


Lay, H. , Z. Physl k 91, 533 (1934). 

Lewin, W. H. G. , Clark, G. W. , and Smith, W. B. , 

Lewin, W. H. G. , Clark, G. W. , and Smith, W. B. , 

(1968a). 

LG vm, W. H. G. , Clark, G. W. , and Smith, W. B. , 

Lewin, W. H. G. , Clark, G. W. , and Smith, W. B. , 

Manley, O. P. , Ap. J. 144, 1253 (1966). 

Matthews, T. A. and Sandage, A. B. , Ap. J. 138, 30 (1963). 

Maxon, M. S. and Corman, E. G. , Phys. Rev. 163, 156 (1967). 

McDonald, F. B. and Webber, W. R. , Phys. Rev. 115 , 194 (1959). 

Metzger, A. E. , Anderson, E. C. , Van Dilla, M. A. , and Arnold, J. R. , 

Nature 204, 766 (1964). 

Morrison, P. , Nuovo Cimento 7. 858 (1958). 

Nagoya University, Research Institute of Atmospherics, Month. Rept. Solar 
Radio Emission , Nagoya University, Toyokawa (1967). 

Oda, M. , Proc. kit. Conf. Cosmic Rays (The Institute of Physics and 
The Physical Society, London, 1965), p. 68. 

O’Dell, C. R. , Ap. J. 147, 855 (1967). 

Ormes, J. , Proton- and Helium-Nuclei Cosmic-Ray Spectra and Modulation between 
100 and 2000 MeV/Nucleon, Ph. D. Thesis . Univ. of Minnesota (1967). 

Overbeck, J. W. and Tananbaum, H. D. , Ap. J. 153, 899 (1968). 

Pearson, K, , Phil. Mag. 50, 157 (1900). 

Penzias, A. A. and Wilson, R. W. , Ap. J, 146, 666 (1966). 

Perola, G. C. and Scarsi, L. , Nuovo Cimento 46, 718 (1966). 

Peterson, L. E. , J. Geophys. Res. 68, 979 (1963). 

Peterson, L. E. , J)acobson, A. S. , and Pelling, R. M. , Phys. Rev. Letters 16, 

142 (1966a). ' 


Ap. J. 150, L153 (1967). 
Can. J. Phys, 46, 409 

Ap. J. 152, L49 (1968b). 
Ap. J. 152 , L55 (1968c). 


129 


Peterson, L. E, and Jacobson, A. S. , Ap. J, 145, 962 (1966b). 

Peterson, L. E. , Jerde, R. L. , and Jacobson, A. S. , AIAA Journal 5, 1921 
(1967a). 

Peterson, L. E. , Gamma-Ray Production By Cosmic Rays Observed on OSQ-1 
(Univ. of California, San Diego, 1967b). 

Peterson, L. E., Jacobson, A. S. , Peliing, R. M. and Schwartz, D. A., 

Can. J, Phys. 46, 437 (1968). 

Prendergast, K. H. and Burbidge, G. R. , Ap, J. 151, L83 (1968). 

Rossi, B. , High-Energy Particles (Prentice-Hall, Inc. , Englewood Cliffs, 

N. J., 1952). 

Rothenflug, R. , Rocchia, R. , and Koch, L. , Proc. Int. Conf. Cosmic Rays 
(The Institute of Physics and The Physical Society, London, 1965), p. 446, 

Sandage, A. R. , Osmer, P. , Giacconi, R. , Gorenstein, P. , Gursky, H. , 
Waters, J. , Bradt, H. , Garmire, G. , Sreekantan, B. V. , Oda, M. , 
Osawa, K. , and Jugaku, J. , Ap. J. 146, 316 (1966). 

Sartorx, L. and Morrison, P. , Ap. J. 150, 385 (1967a). 

Sartori, L. and Morrison, P. , presented at the Conference on Supernovas, 
held at NASA/Goddard Institute of Space Studies, New York (1967b). 

Schmidt, M. , Ap. J. 141 ? 1 (1965). 

Seward, F. , Chodil, G. , Mark, H. , Swift, C. , and Toor, A. , Ap. J. 150 , 

845 (1967). 

Shea, M. A. , Smart, D. F. , and McCall, J. R. , Can. J. Phys. 46, 1098 
(1968). 

Shklovsky, I. S. , Doklady Akad. Nauk USSR 94, 1020 (1954). 

Shklovsky, I. S. , Ap. J. 148, LI (1967). 

Stein, J. A. and Le win, W. H, G. , J. Geophys. Res. 72^ 383 (1967). 

Svensson, G. , Arkiv f . Fysik 13, 347 (1958). 

Teegarden, B. L. , A Study of Low Energy Galactic Cosmic Rays from 1961 to 
1965,, Ph. D. Thesis, Univ. of Maryland (1967). 


130 


Tucker, W. , Ap. J. 148, 745 (1067). 

Verma, S. D. , A Calculation of the Flux and Energy Spectrum of Secondary 
Electrons at me*. Altitudes in the Atmosphere (University of Chicago, 

EFINS 66-44, 19«6). 

Wallerstein, C. R. , Ap. J. 151, L121 (1968). 

Walt, M. , J. Geophys. Res. 69, 3947 (1964). 

Webber, W. R. , X-Aay Astronomy - 1968 Vintage (Univ. of Adelaide, 1968). 
Westphal, J. A., Sandage, A., and Kristian, J. , Ap. J. 154 , 139 (1968). 

Wharton, C, B. , Plasma Physics, J. E. Drummond, Ed. (McGraw-Hill Co. , 
New York, N. Y. , 1961). 

Williams, J. H. , Phys. Rev. 44 , 146 (1933). 

Wyndham, J. D. , Ap. J. 144 , 459 (1966). 

Zwicky, F. , presented at the Conference on Supernovas, held at NASA/Goddard 
Institute for Space Studies, New York, 1967. 




TT 







Effective 
Relative energy 
Energy line Kr 

Source Tw (keV) intensity (keV) 


Co- 57 267(1. 14.4 

Cd-109 470d. 22.2 

25.0 

87.5 

Sn-119 250d. 24.0 

25.3 

28.5 

Ba-133 7. 2y. 30.6 

31.0 

35.0 

35.8 

81.0 

Eu-155 1. 81y. 39. 5 

42.3 

43.0 

45.0 

48.7 

50.0 

86.5 

Au-195 192d. 65.1 

66. 8 

75.7 

77.9 


Bi-207 304d. 


72.8 
75.0 

84.9 
87.3 


:l } 

.048 


0.12 
0.70 
0.17 
0.05 
0. 55 

0.08 
0.12 
0.26 
0.026 
0.074 
0. 017 
1.0 


22.41 


24.66 


31.45 


43.0 



66.5 


74.5 


Effective Resolution Resolution 
energy spreading spreading 
Csl Kr Csl 

(keV) (keV) (keV) 


22.38 


24.61 


31. 44 


42.85 


67.7 


75.85 















TABLE II. Pulse-height analysis, logical requirements 
and data recording scheme. 


Mode 

Function 

Criteria 

Recorded information 

0 

Csl analysis 

B (20 keV < E < 100 keV) 

PIIA of B in 64 channels 



A (no E > 10 keV) 

indication of C threshold 
fired (4 levels) 

1 

Kr analysis 

A (10 keV< E <30 keV) 
B ( no E >20 keV) 

C ( no E >10 keV) 

PHA of A in 64 channels 

2 

A and B 


scaled number of A events 


bitrate 


in preceding 5 sec interval 
in 64 channels 




scaled number of B events 
in preceding 5 sec interval 
in 64 channels 

3 

B and C 
overload 

B (E > 100 keV) ) 

C ( no E > 10 keV) j 

number of these events 
during 100 msec in 64 
channels 



\ 

number of these events 



C (E > 100 keV) } 

during 100 msec in 64 
channels 

4 

Csl analysis 

same as mode 0, but with r 

nagnetometer indication 

5 

Kr analysis 

same as mode 1, but with magnetometer indication 

6 

A and B 

same as mode 2, but with magnetometer indication 


bitrate 



7 

Timekeeping 


number of occurrences of 
mode 7 (every 12 min) in 
64 channels 
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TABLE HL Geometric factor and full-width angle at half-maximum 

of the effective exposure area. 
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TABLE V, Summary of flight information. 



Eaperimaat 


1 

3 

4 

0 

T 

9 

10 

11 

Flight Date 

6 Doc, G5 

13 Jan. 66 

27 Apr. 66 

2 Sept, G6 

18 Nov. 96 

15 Dec. 86 

18 Dec, 66 

15 Jul, 67 

Launch place 1 

HAFB 

HAFB 

HAFB 

HAFB 

HAFB 

MILDURA 

MILDURA 

PALEST, 

Beginning of ceiling 
coverage (UT) 

1200 

1747 

1220 

1007 

1300 

0130 

0050 

1400 

End of celling 
ooverage (UT) 

1910 

2213 

2230 

2608 

2230 

0338 

0413 

2004 

Residual atmospheric 
depth (g/cm 2 ) 

2. 6-2. 8 

2. 4-2, 6 

2, 7-3. 3 

2.9 

3. 0-3.4 

3. 6-3.1 

3, i-3. a 

. 

3. 2-3. 4 

Rigidity cutoff (GV) 

4. 9-4. 5 

4, 8-3. 2 

4, 7-3,4 

4. 9-5. 1 

4.8-3, 4 

4, 2-4,4 

4,3-4.16 

4. 7-5. 5 

Planetary magnetic 
index K 

P 

0+ to 1- 

lo to 0+ 

Oo to 1- 

6 to 9- 

3- to 4 

6 to 3 + 

1+ to*- 

2+ to 4- 

Detector symbol 

t 

1 

W 

B 3 * 

D 

E 

E 

E 

o 

Scanning method 

VERT+MOD 

VERT+MOD 

VERT+MOD 

VERT+MOD 

VERT 

OH+MOD 

VERT+MOD 

VERT+MOD 

Detector symbol 



A 

W 



G 

G 

2 

Scanning method 



VERT 

ROTATING 



ROTATING 

ROTATING 

4 

Solar flares 

1840 (1-) 

1925 (-B) 
2200 (IN) 

2116 (-F) 

1649 (1-) 
1758 (1-) 
1825 (1-) 
2137 (1-) 
2235 (l-> 
2335 (1-) 

1308 (1-) 
1310 <l-f 
1905 (1-) 
2054 (1-) 
2130 (1-) 
2201 (1-) 
2233 (1-) 

NONE 

0135 (1-) 
0352 (1-) 

1531 (-N) 
1612 (-N) 
1801 (-F) 

Celestial object 
or feature 5 

NGP 

Cyg XR-1 

Cyg xn-i 

NGP 

NGP 

Sco X-l 

Gal, oen. 

Tau XR-1 

UT of meridian 
transit 

1418 

1858 

1230 

2120 

1537 

0116 

0231 

1636 

Celestial object 



Tau XR-1 


Cyg XR-1 




UT of meridian 
transit 



2149 


2156 





1 HAFB = Holloman AFB, Alamogordo, N, M, , 1 - -105. 9°, b = +33. 0° • MILDURA, Australia, 1 = +142, 1° , b = -34. 2°; 
PALESTT. = Palestine, Texas, l = -95.7°, b*31.7°. 

2 VEST = vertical detector, no modulatien, VERT+MOD = vertical detector with w odutaSi o u disk and ring. 

OR+MOD * oriented (aligned for due nerth, fixed zenith angle) ndth rmaJalatkai dtsfc aad ring, ROTATING = rotating 
detector with varying zenith angle. 

Tjstector malfunction, not analyzed, 

* Universal time of maximum phase is shown. The importance of the flare is shown la purewthesls. 

'NGP = North galactic pole; Gal. cea. = Galactic center. 
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TABLE VI. Spectral parameters for the Crab Nebula. Kesults 












TABLE VII. Spectral parameters for the vicinity of the galactic center. 

Kesults of Experiment 10, detector G. 


— — . ■ ■■ ■ 

Procedure 1 

Procedure 2 

Power-law spectrum* 

A = 85.0 

A = 87.0 

J(E) = A E" n 

n = 2.61 

-0.11 

n ~ 2. 60± 0. 92 


Q = 0.11 


Exponential spectrum* 

A = 2.4 

i 

A = 2.32 

J(E) = A exp (-E/n)/E 

+6. 0 

n = 15 '°-4.0 

n = i6 -°-1:o 


Q = 0.043 


Full thermal spectrum* 

A = 0.012 


J( E ) =^exp (-|)k o (|) 

. +13 

n = 40 _ u 


[ 1 + ( 2+0 - 64 1)ht|i 

Q - 0.22 



*J in units of photons/cm -sec-keV, e\ r rgies in keV. 








TABLE vm. X-ray sources in the vicinity of the galactic center. 
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TABLE DX. Photon flux from the Coma cluster of galaxies. Balloon-borne 
measurements at 25 keV, extrapolated to the top of the atmosphere. 


Observation date 

Photon flux 
photons/ cm^-s ec-ke V 

Reference 

25 Nov. 1964 

<9.4 • 10" 4 * 

-4 

Friedman and Byram (1967a) 

6 Dec. 1965 

(140 ± 42) • 10 
-4 

Experiment 1 (Boldt et al. , 1966) 

2 Sept. 1966 

<72 • 10 

-4 

Experiment 6 

18 Nov. 1966 

<50 • 10 

Experiment 7 

13 Feb. 1967 

<15 • 10" 4 

Lewin et al. (1968a) 


*extrapolation from lower-energy rocket measurements assuming a power-law 
of index -2. 














TABLE X. Characteristic line energies and relative intensities of 

Cs, I and Kr. 



Material 

Kq x 

K « 2 

k 3i 

% 

Average 

Cs \ 

energy (keV) 
.relative intensity 

30. 970 
100 

30.623 

50 

34. 984 
30 

35.819 

8 

31.72 

1 1 

energy (keV) 
^relative intensity 

28. 610 
100 

28.315 

50 

32. 292 
30 

33.016 

7.5 

29. 297 

Kr < 

energy (keV) 
^relative intensity 

12. 648 
100 

12.597 

50 

14. 112 
25 

14.313 

3 

12. 867 
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Figure 1 — Essential features of the X-ray telescope 




ENERGY (keV) 



Figure 3 — Energy calibration for the proportional counter of de- 


tector G. Arrows indicate the lower and upper cutoffs. j 

The solid line connects the calibration points * The 

broken line represents a linear approximation to the [ 

curve . [ 
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Figure 5 — Energy calibration for the scintillation crystal 


of detector E. Arrows indicate the upper and 
lower cutoffs 
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Figure 7— Block diagram for the X-ray telescope 
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Figure 9 — Energy dependence of the transfer function 








OBSERVED COUNTING 


G 



LOCAL CUTOFF RIGIDITY (GV) 

Figure 10 — Rigidity dependence of the atmospheric and de- 
tector background for the scintillation crystal, 
22 - 60 keV. Circles: unobstructed detectors. 
Crosses: detectors with occul.tation disk and 
ring. The line connects data points for unob- 
structed detectors 
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4000 H 



UNIVERSAL TIME (HOURS) 


Figure 12 — Observed counting rate in the scintillation 
crystal of detector W (22 - 60 keV) in exper 
iment 4. The exposure to Cyg XR-1 and to 
the Crab Nebula is indicated 
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Figure 14 — Counting rate observed by detector E in experiment 11. 

Top curve: Csl crystal, 22 - 60 keV. Bottom curve: 
proportional counter, 12 - 28 keV. Middle curve: ex- 
posure to a source at the position of the Crab Nebula 
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40 


Figure 15 — Determination of the dominant source declina- 
tion based on the counting rate of detector E in 
experiment 11. Upper curve: Csl crystal data, 
22 - 60 keV. Lower curve: with simulated data 
assuming a source at the position of Tau XR-1 





0 20 30 40 50 60 70 80 90 1 00 

ENERGY ( keV ) 

Figure 16 — Net photon spectrum of the Crab Nebula (ex- 
periment 11), Diamonds: proportional counter 
of detector E. Double-bar crosses: scintil- 
lation crystal of detector E. Single-bar 
crosses: scintillation crystal of detector G 
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Figure 17— Superposition of epochs for the Csl crystal 
data (22 - 60 keV) for experiment 10 (A) and 
experiment 11 (C, E). The calculated expo- 
sures corresponding to A, C, and E are 
shown in B, D, and F, respectively. The 
brackets indicate la-levels 
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Figure 18 — Superposition of epochs for the proportional counter 
data (12-28 keV) for experiment 10 (A) and experi- 
ment 11 (C, E). The calculated exposures corre- 
sponding to A, C, and E are shown in B, D, and F, 
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Figure 19 — Net photon spectrum of Cyg XR-1 (detector A, ex- J 

periment 4). Diamonds: proportional counter. j 

Crosses: Csl crystal. Solid line: measurements 
by Bleeker et al. (1967). Broken line: measure- j 

ments by Peterson et al. (1968) ' u I 
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ENERGY (keV) 


Figure 21 — Net energy spectrum of Sco X-l. Diamonds: pro- 
portional counter results of experiment 9. Crosses 
Csl crystal results of experiment 9. Solid line: 
schematic representation of the results of Peterson 
and Jacobson (1966b), 18 June 1965. Broken line: 
thermal spectrum for T = 5 • 10 7 °K. Dash-dotted 
line: best-fit power- law to the data of Buselli et al. 
(1968), 29 February 1968 
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Figure 21a — The X-ray spectrum of Sco X-l. Crosses: this 
experiment (15 Dec. 1966). Double-bar crosses: 
Lewin et al. (1967) (13 Feb. 1967). Curve 1: 

Chodil et al. (1968) (2 Sept. 1967). Curve 2: 
Gorenstein et al. (1968) (11 Oct. 1966). Curve 3: 
bremsstrahlung from isovelocity plasma. Curve 4: 
bremsstrahlung from isovelocity plasma with 
higher level of ionization. 
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Figure 22 — Trajectory of the axis of detector G on the celestial 
sphere during experiment 10. Six passes from 
galactic north and 6 passes from galactic south are 
drawn. Some relevant sources of low-energy X- 
radiation are shown 
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(22-60 keV) 




Kr 
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Figure 23 — Superposition of epochs for 6 passes from galactic 
north (left side) and from galactic south (right 
side) during experiment 10. The azimuth angles of 
Sgr XR-1 and Sco X-l are indicated. Upper curves: 
Csl crystal, 22 - 60 keV. Lower curves: propor- 
tional counter, 12 - 28 keV. 
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ENERGY (keV) 

Figure 24 — Net photon flux from the vicinity of the galactic center 
(experiment 10, detector G). Diamonds: proportional 
counter. Crosses: scintillation detector: Dotted line: 
best-fit power- law, n = 2. 69. Solid line: best -fit ex- 
ponential law, kT = 16 keV. Broken line: best-fit ex- 
ponential spectrum with kT = 11 keV for the data of 
Lewin et al. (1968b). Dash-dotted line: best-fit power- 
law with n = 2 . 0 for the observation of GX 3+1 by Buselli 
et al. (1968). Circle: expected flux at 4 keV 
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Figure 25 — Determination of the dominant source declination 
based on the counting rate of detector E in experi- 
ment 10. Upper curve: Csl crystal data, 22 - 60 
keV, Lower part: with simulated data assuming a 
source at the position of Sgr XR-1 (solid line and 
the sun) (broken line) 
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RIGHT ASCENSION 

Figure 26 — Map of the sky in the vicinity of the galactic center. 
Broken lines: trajectory of the detector axis of 
telescope G during the first and last rotation of data 
accumulation in experiment 10. Source positions: 
solid lines: Gursky et al. (1967), diamonds: Fried- 
man et al. (1967b), small circles : Bradt et al. (1968a), 
triangles: Fisher et al. (1968), crosses: Buselli et al. 
(1968), large circles: source positions adopted here 
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Figure 28 
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Figure 29 — Superposition of epochs phased in azimuth angle relative 1 

to the direction of the north galactic pole (experiment 6) . j 

Top curve: Csl crystal, 22 - 60 keV. Middle curve: j 

proportional counter , 12-28 keV. Bottom curve: ex- j) 

posure to a point source at the position of the north 

galactic pole !j 

if 




Figure 30 3<r- limiting signals to the photon flux from the vicinity of 

the north galactic pole (experiment 6). Single bars: Csl 
crystal. Double bars: proportional counter 






3 <t- limiting signals to the photon flux from the vicinity of, 
the north galactic pole (experiment 7) . The observed net 
counts of experiment 1 are also given 
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Figure 36 — The atmospheric X-ray spectrum. Curve 1: Peterson 

et al. (1967a), E c = 3600 MeV, d = 3. 5 g/cm 2 . Curve 2: 

calculated for E c = 400 MeV, d = 4 g/cm 2 . Curve 3: 

calculated for E_ = 3600 MeV, d = 4 g/cm^ 
c 
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Figure 37— The atmospheric transition profile for 1 MeV-photons 
Curve 1: Peterson (1967b). Curve 2: solution of a 
one-dimensional photon transport equation (E = 1 MeV 




